Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,807,638 T→G 5.7% E154A (GAA→GCA)  waaQ ← lipopolysaccharide core biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,807,6380TG 5.7% 2796.5 / 128.4 820E154A (GAA→GCA) waaQlipopolysaccharide core biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base T (368/402);  new base G (25/23);  total (393/427)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.54e-01

Only 100 of 841 total aligned reads displayed.

CCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGTAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGCTTTTTTTCCAAATACCATGCTGCCGATGACCATAAAGTTGCGATATTTTCATCCGTGCAGGTAAACAACGTACCAGCAGCGCCACCATCCACTGATCCGTAAGATTAATGACCAGGTCATAATTATTCGCACGCAGAGTTTTTATCAACGAAAGCACATT  >  NC_000913/3807365‑3807912
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                  
          gAAAACGTGTTTTACCTGCAAGGCCAGTAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTc                                                                                                                                                                                                                                                                                   >  2:8680/1‑265 (MQ=255)
                                                                               atgCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTcaacaa                                                                                                                                                                                                      <  2:326322/271‑1 (MQ=255)
                                                                                           cTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAACTCGGTAATACCTAATGGCTCACGGACCGATAAATTACGCTCACTAATCTGTGTACCttgt                                                                                                                                                                                        >  1:11582/1‑275 (MQ=255)
                                                                                                            gTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCt                                                                                                                                                                               <  2:170548/267‑1 (MQ=255)
                                                                                                                      ccTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGc                                                                                                                                                                                                             >  4:237978/1‑227 (MQ=255)
                                                                                                                      ccTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCt                                                                                                                                                                                                            >  2:50277/1‑228 (MQ=255)
                                                                                                                      ccTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTa                                                                                                                                                                                       >  2:124991/1‑249 (MQ=255)
                                                                                                                      ccTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTa                                                                                                                                                                                       <  4:124991/249‑1 (MQ=255)
                                                                                                                      ccTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGCtttt                                                                                                                                                              >  4:298695/1‑274 (MQ=255)
                                                                                                                        tGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGaa                                                                                                                                                                                                                                                     <  4:73725/185‑1 (MQ=255)
                                                                                                                              ccTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTc                                                                                                                                                                                                                                                                                   >  5:498213/1‑149 (MQ=255)
                                                                                                                               cTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTGCTTCCAGCAATCTTCg                                                                                                                                                                                                                                                                                  >  7:1245095/1‑149 (MQ=255)
                                                                                                                                tCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCgg                                                                                                                                                                                                                                                                                 <  5:227442/149‑1 (MQ=255)
                                                                                                                                tCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCgg                                                                                                                                                                                                                                                                                 >  7:867001/1‑149 (MQ=255)
                                                                                                                                tCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCgg                                                                                                                                                                                                                                                                                 >  6:1287554/1‑149 (MQ=255)
                                                                                                                                tCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAAAACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTGCgg                                                                                                                                                                                                                                                                                 >  8:400599/1‑149 (MQ=255)
                                                                                                                                  gCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCg                                                                                                                                                                                                                                                                               >  7:1368761/1‑149 (MQ=255)
                                                                                                                                         cAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTGCGGCGTAACTCa                                                                                                                                                                                                                                                                        <  6:686614/149‑1 (MQ=255)
                                                                                                                                           gcgcATCGCTAACCTTAGAAAATTTATCGTTATCCCAACACCTAAATATCTGACGCGCTGTCGGTTGGGTGACAACATAATGATGTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCGTCTTCCAGCAATCTCCGGCGTAACTCAtt                                                                                                                                                                                                                                                                      <  8:202067/149‑1 (MQ=255)
                                                                                                                                           gcgcATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCAtt                                                                                                                                                                                                                                                                      <  8:759307/149‑1 (MQ=255)
                                                                                                                                           gcgcATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCAtt                                                                                                                                                                                                                                                                      <  6:198986/149‑1 (MQ=255)
                                                                                                                                           gcgcATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCAtt                                                                                                                                                                                                                                                                      <  5:158274/149‑1 (MQ=255)
                                                                                                                                             gcATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCAttgt                                                                                                                                                                                                                                                                    >  6:861864/1‑149 (MQ=255)
                                                                                                                                                   aTAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCg                                                                                                                                                                                                                                                                                  >  8:718693/1‑129 (MQ=255)
                                                                                                                                                   aTAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCAttgtt                                                                                                                                                                                                                                                                   <  7:1263882/144‑1 (MQ=255)
                                                                                                                                                    tAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCg                                                                                                                                                                                                                                                             <  5:538844/149‑1 (MQ=255)
                                                                                                                                                      aCCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGgt                                                                                                                                                                                                                                                           <  7:188226/149‑1 (MQ=255)
                                                                                                                                                        cTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATGTTCGGCGTAACTCATTGTTGTGTCGgtgt                                                                                                                                                                                                                                                         <  8:107553/149‑1 (MQ=255)
                                                                                                                                                        cTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGgtgt                                                                                                                                                                                                                                                         >  5:1219951/1‑149 (MQ=255)
                                                                                                                                                         ttAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTa                                                                                                                                                                                                                                                        >  5:273862/1‑149 (MQ=255)
                                                                                                                                                         ttAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTa                                                                                                                                                                                                                                                        <  5:1352179/149‑1 (MQ=255)
                                                                                                                                                           aGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAg                                                                                                                                                                                                                                                       >  4:377300/1‑148 (MQ=255)
                                                                                                                                                                  ttATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTGCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTca                                                                                                                                                                                                          >  2:21148/1‑186 (MQ=255)
                                                                                                                                                                   tATCGTTATCCCAAGACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTa                                                                                                                                                                                                                                              >  7:110483/1‑149 (MQ=255)
                                                                                                                                                                       gTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGTGGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATAc                                                                                                                                                                                                                                          >  7:624921/1‑149 (MQ=255)
                                                                                                                                                                           tCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTaa                                                                                                                                                                                                                                      <  6:335744/149‑1 (MQ=255)
                                                                                                                                                                             ccAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATg                                                                                                                                                                                                                                    >  5:417549/1‑149 (MQ=255)
                                                                                                                                                                              cAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTGCGGCGTAACTCAttgtt                                                                                                                                                                                                                                                                   >  8:1374619/1‑117 (MQ=255)
                                                                                                                                                                               aaCACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGc                                                                                                                                                                                                                                  <  6:57325/149‑1 (MQ=255)
                                                                                                                                                                                  acTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCa                                                                                                                                                                                                                               >  7:1425502/1‑149 (MQ=255)
                                                                                                                                                                                    ttAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCg                                                                                                                                                                                                                                                 >  5:1288664/1‑129 (MQ=255)
                                                                                                                                                                                    ttAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAg                                                                                                                                                                                                                             <  5:587408/149‑1 (MQ=255)
                                                                                                                                                                                      aaaTATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCgtg                                                                                                                                                                                         >  2:353116/1‑183 (MQ=255)
                                                                                                                                                                                          atCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGa                                                                                                                                                                                                                       <  6:481587/149‑1 (MQ=255)
                                                                                                                                                                                           tCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATtgttgt                                                                                                                                                                                                                                                                 >  6:812164/1‑106 (MQ=255)
                                                                                                                                                                                            cTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTaa                                                                                                                                                                                                                                             <  8:841177/125‑1 (MQ=255)
                                                                                                                                                                                            cTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATa                                                                                                                                                                                                                     >  7:762032/1‑149 (MQ=255)
                                                                                                                                                                                              gACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATaaa                                                                                                                                                                                                                   >  5:857630/1‑149 (MQ=255)
                                                                                                                                                                                               aCGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAAt                                                                                                                                                                                                                  <  6:740002/149‑1 (MQ=255)
                                                                                                                                                                                                  cgcTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTGCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTAc                                                                                                                                                                                                               <  8:431183/149‑1 (MQ=255)
                                                                                                                                                                                                    cTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGc                                                                                                                                                                                                             >  7:54455/1‑149 (MQ=255)
                                                                                                                                                                                                     tGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCt                                                                                                                                                                                                            <  6:435278/149‑1 (MQ=255)
                                                                                                                                                                                                     tGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCt                                                                                                                                                                                                            >  5:1192418/1‑149 (MQ=255)
                                                                                                                                                                                                             ggATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTACGGCGTAACACATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAata                                                                                                                                                                                                    <  6:1080403/149‑1 (MQ=255)
                                                                                                                                                                                                               aTGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATAtg                                                                                                                                                                                                  <  8:291033/149‑1 (MQ=255)
                                                                                                                                                                                                                 gACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATAtgtg                                                                                                                                                                                                >  6:1392917/1‑149 (MQ=255)
                                                                                                                                                                                                                 gACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGCTTTTTTTCCAAATACCATTCTGTCGATGCCCAGAAAGTTGCGATATTTCCATGCTTGCAGGAAAACAACGTAACAGCACCGACCCCATCTACTTAt                                                                  >  4:174672/1‑275 (MQ=255)
                                                                                                                                                                                                                    aacaTAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTAc                                                                                                                                                                                             >  8:587546/1‑149 (MQ=255)
                                                                                                                                                                                                                       aTAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCgt                                                                                                                                                                                          <  7:352097/149‑1 (MQ=255)
                                                                                                                                                                                                                        tAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTGCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCgtg                                                                                                                                                                                         >  8:677019/1‑149 (MQ=255)
                                                                                                                                                                                                                         aaTGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCaa                                                                                                                                                                                                                              <  6:1130839/111‑1 (MQ=255)
                                                                                                                                                                                                                           tGATCTTTTACGCCCAGGGCATCTAATTCGCGGCGGATCTTCTTCCAGCAATCTTAGGCGTAACTCATCGTTGTGTCGGGGTAGAGATCAGTAATACCTAATGGCTCAAGGACCGATAAGTTGAGCTCAACAATATGTGTACCGTGTAt                                                                                                                                                                                      >  8:905004/1‑149 (MQ=255)
                                                                                                                                                                                                                             aTCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAAtt                                                                                                                                                                                                                 >  2:252264/1‑120 (MQ=255)
                                                                                                                                                                                                                             aTCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTg                                                                                                                                                                                    >  6:1200752/1‑149 (MQ=255)
                                                                                                                                                                                                                                ttttACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTcaa                                                                                                                                                                                                         <  7:712356/125‑1 (MQ=255)
                                                                                                                                                                                                                                ttttACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTg                                                                                                                                                                                    >  7:663295/1‑146 (MQ=255)
                                                                                                                                                                                                                                ttttACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGgcg                                                                                                                                                                                 >  6:1028417/1‑149 (MQ=255)
                                                                                                                                                                                                                                ttttACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGgc                                                                                                                                                                                  >  8:492621/1‑148 (MQ=255)
                                                                                                                                                                                                                                ttttACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGgc                                                                                                                                                                                  <  7:492621/148‑1 (MQ=255)
                                                                                                                                                                                                                                 tttACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGNGTATTGGCGCTAAGGGTGTGAAGCTTTTTTTCCAACTACCATGCTGCCGATTACCATAAAGTTGCGATATTTTTNTCCGTGCAGTGAAACCACGTACCACAACCANNCACATCCACTGAACCGTAAGaatcaagaa                                                  >  3:48120/1‑267 (MQ=255)
                                                                                                                                                                                                                                      gCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAata                                                                                                                                                                                                    >  5:766805/1‑124 (MQ=255)
                                                                                                                                                                                                                                      gCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAgt                                                                                                                                                                           <  5:1420675/149‑1 (MQ=255)
                                                                                                                                                                                                                                        ccAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAgtgt                                                                                                                                                                         <  8:1425502/149‑1 (MQ=255)
                                                                                                                                                                                                                                        ccAGGGCATCCAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAgtgt                                                                                                                                                                         >  6:636007/1‑149 (MQ=255)
                                                                                                                                                                                                                                          aGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAgtgtgt                                                                                                                                                                       <  7:728110/149‑1 (MQ=255)
                                                                                                                                                                                                                                           gggCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTg                                                                                                                                                                                    <  7:157030/135‑1 (MQ=255)
                                                                                                                                                                                                                                           gggCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTg                                                                                                                                                                                    >  8:157030/1‑135 (MQ=255)
                                                                                                                                                                                                                                           gggCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGtgtgtg                                                                                                                                                                      >  6:251194/1‑149 (MQ=255)
                                                                                                                                                                                                                                               aTCTAATTCTCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTTTTGCGTCCATTTAGAAATCTGTAATACCTAATGGATCAAGTTCCGATAAATTACGCTCAACAATATCTGTATCGTGTTTTGGCGCTAAGTGTGTGCAGc                                                                                                                                                                  >  8:353255/1‑149 (MQ=255)
                                                                                                                                                                                                                                               aTCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGc                                                                                                                                                                  >  5:34920/1‑149 (MQ=255)
                                                                                                                                                                                                                                               aTCTAATTCCCGGCGCATCTTCTTCCAGCAATCTGCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGc                                                                                                                                                                  >  8:1250670/1‑149 (MQ=255)
                                                                                                                                                                                                                                                tCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCAt                                                                                                                                                                                                                                                                       <  5:1079997/47‑1 (MQ=255)
                                                                                                                                                                                                                                                tCTAATGCACGGCGCAGCTTCTTCAAGTACTCTTCGGCGGAACTCATTGTTGTGTCGGAGTAGAAATCGGTAATACCTAATTGCTCAAGGACCGATAAATTACGCTCATCAATATGTGTACCGTGTATTGGCGCTAAGTGTTTGAAGCTTTTTTTCCAAATACCATGCTGCCGATGACCATAAAGTTGCGATATTTCCATCCGTGCAGGTAAACAACGTACCAGCAGCGCCACCATCCACTGATCCGTAAGATTAATGACCAGGTCATAattatt                                   <  1:66142/275‑1 (MQ=255)
                                                                                                                                                                                                                                                  tAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCaa                                                                                                                                                                                                                              <  7:222779/86‑1 (MQ=255)
                                                                                                                                                                                                                                                  tAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGCttt                                                                                                                                                               >  6:1439897/1‑149 (MQ=255)
                                                                                                                                                                                                                                                   aaTTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATAtgt                                                                                                                                                                                                 >  8:1134782/1‑114 (MQ=255)
                                                                                                                                                                                                                                                    aTTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCgt                                                                                                                                                                                          <  7:912708/120‑1 (MQ=255)
                                                                                                                                                                                                                                                    aTTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGCttttt                                                                                                                                                             <  8:224556/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGAc                                                                                                                                                                                                                          >  7:384043/1‑87 (MQ=255)
                                                                                                                                                                                                                                                     ttCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTaa                                                                                                                                                                             <  8:639208/132‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttCCCGGCGCATCTTCTTCCAGCAATCTGCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGCtttttt                                                                                                                                                            <  5:943268/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       cccGGCGCATCTTCTTCCAGCAATCTGCGGCGTCACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGCTTTTTTTc                                                                                                                                                          <  8:1281822/149‑1 (MQ=255)
                                                                                                                                                                                                                                                         cGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGCTTTTTTTCCa                                                                                                                                                        <  7:299227/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             gcATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGCTTTTTTTCCAAATa                                                                                                                                                    <  5:1028417/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              cATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGCTTTTTTTCCAAATAc                                                                                                                                                   <  5:694320/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               aTCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGCTTTTTTTCCAAATAcc                                                                                                                                                  >  7:1172306/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               aTCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGCCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGCTTTTTTTCCAAATAcc                                                                                                                                                  >  6:908806/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         aGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATAtg                                                                                                                                                                                                  >  6:979174/1‑91 (MQ=255)
                                                                                                                                                                                                                                                                                ttCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGCTTTTTTTCCAAATACCATGCTGCCGATGACCAt                                                                                                                                 >  8:118458/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                ttCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGCTTTTTTTCCAAATACCATGCTGCCGATGACCAt                                                                                                                                 >  5:1139557/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                  
CCCAGTTCAGGAAAACGTGTTTTACCTGCAAGGCCAGTAATGGGTTTTGTTTCGCAACCTCGTGCAATCTCATCTACACAAGCGAGATCATCTGCCGAGGGCCCACAGGTTAGCACAACCTGATAGCCTCGCTGTTGCAGCGCATCGATAACCTTAGAAAATTTATCGTTATCCCAACACTTAAATATCTGACGCGCTGTCGGTTGGATGACAACATAATGATCTTTTACGCCCAGGGCATCTAATTCCCGGCGCATCTTCTTCCAGCAATCTTCGGCGTAACTCATTGTTGTGTCGGTGTAGAAATCGGTAATACCTAATGGCTCAAGGACCGATAAATTACGCTCAACAATATGTGTACCGTGTATTGGCGCTAAGTGTGTGAAGCTTTTTTTCCAAATACCATGCTGCCGATGACCATAAAGTTGCGATATTTTCATCCGTGCAGGTAAACAACGTACCAGCAGCGCCACCATCCACTGATCCGTAAGATTAATGACCAGGTCATAATTATTCGCACGCAGAGTTTTTATCAACGAAAGCACATT  >  NC_000913/3807365‑3807912

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