Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,113,177 A→G intergenic (‑116/+257) rfbB ← / ← galF dTDP‑glucose 4,6‑dehydratase 1/UTP:glucose‑1‑phosphate uridylyltransferase, low activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,113,1770AG97.6% 140.3 / NA 41intergenic (‑116/+257)rfbB/galFdTDP‑glucose 4,6‑dehydratase 1/UTP:glucose‑1‑phosphate uridylyltransferase, low activity
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (18/22);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CAAGTGATTCCCGGTTTCCGGCGTACGTTAATTTATCGACATTAACAACACTATCCTGCGTATTATTTATAATGTGACGAACTACAGCTGAACCAATAAATCCTGCGCCACCAGTAACAAGTATTTTCACTTAATTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑ATGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCT  >  NC_000913/2112932‑2113403
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
cttggacattacatactcctgtcaggtgtgatcaccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTgc‑gtgcat                                                                                                                                                                                                                               >  2:275014‑M2/136‑251 (MQ=255)
  tggacattacatactcctgtcaggtgtgatccccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggttttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTgc‑gtgcat                                                                                                                                                                                                                               <  2:13993‑M2/116‑1 (MQ=255)
    gacattacatactcctgtcaggtgtgatcaccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCt                                                                                                                                                                                                                           <  1:493341‑M2/120‑1 (MQ=255)
          acatactcctgtcaggtgtgatcaccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTatgatg                                                                                                                                                                                                                     >  1:13016‑M2/126‑251 (MQ=255)
                    gtcaggtgtgatcaccgcgctgaatgctcgcggcctttaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgtttacggctttatTTATTCCATATTACTTCGGAGCATGCTTTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCATGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTTCATGCCTATGATGTTACGAGCGt                                                                                                                                                                                                           <  2:45499‑M2/136‑1 (MQ=255)
                     tcaggtgtgatcaccgcgctgcatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGTGGTATTAATGCACGCTACCGACCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGtt                                                                                                                                                                                                          >  2:255133‑M2/115‑251 (MQ=255)
                        ggcgtgatcaccgcgctgaatgctcgcggcctggaatgttcgcaaaaaagtatacaaagtacttagtaaagtataaatagcttagccatgcttatattcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTAGTAAGGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTAcc                                                                                                                                                                                                       <  2:253232‑M2/140‑1 (MQ=255)
                                caccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTcgc                                                                                                                                                                                                <  1:236531‑M2/147‑1 (MQ=255)
                                    gcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTTCTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTaa                                                                                                                                                                                           <  1:572010‑M2/152‑1 (MQ=255)
                                                         gaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTa                                                                                                                                                                      <  1:259800‑M2/173‑1 (MQ=255)
                                                           atgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAAt                                                                                                                                                                    <  1:275256‑M2/175‑1 (MQ=255)
                                                               tcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGt                                                                                                                                                                 <  1:176023‑M2/178‑1 (MQ=255)
                                                                 gcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCtt                                                                                                                                                              <  1:275014‑M2/181‑1 (MQ=255)
                                                                       aagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAAt                                                                                                                                                        >  1:305440‑M2/65‑251 (MQ=255)
                                                                           atacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCaatcagc                                                                                                                                                    >  2:423167‑M2/61‑248 (MQ=255)
                                                                                       tattaaactataaatagcttatccctgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGCGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCACTCCGCTCTGGAAACa                                                                                                                                          >  2:247209‑M2/49‑249 (MQ=255)
                                                                                        attaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAgg                                                                                                                                       >  2:210268‑M2/48‑251 (MQ=255)
                                                                                        attaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAgg                                                                                                                                       <  2:148192‑M2/204‑1 (MQ=255)
                                                                                           aaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCACTCCGCTCTGGAACCAAGGaa                                                                                                                                     >  2:416809‑M2/45‑250 (MQ=255)
                                                                                             cgtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTAATTCAGAGCATGCTGAGAACTAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCCCCCCTGGCTTTACAGCTACCAGTGCACTGC‑GTGCATGCCTATGATGTGACGAGCGGTACCCACTCGCGCTAACCCCGAACAACTCACCCGTTCATTGGCTTACCAATCCGCTCTGGACacaaggcata                                                                                                                                   >  2:348805‑M2/43‑246 (MQ=255)
                                                                                                   aatagcttatccatgcttattggcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCATTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGa                                                                                                                            <  2:62325‑M2/215‑1 (MQ=255)
                                                                                                         ttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACtt                                                                                                                      >  1:533935‑M2/31‑251 (MQ=255)
                                                                                                         ttatccatgattatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACtt                                                                                                                      >  1:66219‑M2/31‑251 (MQ=255)
                                                                                                          tatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGTTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACtt                                                                                                                      <  2:293006‑M2/221‑1 (MQ=255)
                                                                                                                    tatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAAcc                                                                                                           >  2:374560‑M2/20‑251 (MQ=255)
                                                                                                                       atgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCt                                                                                                         <  1:210268‑M2/234‑1 (MQ=255)
                                                                                                                       atgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCt                                                                                                         >  1:580467‑M2/17‑250 (MQ=255)
                                                                                                                         gcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATACAGGTATTAATGCACTCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTCCCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGATTAAAACCCTAc                                                                                                       <  1:47624‑M2/236‑1 (MQ=255)
                                                                                                                                     aTTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCgtt                                                                                           <  1:276190/250‑1 (MQ=255)
                                                                                                                                                 ttACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATgg                                                                              >  2:179468/1‑251 (MQ=255)
                                                                                                                                                   aCTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGtt                                                                            >  1:35843/1‑251 (MQ=255)
                                                                                                                                                       cAGAGCATGCTGTGAAATACGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACAGCCACTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTtat                                                                         >  2:547479/1‑250 (MQ=255)
                                                                                                                                                            cATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATaaaa                                                                   <  1:212982/251‑1 (MQ=255)
                                                                                                                                                                     aaaTAAGCGGCTCTCAGTTTGTTTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAg                                                          <  2:558450/251‑1 (MQ=255)
                                                                                                                                                                                         gagtaacagtggTATTAATGCACGCTCCCGCCCCTGGTTTTACAGATACCAGAGCCCTAC‑GTGCATGCCTATGATGTTGCGAGCATTACCCACTCGCGGTAAAACCGAAAAATTCAAACGCTAATTGGCTTACCAATCCGCTATGTAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACGGCTGACTCGTTGTTATTCTGATGTTTTATATAAAACACCGGCAGGAAGACTCGaaacacaatac                                      <  2:330298/241‑7 (MQ=255)
                                                                                                                                                                                          aTTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAACCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATTTAAAACAACGGCAGGAAGATTCGCAACAAATTACt                                     >  2:584817/1‑251 (MQ=255)
                                                                                                                                                                                                agagGTATTAATGCCCTCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTtgct                               <  2:362878/251‑1 (MQ=255)
                                                                                                                                                                                                 gagGTATTAAGGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTCCCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTtgctg                              <  1:416809/251‑1 (MQ=255)
                                                                                                                                                                                                  agGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTtgctgc                             <  2:523606/251‑1 (MQ=255)
                                                                                                                                                                                                      aTTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAAt                         <  2:441674/251‑1 (MQ=255)
                                                                                                                                                                                                        tAATGCACGCTACCACCCCTGTCTTTACAGCTGCCAGAGCACTGCGAT‑CATGCCTATGATGTGACGAGCGTTACCCACTCGCACTAAACCCGCATAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGTAAACTAGGAAAATCCTGGAAAACCTTGACTAAACCCCTACTGCTAACTCGTTGTGATTCTGCTGGTTTATATAAAACAACGGCAGTACAATTCCCAACAAATTACtttttcttagaatat                       >  2:467781/1‑240 (MQ=255)
                                                                                                                                                                                                         aaTGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGCTTCGCAACAAATTACTTTTGCTGCGAAtttt                      >  2:12986/1‑251 (MQ=255)
                                                                                                                                                                                                                             gCTTTACAGCTACCAGAGCACTGC‑GTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTTCGAATTTTCACTGCCGTTATAATTTTCt  >  1:312638/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
CAAGTGATTCCCGGTTTCCGGCGTACGTTAATTTATCGACATTAACAACACTATCCTGCGTATTATTTATAATGTGACGAACTACAGCTGAACCAATAAATCCTGCGCCACCAGTAACAAGTATTTTCACTTAATTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGC‑ATGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCT  >  NC_000913/2112932‑2113403

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