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

CAAGTGATTCCCGGTTTCCGGCGTACGTTAATTTATCGACATTAACAACACTATCCTGCGTATTATTTATAATGTGACGAACTACAGCTGAACCAATAAATCCTGCGCCACCAGTAACAAGTATTTTCACTTAATTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCATGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCT  >  NC_000913/2112932‑2113403
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
cttggacattacatactcctgtcaggtgtgatcaccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTgcgtgcat                                                                                                                                                                                                                               >  2:275014‑M2/136‑251 (MQ=255)
  tggacattacatactcctgtcaggtgtgatccccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggttttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTgcgtgcat                                                                                                                                                                                                                               <  2:13993‑M2/116‑1 (MQ=255)
    gacattacatactcctgtcaggtgtgatcaccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCt                                                                                                                                                                                                                           <  1:493341‑M2/120‑1 (MQ=255)
          acatactcctgtcaggtgtgatcaccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTatgatg                                                                                                                                                                                                                     >  1:13016‑M2/126‑251 (MQ=255)
                    gtcaggtgtgatcaccgcgctgaatgctcgcggcctttaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgtttacggctttatTTATTCCATATTACTTCGGAGCATGCTTTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCATGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTTCATGCCTATGATGTTACGAGCGt                                                                                                                                                                                                           <  2:45499‑M2/136‑1 (MQ=255)
                     tcaggtgtgatcaccgcgctgcatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGTGGTATTAATGCACGCTACCGACCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGtt                                                                                                                                                                                                          >  2:255133‑M2/115‑251 (MQ=255)
                        ggcgtgatcaccgcgctgaatgctcgcggcctggaatgttcgcaaaaaagtatacaaagtacttagtaaagtataaatagcttagccatgcttatattcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTAGTAAGGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTAcc                                                                                                                                                                                                       <  2:253232‑M2/140‑1 (MQ=255)
                                caccgcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTcgc                                                                                                                                                                                                <  1:236531‑M2/147‑1 (MQ=255)
                                    gcgctgaatgctcgcggcctggaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTTCTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTaa                                                                                                                                                                                           <  1:572010‑M2/152‑1 (MQ=255)
                                                         gaatgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTa                                                                                                                                                                      <  1:259800‑M2/173‑1 (MQ=255)
                                                           atgttcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAAt                                                                                                                                                                    <  1:275256‑M2/175‑1 (MQ=255)
                                                               tcgcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGt                                                                                                                                                                 <  1:176023‑M2/178‑1 (MQ=255)
                                                                 gcaaataagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCtt                                                                                                                                                              <  1:275014‑M2/181‑1 (MQ=255)
                                                                       aagtatacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAAt                                                                                                                                                        >  1:305440‑M2/65‑251 (MQ=255)
                                                                           atacaaagtacttattaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCaatcagc                                                                                                                                                    >  2:423167‑M2/61‑248 (MQ=255)
                                                                                       tattaaactataaatagcttatccctgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGCGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCACTCCGCTCTGGAAACa                                                                                                                                          >  2:247209‑M2/49‑249 (MQ=255)
                                                                                        attaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAgg                                                                                                                                       >  2:210268‑M2/48‑251 (MQ=255)
                                                                                        attaaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAgg                                                                                                                                       <  2:148192‑M2/204‑1 (MQ=255)
                                                                                           aaagtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCACTCCGCTCTGGAACCAAGGaa                                                                                                                                     >  2:416809‑M2/45‑250 (MQ=255)
                                                                                             cgtataaatagcttatccatgcttatatgcttacggctttatTTATTCCATATTAATTCAGAGCATGCTGAGAACTAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCCCCCCTGGCTTTACAGCTACCAGTGCACTGCGTGCATGCCTATGATGTGACGAGCGGTACCCACTCGCGCTAACCCCGAACAACTCACCCGTTCATTGGCTTACCAATCCGCTCTGGACacaaggcata                                                                                                                                   >  2:348805‑M2/43‑246 (MQ=255)
                                                                                                   aatagcttatccatgcttattggcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCATTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGa                                                                                                                            <  2:62325‑M2/215‑1 (MQ=255)
                                                                                                         ttatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACtt                                                                                                                      >  1:533935‑M2/31‑251 (MQ=255)
                                                                                                         ttatccatgattatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACtt                                                                                                                      >  1:66219‑M2/31‑251 (MQ=255)
                                                                                                          tatccatgcttatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGTTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACtt                                                                                                                      <  2:293006‑M2/221‑1 (MQ=255)
                                                                                                                    tatatgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAAcc                                                                                                           >  2:374560‑M2/20‑251 (MQ=255)
                                                                                                                       atgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCt                                                                                                         <  1:210268‑M2/234‑1 (MQ=255)
                                                                                                                       atgcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCt                                                                                                         >  1:580467‑M2/17‑250 (MQ=255)
                                                                                                                         gcttacggctttatTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATACAGGTATTAATGCACTCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTCCCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGATTAAAACCCTAc                                                                                                       <  1:47624‑M2/236‑1 (MQ=255)
                                                                                                                                     aTTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCgtt                                                                                           <  1:276190/250‑1 (MQ=255)
                                                                                                                                                 ttACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATgg                                                                              >  2:179468/1‑251 (MQ=255)
                                                                                                                                                   aCTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGtt                                                                            >  1:35843/1‑251 (MQ=255)
                                                                                                                                                       cAGAGCATGCTGTGAAATACGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACAGCCACTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTtat                                                                         >  2:547479/1‑250 (MQ=255)
                                                                                                                                                            cATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATaaaa                                                                   <  1:212982/251‑1 (MQ=255)
                                                                                                                                                                     aaaTAAGCGGCTCTCAGTTTGTTTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAg                                                          <  2:558450/251‑1 (MQ=255)
                                                                                                                                                                                         gagtaacagtggTATTAATGCACGCTCCCGCCCCTGGTTTTACAGATACCAGAGCCCTACGTGCATGCCTATGATGTTGCGAGCATTACCCACTCGCGGTAAAACCGAAAAATTCAAACGCTAATTGGCTTACCAATCCGCTATGTAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACGGCTGACTCGTTGTTATTCTGATGTTTTATATAAAACACCGGCAGGAAGACTCGaaacacaatac                                      <  2:330298/241‑7 (MQ=255)
                                                                                                                                                                                          aTTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAACCCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATTTAAAACAACGGCAGGAAGATTCGCAACAAATTACt                                     >  2:584817/1‑251 (MQ=255)
                                                                                                                                                                                                agagGTATTAATGCCCTCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTtgct                               <  2:362878/251‑1 (MQ=255)
                                                                                                                                                                                                 gagGTATTAAGGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTCCCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTtgctg                              <  1:416809/251‑1 (MQ=255)
                                                                                                                                                                                                  agGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTtgctgc                             <  2:523606/251‑1 (MQ=255)
                                                                                                                                                                                                      aTTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAAt                         <  2:441674/251‑1 (MQ=255)
                                                                                                                                                                                                        tAATGCACGCTACCACCCCTGTCTTTACAGCTGCCAGAGCACTGCGATCATGCCTATGATGTGACGAGCGTTACCCACTCGCACTAAACCCGCATAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGTAAACTAGGAAAATCCTGGAAAACCTTGACTAAACCCCTACTGCTAACTCGTTGTGATTCTGCTGGTTTATATAAAACAACGGCAGTACAATTCCCAACAAATTACtttttcttagaatat                       >  2:467781/1‑240 (MQ=255)
                                                                                                                                                                                                         aaTGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGCTTCGCAACAAATTACTTTTGCTGCGAAtttt                      >  2:12986/1‑251 (MQ=255)
                                                                                                                                                                                                                             gCTTTACAGCTACCAGAGCACTGCGTGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTTCGAATTTTCACTGCCGTTATAATTTTCt  >  1:312638/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  
CAAGTGATTCCCGGTTTCCGGCGTACGTTAATTTATCGACATTAACAACACTATCCTGCGTATTATTTATAATGTGACGAACTACAGCTGAACCAATAAATCCTGCGCCACCAGTAACAAGTATTTTCACTTAATTTATTCCATATTACTTCAGAGCATGCTGTGAAATAAGCGGCTCTCAGTTTGATTAATAGAGGTATTAATGCACGCTACCGCCCCTGGCTTTACAGCTACCAGAGCACTGCATGCATGCCTATGATGTGACGAGCGTTACCCACTCGCGCTAAACCCGAAAAATTCAAACGCTAATTGTCTTACCAATCCGCTCTGGAAACAAGGAAAATCCTGGAAAACTTTGAATAAAACCCTACTGCTAACTCGTTGTTATTCTGATGGTTTATATAAAACAACGGCAGGAAGATTCGCAACAAATTACTTTTGCTGCGAATTTTCACTGCCGTTATAATTTTCT  >  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: