Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,789,083 A→C 5.5% M217L (ATG→CTG)  aroH → 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,789,0830AC 5.5% 454.1 / 15.8 128M217L (ATG→CTG) aroH3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible
Reads supporting (aligned to +/- strand):  ref base A (64/57);  new base C (3/4);  total (67/61)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.07e-01

Only 100 of 136 total aligned reads displayed.

CGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAAGCAGCACCGTCGCCAGTTAGAAGTTTGTGAGGATATTTGTCAGCAAATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAAC  >  NC_000913/1788810‑1789357
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                  
cccAAATTACTTTTTCGGGGAAATTAGCTGGGCGGCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGAGATTTGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGAt                                                                                                                                                                                                                                                                                   <  4:221023/273‑1 (MQ=255)
     aTTACTTTGACAGTAAATGGGCTGGGCGGCCAAACCGCGACCGAGTTCCGCGATATGGTGCCCGGTCAGTTTGTTGCTGATTTAATCGGTTGTGGCGCGATTGGCGCACGTACTACCGACAGTCAGAGCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAACAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTAACCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCa                                                                                                                                                                                                                                                                              <  4:75408/274‑1 (MQ=255)
                              ggCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                                                                                                                                                                                                                                                     >  2:389324/1‑275 (MQ=255)
                                    ccAACCGCGACCGAGTTCCTGGATATGGTGACCGGTCAGTTTAATGCTGCTTTAATCAGTTGGGGCGCGATTGGCGAACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCca                                                                                                                                                                                                                                               <  4:449164/275‑1 (MQ=255)
                                      aaccttgACCGAGTTCTTCGATATGGTTGCCGGTCAGTTTATTGCTCATTGTATCATTCGGGGCGGGATTGGCGAGCGTACAACCGAAAGGCTGATCCACCGGGATATGGCTTCGGCACTGCCCTCTCCGGAAGGTTTTAAGAATGGATCCGGTGGCAAGACGTGGATTGTTGTGGATGCTAGCCGCGCAGCCCGCTGCAGCCATTTGTTCCTCTGGCCTGGCAAAAAGGGTCAGATGACCATCTATCAGACCAGCGGCAACCTGTATGGcta                                                                                                                                                                                                                                               <  4:150653/267‑3 (MQ=255)
                                                gnattcgtgattatGGTGACCGGTCAGGTTATTGATTATTTAATCAGTTGGGGCTGGATTGTCGCACGTACTACCGAAAGTCAGAGCCACCGCGAAATGGCTTCGGCANTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCg                                                                                                                                                                                                                                                               <  3:108744/236‑1 (MQ=255)
                                                                     ggTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                                                                                                                                                                                                                                                    <  1:95280/237‑1 (MQ=255)
                                                                                     aTTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATc                                                                                                                                                                                              >  1:67159/1‑275 (MQ=255)
                                                                                     aTTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATc                                                                                                                                                                                              >  2:229472/1‑275 (MQ=255)
                                                                                     aTTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATc                                                                                                                                                                                              >  2:248058/1‑275 (MQ=255)
                                                                                                 ggggCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGc                                                                                                                                                                                  >  2:513367/1‑275 (MQ=255)
                                                                                                   ggCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGa                                                                                                                                                                                >  2:490498/1‑275 (MQ=255)
                                                                                                             gggcACGTCTTACCGAAAGTCAGATCAACGGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCa                                                                                                                                                                       <  4:490497/272‑1 (MQ=255)
                                                                                                                        acCGAAGGTCAGATCAACCGCGAAAGGGGTTCGGCAACGTCCGGTCCGGTAGGTTTTAAAAATGGTACAGATGGAAATACGCGGATGGCTGTGGAGGCTATACGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGGCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGAtt                                                                                                                                                            <  4:121381/274‑1 (MQ=255)
                                                                                                                            aaaGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATTATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCt                                                                                                                                                       >  2:18032/1‑275 (MQ=255)
                                                                                                                                   gATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGAcc                                                                                                                                                                                                                                                                               >  5:860912/1‑148 (MQ=255)
                                                                                                                                   gATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCTTGCGATACGCTGACCGAGTTTGAGTTACCttgaactc                                                                                                                                                >  4:474108/1‑268 (MQ=255)
                                                                                                                                     tCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCAt                                                                                                                                                                                                                                                                             >  6:185156/1‑148 (MQ=255)
                                                                                                                                      ccACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATc                                                                                                                                                                                                                                                                            <  6:860912/148‑1 (MQ=255)
                                                                                                                                       cACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCATCTGACCAttt                                                                                                                                                                                                                                                                           <  8:92131/148‑3 (MQ=255)
                                                                                                                                          cgcgAAATGGTTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtggtggt                                                                                                                                         <  4:389323/275‑1 (MQ=255)
                                                                                                                                           gcgAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCa                                                                                                                                                                                                                                                                       >  5:290939/1‑148 (MQ=255)
                                                                                                                                           gcgAAATGGCTTCGGCACTCTCATGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCGATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtggtggt                                                                                                                                         <  1:45911/274‑1 (MQ=255)
                                                                                                                                             gAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGa                                                                                                                                                                                                                                                                     <  5:14847/148‑1 (MQ=255)
                                                                                                                                             gAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGa                                                                                                                                                                                                                                                                     <  7:259477/148‑1 (MQ=255)
                                                                                                                                              aaaTGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACCTCGGGTGGTGGa                                                                                                                                      >  2:352190/1‑274 (MQ=255)
                                                                                                                                                  ggCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                                                                                                                                                                                                                                                    >  8:732482/1‑128 (MQ=255)
                                                                                                                                                  ggCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                                                                                                                                                                                                                                                    <  7:732482/128‑1 (MQ=255)
                                                                                                                                                           aCTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCg                                                                                                                                                                                                                                                       <  6:290939/148‑1 (MQ=255)
                                                                                                                                                              ctcCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTAt                                                                                                                                                                                                                                                    >  5:50025/1‑148 (MQ=255)
                                                                                                                                                                      cGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAag                                                                                                             >  2:213640/1‑275 (MQ=255)
                                                                                                                                                                       ggTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACatt                                                                                                                                                                                                                                           <  6:50025/148‑1 (MQ=255)
                                                                                                                                                                         tAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCCGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCACCGGCCACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGAAGATATCGCCGCCGCCTGCGTTACCCTGCACGAGTTTGATTTACCTGACCCACTGGTGGGGGTTTTCAGCCCCGGTAACTGccctcagcc                                                                                                           >  4:528686/1‑267 (MQ=255)
                                                                                                                                                                         tAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGTGTTTGATTTACCTGAACATCGGGTGGTTGATTTCAGCCACGGTAACCGCCAGAaacag                                                                                                          >  4:236491/1‑275 (MQ=255)
                                                                                                                                                                         tAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAagcag                                                                                                          >  1:110973/1‑275 (MQ=255)
                                                                                                                                                                             ttttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTAt                                                                                                                                                                                                                                                                         >  8:864148/1‑112 (MQ=255)
                                                                                                                                                                              tttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                                                                                                                                                                                                                                                    >  7:559221/1‑100 (MQ=255)
                                                                                                                                                                              tttAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                                                                                                                                                                                                                                                    <  8:559221/100‑1 (MQ=255)
                                                                                                                                                                                 aaaaaTGTGACCGATGGCAATACGCGGATTGCTGTGGAGGCTTTCAGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATTATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAAGCAGCACCGTCg                                                                                                  <  4:18032/275‑1 (MQ=255)
                                                                                                                                                                                 aaaaaTGGTACCGATGGCAATACGCGGAGTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAAGCAGCACCGTCg                                                                                                  <  2:474109/275‑1 (MQ=255)
                                                                                                                                                                                    aaTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGt                                                                                                                         <  2:170086/249‑1 (MQ=255)
                                                                                                                                                                                     aTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                                                                                                                                                                                                                                                     >  5:480596/1‑124 (MQ=255)
                                                                                                                                                                                     aTGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                                                                                                                                                                                                                                                     <  6:480596/124‑1 (MQ=255)
                                                                                                                                                                                      tGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATTACCATCTATCAGACCAGCGGCaa                                                                                                                                                                                                                                                           >  5:193120/1‑117 (MQ=255)
                                                                                                                                                                                      tGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCaa                                                                                                                                                                                                                                                           <  6:193120/117‑1 (MQ=255)
                                                                                                                                                                                       ggTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGCTGATATCGCGGCAGCCTGCGATCCCCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGTTTTCAGCCACGGTAACTGCCAGAGGCATCCCCGTCGCCAGtt                                                                                            >  4:193480/1‑275 (MQ=255)
                                                                                                                                                                                         tACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACAttatta                                                                                                                                                                                                                                       <  8:14394/134‑1 (MQ=255)
                                                                                                                                                                                         tACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACAttatta                                                                                                                                                                                                                                       >  7:14394/1‑134 (MQ=255)
                                                                                                                                                                                              aTGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAAcc                                                                                                                                                                                                                    >  6:473837/1‑148 (MQ=255)
                                                                                                                                                                                                    aTACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTa                                                                                                                                                                                                              >  5:191203/1‑148 (MQ=255)
                                                                                                                                                                                                     tACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCg                                                                                                                                                                                                                                                               >  2:544477/1‑98 (MQ=255)
                                                                                                                                                                                                      aCGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAg                                                                                                                                                                                                                                                                      <  7:304460/90‑1 (MQ=255)
                                                                                                                                                                                                      aCGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAg                                                                                                                                                                                                                                                                      >  8:304460/1‑90 (MQ=255)
                                                                                                                                                                                                       cgcgGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCaaaaa                                                                                                                                                                                                                       >  5:92066/1‑136 (MQ=255)
                                                                                                                                                                                                       cgcgGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCaaaaa                                                                                                                                                                                                                       <  6:92066/136‑1 (MQ=255)
                                                                                                                                                                                                         cgGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCaa                                                                                                                                                                                                                                                           <  4:544476/98‑1 (MQ=255)
                                                                                                                                                                                                         cgGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATg                                                                                                                                                                                                         <  5:473837/148‑1 (MQ=255)
                                                                                                                                                                                                           gATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTAt                                                                                                                                                                                                                                                                         >  7:588544/1‑82 (MQ=255)
                                                                                                                                                                                                           gATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTAt                                                                                                                                                                                                                                                                         <  8:588544/82‑1 (MQ=255)
                                                                                                                                                                                                            aTTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCg                                                                                                                                                                                                      <  6:191203/148‑1 (MQ=255)
                                                                                                                                                                                                             ttGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCga                                                                                                                                                                                                     <  5:848269/148‑1 (MQ=255)
                                                                                                                                                                                                                    ggATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                                                                                                                                                                                                                                                    >  5:1125037/1‑62 (MQ=255)
                                                                                                                                                                                                                    ggATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGa                                                                                                                                                                                                                                                                                    <  6:1125037/62‑1 (MQ=255)
                                                                                                                                                                                                                     gATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCAAGGTAACTGCCAGAAGCAGCACCGTCGCCCGTTAGCAGTTTGTGAGGATATTTGTCAGCaaa                                                               >  2:118403/1‑274 (MQ=255)
                                                                                                                                                                                                                      aTGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCATCTGACCAtt                                                                                                                                                                                                                                                                            >  7:978552/1‑67 (MQ=255)
                                                                                                                                                                                                                       tGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGatat                                                                                                                                                                                               <  7:71323/144‑1 (MQ=255)
                                                                                                                                                                                                                          tATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACattatt                                                                                                                                                                                                                                        <  7:483713/100‑1 (MQ=255)
                                                                                                                                                                                                                             ccGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATg                                                                                                                                                                                                                                     >  6:492177/1‑100 (MQ=255)
                                                                                                                                                                                                                                cgcAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                                                                                                                                                                                                                                                     <  6:1037713/81‑1 (MQ=255)
                                                                                                                                                                                                                                cgcAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTa                                                                                                                                                                                                                                                     >  5:1037713/1‑81 (MQ=255)
                                                                                                                                                                                                                                        gcgcCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAAGCAGCACCGTCGCCAGTTAGAAGTTTGTGAGGATATTTGTCAGCAAATCCGCAATGGCTCTACGGc                                            <  2:412448/274‑1 (MQ=255)
                                                                                                                                                                                                                                        gcgcCAGCCATAAGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCaaaaa                                                                                                                                                                                                                       <  7:963651/103‑1 (MQ=255)
                                                                                                                                                                                                                                        gcgcCAGCCATAAGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCaaaaa                                                                                                                                                                                                                       >  8:963651/1‑103 (MQ=255)
                                                                                                                                                                                                                                              gccaGATGTTCCTCTAGCCAAACAAAAATAGACAAAGGACCAGTTAAAGGAACAGGGGCAACCCGTATGGCAACATTATTAGGCGTGGCGGCAAAAAACCGAATTATCATGCGGATGATATAGCCGCAGGCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAAGCAGCACCGTCGCCAGTTAGAAGTTTGTGAGGATATTGGTCAGCAAATCCGCAATGGCTCTACGGCGATTGc                                      <  4:155581/274‑1 (MQ=255)
                                                                                                                                                                                                                                                      ttCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGa                                                                                                                                                                                                                  >  8:200545/1‑94 (MQ=255)
                                                                                                                                                                                                                                                      ttCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGa                                                                                                                                                              <  6:856742/146‑1 (MQ=255)
                                                                                                                                                                                                                                                           ctcGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCa                                                                                                                                                                                                           >  1:13569/1‑96 (MQ=255)
                                                                                                                                                                                                                                                           ctcGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATcg                                                                                                                                                                                             >  8:361876/1‑110 (MQ=255)
                                                                                                                                                                                                                                                               ccAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATg                                                                                                                                                                                                         <  8:78691/94‑1 (MQ=255)
                                                                                                                                                                                                                                                               ccAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATg                                                                                                                                                                                                         >  7:78691/1‑94 (MQ=255)
                                                                                                                                                                                                                                                               acAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGcc                                                                                                                                                                                                       <  3:13569/95‑1 (MQ=255)
                                                                                                                                                                                                                                                                cAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACa                                                                                                                                                  >  7:1099813/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                 aGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCa                                                                                                                                                                                                           <  6:706497/90‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCa                                                                                                                                                                                                           >  5:706497/1‑90 (MQ=255)
                                                                                                                                                                                                                                                                 aGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGa                                                                                                                                                     >  6:296232/1‑144 (MQ=255)
                                                                                                                                                                                                                                                                 aGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGa                                                                                                                                                     <  5:296232/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                  gACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATtattat                                                                                                                                                                                                                                      <  5:485842/62‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aCAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTAcc                                                                                                                                                        <  5:578548/139‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aCAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCt                                                                                                                                               >  8:636058/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                        gaTGGTCATCTGACCATCTATCAGACCAGCGGCAACCCGTATTGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTg                                                                                                                                                               >  5:502769/2‑127 (MQ=255)
                                                                                                                                                                                                                                                                        gaTGGTCATCTGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAATCCGAAGTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTg                                                                                                                                                               <  6:502769/126‑1 (MQ=255)
                                                                                                                                                                                                                                                                        gaTGGTCATCTGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCgatga                                                                                                                                                                                                  >  6:1024710/2‑92 (MQ=255)
                                                                                                                                                                                                                                                                        aaTGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCgatg                                                                                                                                                                                                   >  6:391235/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                        aaTGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCgatg                                                                                                                                                                                                   <  5:391235/91‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aaTGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtggtgg                                                                                                                                          >  5:986174/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                         aTGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCa                                                                                                                                                                       >  7:879740/1‑118 (MQ=255)
                                                                                                                                                                                                                                                                         aTGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtggtggt                                                                                                                                         >  8:288159/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                          tGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCtggtggtg                                                                                                                                        >  6:633271/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                gATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTa                                                                                                                                                          >  6:699402/1‑124 (MQ=255)
                                                                                                                                                                                                                                                                                 aTGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCa                                                                                                                                 <  7:288159/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                  
CGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGTTTTAAAAATGGTACCGATGGCAATACGCGGATTGCTGTGGATGCTATCCGCGCAGCCCGCGCCAGCCATATGTTCCTCTCGCCAGACAAAAATGGTCAGATGACCATCTATCAGACCAGCGGCAACCCGTATGGCCACATTATTATGCGTGGCGGCAAAAAACCGAATTATCATGCCGATGATATCGCCGCAGCCTGCGATACGCTGCACGAGTTTGATTTACCTGAACATCTGGTGGTGGATTTCAGCCACGGTAACTGCCAGAAGCAGCACCGTCGCCAGTTAGAAGTTTGTGAGGATATTTGTCAGCAAATCCGCAATGGCTCTACGGCGATTGCTGGAATTATGGCGGAAAGTTTCCTGCGCGAAGGAAC  >  NC_000913/1788810‑1789357

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