Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,409,402 Δ1 bp 15.1% coding (128/732 nt) rlmB → 23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,409,4020T.15.1% 413.4 / 68.3 126coding (128/732 nt)rlmB23S rRNA mG2251 2'‑O‑ribose methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (55/52);  new base . (10/9);  total (65/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 131 total aligned reads displayed.

GTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTC  >  NC_000913/4409256‑4409550
                                                                                                                                                  |                                                                                                                                                     
gTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTg                                                                                                                                                     >  4:9016/1‑149 (MQ=255)
gTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTg                                                                                                                                                     >  1:375629/1‑149 (MQ=255)
   aacaacGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGt                                                                                                                                                  <  2:478980/149‑1 (MQ=255)
    acaacGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGtc                                                                                                                                                 >  4:1463843/1‑149 (MQ=255)
     caacGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGtct                                                                                                                                                >  4:1004997/1‑149 (MQ=255)
       acGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑a                                                                                                                                             <  2:375629/149‑1 (MQ=255)
        cGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑Ag                                                                                                                                            <  4:627831/149‑1 (MQ=255)
             aCATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCtgagactc                                                                                                                                                <  1:803469/141‑8 (MQ=255)
               aTTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTgg                                                                                                                                     <  3:724079/149‑1 (MQ=255)
               aTTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTgg                                                                                                                                     <  2:419594/149‑1 (MQ=255)
               aTTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTgg                                                                                                                                     >  1:419594/1‑149 (MQ=255)
                  aaTGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑Agg                                                                                                                                           <  3:1306235/140‑1 (MQ=255)
                  aaTGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑Agg                                                                                                                                           >  4:1306235/1‑140 (MQ=255)
                         gAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGa                                                                                                                                                    <  2:1252082/125‑1 (MQ=255)
                         gAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGa                                                                                                                                                    >  3:240001/1‑125 (MQ=255)
                         gAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGa                                                                                                                                                    <  4:240001/125‑1 (MQ=255)
                         gAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGa                                                                                                                                                    >  1:1252082/1‑125 (MQ=255)
                         gAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGtt                                                                                                                           <  2:1019127/149‑1 (MQ=255)
                         gAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGtt                                                                                                                           >  2:722367/1‑149 (MQ=255)
                           aaTGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGa                                                                                                                                                    >  1:1025318/1‑123 (MQ=255)
                           aaTGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGa                                                                                                                                                    <  2:1025318/123‑1 (MQ=255)
                           aaTGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTgg                                                                                                                         <  3:1463843/149‑1 (MQ=255)
                               aTTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCC‑TGAGACTC‑Aagg                                                                                                                                          <  2:992851/127‑1 (MQ=255)
                               aTTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCC‑TGAGACTC‑Aagg                                                                                                                                          >  1:992851/1‑127 (MQ=255)
                                tttACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTg                                                                                                                                      >  4:370327/1‑131 (MQ=255)
                                tttACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTg                                                                                                                                      >  4:1024187/1‑131 (MQ=255)
                                tttACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTg                                                                                                                                      <  3:370327/131‑1 (MQ=255)
                                 ttACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGtt                                                                                                                           >  4:27654/1‑141 (MQ=255)
                                        aTCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAata                                                                                                            >  4:458773/1‑149 (MQ=255)
                                        aTCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCC‑TGAGACTCAAGGGCGt                                                                                                                                       >  3:481548/1‑122 (MQ=255)
                                        aTCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCC‑TGAGACTCAAGGGCGt                                                                                                                                       <  4:736843/122‑1 (MQ=255)
                                            aCGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGt                                                                                                                                    >  2:212758/1‑121 (MQ=255)
                                                  tGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACgag                                                                                                  <  3:972071/149‑1 (MQ=255)
                                                      ggCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGaaaa                                                                                              <  3:1004997/149‑1 (MQ=255)
                                                         ccTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTa                                                                                                                                  >  3:87669/1‑110 (MQ=255)
                                                         ccTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTa                                                                                                                                  <  4:87669/110‑1 (MQ=255)
                                                                     cgcCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGc                                                                               >  2:1319521/1‑149 (MQ=255)
                                                                         ccTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGt                                                                                                                                    <  2:909117/92‑1 (MQ=255)
                                                                         ccTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGt                                                                                                                                    >  1:909117/1‑92 (MQ=255)
                                                                          cTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAg                                                                          >  1:1024888/1‑149 (MQ=255)
                                                                             aaCGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTg                                                                                                                                      >  3:686096/1‑86 (MQ=255)
                                                                             aaCGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCaa                                                                                                              >  1:60109/1‑110 (MQ=255)
                                                                             aaCGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGc                                                                       >  3:1132577/1‑149 (MQ=255)
                                                                              aCGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCa                                                                      <  1:1066104/149‑1 (MQ=255)
                                                                                 ttttCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg                                                                         >  4:1100343/1‑143 (MQ=255)
                                                                                 ttttCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAgg                                                                         <  3:1100343/143‑1 (MQ=255)
                                                                                 ttttCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTa                                                                   <  2:1024888/149‑1 (MQ=255)
                                                                                   ttCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATc                                                                 <  1:722367/149‑1 (MQ=255)
                                                                                     cAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGc                                                               <  1:1319521/149‑1 (MQ=255)
                                                                                         aaGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCAt                                                                             >  4:442364/1‑131 (MQ=255)
                                                                                         aaGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCAt                                                                             <  3:442364/131‑1 (MQ=255)
                                                                                           gTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGt                                                         >  4:840095/1‑149 (MQ=255)
                                                                                             cTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCC‑TGAGACTCAAGGGCGt                                                                                                                                       <  3:718704/69‑1 (MQ=37)
                                                                                             cTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCC‑TGAGACTCAAGGGCGt                                                                                                                                       >  4:718704/1‑69 (MQ=37)
                                                                                              tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTg                                                                                                                                      <  3:927313/69‑1 (MQ=255)
                                                                                              tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTg                                                                                                                                      >  4:927313/1‑69 (MQ=255)
                                                                                              tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                      >  2:907897/1‑149 (MQ=255)
                                                                                              tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCC‑TGAGACTCAAGGGCGt                                                                                                                                       >  2:697076/1‑68 (MQ=37)
                                                                                              tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCC‑TGAGACTCAAGGGCGTga                                                                                                                                     <  1:1448973/70‑2 (MQ=25)
                                                                                              tttATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCC‑TGAGACTCAAGGGCGTga                                                                                                                                     >  2:1448973/1‑69 (MQ=25)
                                                                                                       aaaGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACg                                             >  1:193265/1‑149 (MQ=255)
                                                                                                               tGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATc                                     >  3:1445373/1‑149 (MQ=255)
                                                                                                                ggaGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCa                                    >  1:311523/3‑149 (MQ=255)
                                                                                                                gAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGc                                                                                                                 >  2:639643/1‑72 (MQ=255)
                                                                                                                gAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGc                                                                                                                 <  1:639643/72‑1 (MQ=255)
                                                                                                                   gATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGa                                                                                          >  1:949658/1‑92 (MQ=255)
                                                                                                                   gATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGa                                 >  3:819947/1‑149 (MQ=255)
                                                                                                                    aTAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGaa                                >  3:223788/1‑149 (MQ=255)
                                                                                                                     tAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGc                                                                                            <  1:179229/88‑1 (MQ=255)
                                                                                                                     tAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGc                                                                                            >  2:179229/1‑88 (MQ=255)
                                                                                                                     tAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGaaa                               >  1:341983/1‑149 (MQ=255)
                                                                                                                      aaaCGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                      >  4:1185130/1‑125 (MQ=255)
                                                                                                                      aaaCGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGaa                                                      <  3:1185130/125‑1 (MQ=255)
                                                                                                                          gTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATc                                                                            >  2:1447136/1‑99 (MQ=255)
                                                                                                                          gTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGAt                          <  3:840095/149‑1 (MQ=255)
                                                                                                                          gTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGAt                          >  2:8604/1‑149 (MQ=255)
                                                                                                                          gTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGAt                          <  1:8604/149‑1 (MQ=255)
                                                                                                                               ttACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGAt                          <  2:1106730/144‑1 (MQ=255)
                                                                                                                                  ccGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTc                                           <  3:87147/124‑1 (MQ=255)
                                                                                                                                   cGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAAc                             >  3:29464/1‑137 (MQ=255)
                                                                                                                                    gCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATc                >  1:110005/1‑149 (MQ=255)
                                                                                                                                    gCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATc                >  1:102050/1‑149 (MQ=255)
                                                                                                                                       gATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGa             <  4:223788/149‑1 (MQ=255)
                                                                                                                                           cACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtc           >  4:624815/1‑147 (MQ=255)
                                                                                                                                           cACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGt                                                                                 >  2:1057202/1‑77 (MQ=37)
                                                                                                                                           cACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGt                                                                                 <  1:1057202/77‑1 (MQ=37)
                                                                                                                                           cACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgc         >  2:282571/1‑149 (MQ=255)
                                                                                                                                           cACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTCGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt                 <  4:483540/141‑1 (MQ=255)
                                                                                                                                           cACGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTCGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt                 >  3:483540/1‑141 (MQ=255)
                                                                                                                                            aCGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATc                                                                            <  3:1338863/81‑1 (MQ=255)
                                                                                                                                            aCGCCC‑TGAGACTCAAGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGAt                 >  1:1372758/1‑140 (MQ=255)
                                                                                                                                              gCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGAtcgcttc      <  2:447891/149‑1 (MQ=255)
                                                                                                                                               cccTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAg                                                 >  4:109738/1‑105 (MQ=255)
                                                                                                                                               cccTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAg                                                 <  3:1005717/105‑1 (MQ=255)
                                                                                                                                               cccTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAg                                                 <  3:109738/105‑1 (MQ=255)
                                                                                                                                                ccTTGAGTCTC‑AGGGCGTGTTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAgg                                                >  1:1006250/1‑105 (MQ=255)
                                                                                                                                                 cTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTc                                           >  1:63785/1‑109 (MQ=255)
                                                                                                                                                 cTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTc                                           <  2:63785/109‑1 (MQ=255)
                                                                                                                                                  ttGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATc                                     >  4:1141847/1‑114 (MQ=255)
                                                                                                                                                  ttGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTc  >  3:556616/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTAAACAACGAGTACATTAATGAGCGAAATGATTTACGGCATCCACGCAGTGCAGGCCCTGCTGGAGCGCGCCCCTGAACGTTTTCAGGAAGTCTTTATTTTAAAAGGCCGTGAAGATAAACGTCTGTTACCGCTGATTCACGCCCTTGAGTCTC‑AGGGCGTGGTTATCCAGTTGGCAAACCGCCAATATCTCGACGAGAAAAGCGACGGTGCCGTGCATCAGGGCATTATCGCCCGCGTGAAGCCAGGACGTCAGTATCAGGAAAACGATCTGCCGGATCTGATCGCTTCGCTC  >  NC_000913/4409256‑4409550

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑