Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,729,475:5 +A 5.5% intergenic (‑31/+141) gltW ← / ← rrsG tRNA‑Glu/16S ribosomal RNA of rrnG operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,729,4755.A 5.5% 253.7 / 12.4 73intergenic (‑31/+141)gltW/rrsGtRNA‑Glu/16S ribosomal RNA of rrnG operon
Reads supporting (aligned to +/- strand):  ref base . (31/38);  new base A (2/2);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

GAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGG  >  NC_000913/2729328‑2729620
                                                                                                                                                        |                                                                                                                                                 
gagaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGata                                                                                                                                                        <  3:314695/148‑1 (MQ=32)
gagaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGata                                                                                                                                                        <  7:356320/148‑1 (MQ=32)
 agaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑tt                                                                                                                                                 >  2:364809/1‑149 (MQ=32)
  gaTATTGAGGTCGGCGACTTTCACTCACAAACCATCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑ttt                                                                                                                                                <  6:76076/149‑1 (MQ=17)
  gaTATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑ttt                                                                                                                                                >  1:114339/1‑149 (MQ=32)
    tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTc                                                                                                                                              <  4:270926/149‑1 (MQ=32)
    tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTc                                                                                                                                              <  1:413744/149‑1 (MQ=32)
    tatTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGATGCGTGATA‑‑‑‑‑TTTTc                                                                                                                                              <  3:176996/149‑1 (MQ=17)
     atTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTc                                                                                                                                              <  2:399446/148‑1 (MQ=32)
      ttGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATGCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCgt                                                                                                                                            <  3:130931/149‑1 (MQ=17)
       tGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCgtg                                                                                                                                           >  6:400410/1‑149 (MQ=31)
                           tctgatgaaaagtgaatagcaaggcgtttacgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagctcaccctacacgaaaatatcacgcgttGCGAGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCCTTaa                                                                                                                        <  2:14857‑M1/35‑1 (MQ=255)
                           tctgatgaaaagtgaatagcaaggcgattacgcgttgggagtcaggctgaagagtataaggccgttcgctttctattaatgaaagctcaccctacacgaaaatatgacgcgttGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTaa                                                                                                                        >  1:14857‑M1/114‑148 (MQ=255)
                            actaaCCAGCTAGTGGCGTCCCCTAGGTGATTCGAACCCCTGTTACCGCCGGGTACGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATATCGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTaa                                                                                                                        <  1:425251/144‑1 (MQ=1)
                                                   acgcgttgggagtgaggctgaagagaataaggccgttcgctttctattaatgaaagatcaccctacacgaaaattgcttatcacgcgttGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATtctc                                                                                                >  6:415359‑M1/90‑148 (MQ=255)
                                                                          tgaagagaataaggccgttcgctttctattaatgaaagctcaccctacacgaaaatatcacgcaacGCGTGata                                                                                                                                                        >  4:181171‑M1/67‑74 (MQ=255)
                                                                              gagaataaggccgttcgctttctattaatgaaagctcaccctacacgaaaatatcacgcgttGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTaa                                                                                                                        >  3:398115‑M1/63‑97 (MQ=255)
                                                                              gagaataaggccgttcgctttctattaatgaaagctcaccctacacgaaaatatcacgcgttGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTaa                                                                                                                        <  4:398115‑M1/35‑1 (MQ=255)
                                                                                         gtgtCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCa                                                         <  6:235253/149‑1 (MQ=32)
                                                                                            tCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtt                                                      <  2:205750/149‑1 (MQ=32)
                                                                                            tCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtt                                                      >  7:62822/1‑149 (MQ=32)
                                                                                            tCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACt                                                       <  1:208386/148‑1 (MQ=32)
                                                                                             ccTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACttt                                                     >  1:405634/1‑149 (MQ=32)
                                                                                              cTGGGCCTCTAGACGAAGGGGACATGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                                                    >  3:320087/1‑149 (MQ=17)
                                                                                              cTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                                                    >  1:410888/1‑149 (MQ=32)
                                                                                               tGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCa                                                         <  7:332479/143‑1 (MQ=32)
                                                                                               tGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCa                                                         >  8:332479/1‑143 (MQ=32)
                                                                                                gggCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtca                                                  >  5:457447/1‑149 (MQ=34)
                                                                                                 ggCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc                                                   <  1:270591/147‑1 (MQ=33)
                                                                                                 ggCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc                                                   <  8:180762/147‑1 (MQ=33)
                                                                                                  gCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatc                                                <  2:343300/149‑1 (MQ=35)
                                                                                                   ccTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatca                                               <  2:95535/149‑1 (MQ=37)
                                                                                                     tctAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGa                                             >  8:438994/1‑149 (MQ=37)
                                                                                                         gACGCAGGGGACACGAAAATTGCTTATCTCGCGTTGCGTGATA‑‑‑‑‑TTTTCGGGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat                                                 <  8:221815/141‑1 (MQ=255)
                                                                                                         gACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                                                    >  2:174668/1‑138 (MQ=33)
                                                                                                         gACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                                                    <  1:174668/138‑1 (MQ=33)
                                                                                                         gACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTGTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCCCTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat                                                 >  7:221815/1‑141 (MQ=25)
                                                                                                          aCGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATc                                        <  5:7558/149‑1 (MQ=37)
                                                                                                            gAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtg                                      >  1:40224/1‑149 (MQ=37)
                                                                                                             aaGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgt                                     <  8:108370/149‑1 (MQ=37)
                                                                                                               ggggACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgt                                   <  2:40224/149‑1 (MQ=37)
                                                                                                               ggggACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgt                                   <  6:457447/149‑1 (MQ=37)
                                                                                                                  gACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAg                                >  6:390303/1‑149 (MQ=37)
                                                                                                                     acGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAAt                                         >  5:47971/1‑137 (MQ=37)
                                                                                                                     acGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAAt                                         <  6:47971/137‑1 (MQ=37)
                                                                                                                     acGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCAc                             >  8:96711/1‑149 (MQ=37)
                                                                                                                       tacacgaaaatatcacgcaacGCGTGATAAGCAATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgtgt                                   <  1:434900‑M1/125‑1 (MQ=255)
                                                                                                                       gAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTg                           <  7:438994/149‑1 (MQ=38)
                                                                                                                        aaaaTTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCCTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGc                          <  5:230965/149‑1 (MQ=38)
                                                                                                                         cacgaaaatatcacgcaacGCGTGATAAGCAATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGa                                 >  2:434900‑M1/20‑146 (MQ=255)
                                                                                                                          aaTTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCa                         >  5:92750/1‑148 (MQ=39)
                                                                                                                           aTTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCaaa                       <  5:390303/149‑1 (MQ=39)
                                                                                                                            ttGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCATCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAg                      >  2:66765/1‑149 (MQ=255)
                                                                                                                            gaaaatatcacgcaacGCGTGATAAGCAATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTTTTCActttc                                                    <  4:375386‑M1/108‑2 (MQ=255)
                                                                                                                            gaaaatatcacgcaacGCGTGATAAGCAATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTTTTCActttc                                                    >  3:375386‑M1/17‑123 (MQ=255)
                                                                                                                             tGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACTCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAg                      >  4:41623/1‑148 (MQ=255)
                                                                                                                               aatatcacgcaacGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc                                                   <  5:245074‑M1/104‑1 (MQ=255)
                                                                                                                               aatatcacgcaacGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc                                                   >  6:245074‑M1/14‑117 (MQ=255)
                                                                                                                                ttATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgt                                     <  1:140954/130‑1 (MQ=38)
                                                                                                                                ttATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtgt                                     >  2:140954/1‑130 (MQ=38)
                                                                                                                                  aTCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGct                >  1:6795/1‑149 (MQ=39)
                                                                                                                                       gcgTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttctt            >  7:415878/1‑148 (MQ=39)
                                                                                                                                        cgTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTta          <  2:6795/149‑1 (MQ=39)
                                                                                                                                         gTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTCATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaa         <  8:296446/149‑1 (MQ=255)
                                                                                                                                         gTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCttt           >  2:41678/1‑147 (MQ=39)
                                                                                                                                         gTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaa         <  3:123936/149‑1 (MQ=39)
                                                                                                                                          ttGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaag        <  8:415878/149‑1 (MQ=39)
                                                                                                                                            GCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaagg       <  5:376942‑M1/148‑1 (MQ=255)
                                                                                                                                               TGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaaggtaag   >  6:92750‑M1/1‑149 (MQ=255)
                                                                                                                                                GATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTaaa                                                                       >  5:344969‑M1/1‑80 (MQ=255)
                                                                                                                                                GATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTaaa                                                                       <  6:344969‑M1/80‑1 (MQ=255)
                                                                                                                                                GATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttctt            <  5:348986‑M1/139‑1 (MQ=255)
                                                                                                                                                GATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttctt            >  6:348986‑M1/1‑139 (MQ=255)
                                                                                                                                                 ATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAACGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTCCAAAGTACGCTTCTTtaaggtaagg  >  1:61717‑M1/1‑148 (MQ=255)
                                                                                                                                                 ATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttc              >  5:340098‑M1/1‑136 (MQ=255)
                                                                                                                                                 ATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGcttc              <  6:340098‑M1/136‑1 (MQ=255)
                                                                                                                                                 ATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTtaaggtaag   <  1:228107‑M1/147‑1 (MQ=255)
                                                                                                                                                   A‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCAc                             >  2:239540‑M1/1‑119 (MQ=255)
                                                                                                                                                   A‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCAc                             <  1:239540‑M1/119‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                                 
GAGATATTGAGGTCGGCGACTTTCACTCACAAACCAGCAAGTGGCGTCCCCTAGGGGATTCGAACCCCTGTTACCGCCGTGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACACGAAAATTGCTTATCACGCGTTGCGTGATA‑‑‑‑‑TTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTGTGTGAGCACTGCAAAGTACGCTTCTTTAAGGTAAGG  >  NC_000913/2729328‑2729620

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