Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,168,350 A→T 8.4% intergenic (+150/‑22) rrsB → / → gltT 16S ribosomal RNA of rrnB operon/tRNA‑Glu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,168,3500AT 8.4% 187.1 / 6.4 71intergenic (+150/‑22)rrsB/gltT16S ribosomal RNA of rrnB operon/tRNA‑Glu
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base T (3/3);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

CTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTC  >  NC_000913/4168202‑4168495
                                                                                                                                                         |                                                                                                                                                 
cttaAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCaa                                                                                                                                                   >  7:354391/1‑149 (MQ=39)
cttaAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCaa                                                                                                                                                   >  6:416725/1‑149 (MQ=39)
 ttaAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAAc                                                                                                                                                  >  4:351825/1‑149 (MQ=39)
 ttaAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAAc                                                                                                                                                  >  2:465154/1‑149 (MQ=39)
   aaaGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgc                                                                                                                                                <  3:72345‑M2/149‑3 (MQ=255)
   aaaGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACg                                                                                                                                                 >  6:7725‑M2/1‑147 (MQ=255)
     agaagCGTACTTTGTAGTGCTCACACAGATTGTTTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgt                                                                                                                                              <  3:234344‑M2/149‑5 (MQ=255)
     agaagCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAGGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgt                                                                                                                                              <  6:74627‑M2/149‑5 (MQ=255)
     agaagCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAGGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgt                                                                                                                                              >  5:74627‑M2/1‑145 (MQ=255)
     agaagCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCCACCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgt                                                                                                                                              <  7:412586‑M2/149‑5 (MQ=255)
     agaagCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgt                                                                                                                                              <  4:175548‑M2/149‑5 (MQ=255)
       aagCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCGTTGCGTGa                                                                                                                                            <  3:255357/149‑1 (MQ=255)
           gTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAg                                                                                                                                        >  1:321382/1‑149 (MQ=39)
           gTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAg                                                                                                                                        <  5:131271/149‑1 (MQ=39)
           gTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAg                                                                                                                                        >  4:75757/1‑149 (MQ=39)
             aCTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGc                                                                                                                                       >  5:20390/1‑148 (MQ=39)
              cTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCaa                                                                                                                                     <  2:321382/149‑1 (MQ=39)
                  gTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAAtttt                                                                                                                                 >  5:25870/1‑149 (MQ=38)
                  gTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAAttt                                                                                                                                  <  2:5587/148‑1 (MQ=38)
                   tAGTGCTCACACAGATTGTCTGATAGGAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTc                                                                                                                                <  4:388794/149‑1 (MQ=38)
                   tAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTc                                                                                                                                <  8:206701/149‑1 (MQ=38)
                         tCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTcc                                                                                                                          >  6:263735/1‑149 (MQ=37)
                          cacacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTccc                                                                                                                         >  1:175845/1‑149 (MQ=37)
                           acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTcccc                                                                                                                        >  6:42152/1‑149 (MQ=37)
                           acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTcccc                                                                                                                        >  1:292581/1‑149 (MQ=37)
                           acacaGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTcc                                                                                                                          <  6:217224/147‑1 (MQ=37)
                               aGATTGTCTGATAGAAAGTGAGAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCg                                                                                                                    <  7:80556/149‑1 (MQ=37)
                                gATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGt                                                                                                                   <  1:465154/149‑1 (MQ=37)
                                 aTTGTCTGATATAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTTAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACTCGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTc                                                                                                                  <  3:80370/149‑1 (MQ=255)
                                     tttGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTaga                                                                                                              <  4:234344/147‑1 (MQ=37)
                                     tCTGATGAAAAGTGAATAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCGTTgcgagatattttcgtgtagggtgagctttccttaa                                                                                                               <  2:14857‑M2/148‑36 (MQ=255)
                                     tCTGATGAAAAGTGAATAGCAAGGCGATTACGCGTTGGGAGTCAGGCTGAAGAGTATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATGACGCGTTgcgtgatattttcgtgtagggtgagctttcattaa                                                                                                               >  1:14857‑M2/1‑113 (MQ=255)
                                     tCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTaga                                                                                                              <  5:256833/149‑1 (MQ=37)
                                       tGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTcc                                                                                                                          <  4:318810/135‑1 (MQ=37)
                                       tGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTcc                                                                                                                          >  3:318810/1‑135 (MQ=37)
                                       tGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAgg                                                                                                            >  8:425397/1‑149 (MQ=37)
                                        gATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGc                                                                                                           >  1:185276/1‑149 (MQ=35)
                                         aTAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTagag                                                                                                             <  8:354391/146‑1 (MQ=35)
                                          tAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGcc                                                                                                          >  4:318258/1‑148 (MQ=34)
                                           aGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgt                                                                                                                                              >  2:55316‑M2/1‑107 (MQ=255)
                                           aGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgt                                                                                                                                              <  1:55316‑M2/111‑5 (MQ=255)
                                            gAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCa                                                                                                        <  6:270338/148‑1 (MQ=32)
                                            gAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAg                                                                                                       <  5:416725/149‑1 (MQ=32)
                                              aaGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGa                                                                                                     <  3:351825/149‑1 (MQ=32)
                                              aaGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGa                                                                                                     <  2:185276/149‑1 (MQ=32)
                                                 tGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCaa                                                                                                                                                   >  3:327910/1‑100 (MQ=33)
                                                 tGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCaa                                                                                                                                                   <  4:327910/100‑1 (MQ=33)
                                                  gAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACAcc                                                                                                 >  1:213992/1‑149 (MQ=32)
                                                   aaaaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACAcc                                                                                                 <  6:25870/148‑1 (MQ=32)
                                                   aaaaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCg                                                                                                <  7:261431/149‑1 (MQ=32)
                                                    aaaGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCTATTTTCGTGTCCCCTTCGTCTTGAGGCCCAGCACAcccc                                                                                               >  2:438817/1‑147 (MQ=9)
                                                          agtcGTTTACGCGTTGGGAGTGAGGCTGAAGATAATAAGGCCGTTCGCTTTCTATTAATTAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTc                                                                                         <  5:96449/146‑1 (MQ=0)
                                                                  aCGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGATCACCCTACACGAAAATTGCTTATCACGCGTTgcgtgatattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctc                                                                                       >  6:415359‑M2/1‑89 (MQ=255)
                                                                    gcgTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgtgat                                                                                                                                           <  2:49124‑M2/89‑8 (MQ=255)
                                                                    gcgTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgtgat                                                                                                                                           >  1:49124‑M2/1‑82 (MQ=255)
                                                                                    tGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACgcgtgata                                                                                                                                          >  4:181171‑M2/1‑66 (MQ=255)
                                                                                        gagaATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCGTTgcgtgatattttcgtgtagggtgagctttcattaa                                                                                                               >  3:398115‑M2/1‑62 (MQ=255)
                                                                                        gagaATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCGTTgcgtgatattttcgtgtagggtgagctttcattaa                                                                                                               <  4:398115‑M2/97‑36 (MQ=255)
                                                                                                             ctcttcaTGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCACCGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGttt                                       <  8:412586/142‑1 (MQ=0)
                                                                                                                                 tACACGAAAAT‑‑‑‑‑ATCACGCAACgcgtgataagcaattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgctattcacttttcatcagacaatctgtgt                     <  1:434900‑M2/146‑126 (MQ=255)
                                                                                                                                   cacGAAAAT‑‑‑‑‑ATCACGCAACgcgtgataagcaattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgctattcacttttcatcagacaatctgtgtga                   >  2:434900‑M2/1‑19 (MQ=255)
                                                                                                                                    acGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTaa               >  7:54138/1‑149 (MQ=32)
                                                                                                                                      gAAAAT‑‑‑‑‑ATCACGCAACgcgtgataagcaattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgcttttcactttc                                      >  3:375386‑M2/1‑16 (MQ=255)
                                                                                                                                      gAAAAT‑‑‑‑‑ATCACGCAACgcgtgataagcaattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgcttttcactttc                                      <  4:375386‑M2/124‑109 (MQ=255)
                                                                                                                                      gAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAata             >  5:375500/1‑149 (MQ=32)
                                                                                                                                      gAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAata             >  7:164199/1‑149 (MQ=32)
                                                                                                                                      gAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAata             >  6:357819/1‑149 (MQ=32)
                                                                                                                                      gAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAata             >  5:399049/1‑149 (MQ=32)
                                                                                                                                         aaT‑‑‑‑‑ATCACGCAACgcgtgatattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgctattcacttttc                                          >  6:245074‑M2/1‑13 (MQ=255)
                                                                                                                                         aaT‑‑‑‑‑ATCACGCAACgcgtgatattttcgtgtagggtgagctttcattaatagaaagcgaacggccttattctcttcagcctcactcccaacgcgtaaacgccttgctattcacttttc                                          <  5:245074‑M2/117‑105 (MQ=255)
                                                                                                                                          at‑‑‑‑‑atCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATAtctc         <  5:357819/149‑1 (MQ=32)
                                                                                                                                                      gCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACt   >  4:365721/1‑148 (MQ=32)
                                                                                                                                                      gCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACt   >  5:412010/1‑148 (MQ=32)
                                                                                                                                                      gCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTc  <  7:180116/149‑1 (MQ=33)
                                                                                                                                                      gCAACGCGTGATAAGCAATTTTAGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTc  <  3:318258/149‑1 (MQ=18)
                                                                                                                                                        aaCGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCaaaa     >  5:187765/1‑144 (MQ=32)
                                                                                                                                                        aaCGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCaaaa     <  6:187765/144‑1 (MQ=32)
                                                                                                                                                         |                                                                                                                                                 
CTTAAAGAAGCGTACTTTGTAGTGCTCACACAGATTGTCTGATAGAAAGTGAAAAGCAAGGCGTTTACGCGTTGGGAGTGAGGCTGAAGAGAATAAGGCCGTTCGCTTTCTATTAATGAAAGCTCACCCTACACGAAAAT‑‑‑‑‑ATCACGCAACGCGTGATAAGCAATTTTCGTGTCCCCTTCGTCTAGAGGCCCAGGACACCGCCCTTTCACGGCGGTAACAGGGGTTCGAATCCCCTAGGGGACGCCACTTGCTGGTTTGTGAGTGAAAGTCGCCGACCTTAATATCTCAAAACTC  >  NC_000913/4168202‑4168495

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: