Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,131,524 A→G 10.3% S295S (AGT→AGC) ‡ yghQ ← putative inner membrane polysaccharide flippase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,131,5240AG10.3% 220.5 / 16.0 78S295S (AGT→AGC) ‡yghQputative inner membrane polysaccharide flippase
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base G (4/4);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-04

TACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGCAGACTCAAACAGATTCAATTTGAAGGCGTTATGGATATTTCGGCGGCGTAATTCGCGCGCGGCAAACCAC  >  NC_000913/3131256‑3131785
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                     
tACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAAGGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAAGCAGGTCATACGCTTCGAGATACTTAACGCCAAACCCCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATGCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                                                                                                                                                                                                 <  2:343262/275‑1 (MQ=255)
        aaaaaaTCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATCCCGATTCCGCCCGCCAGTAAACCGGACTTCAAACCTAACAACCACGATCTGTGGGTGCGCGGATCCAAACGCAACCCCCATGGGTAAAAACTTTTACCCAGCaa                                                                                                                                                                                                                                                         >  4:198988/1‑275 (MQ=38)
                                aatgaaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCCTCGAGATACTTAACGCCAAACACCAGCGAAATGCGCGGCTTGCCCACGATCAACACCGCTAGCGCCACCAGACTACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGATCTGGTGGTTCGCGGATCTAAACGCATCACCCATGGGTAAAAACTTTTACCCAGCAAACCCGCCGGCGTTCCgacggcct                                                                                                                                                                                                                                  >  2:462865/1‑272 (MQ=255)
                                                                       tATTCGCACCTTCCCGAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACAGCGGAAATGAGCGGCTTGCCAACAATCAACACCGCGAGCGCAACCAGAATACCGATTCCGACCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGc                                                                                                                                                                                          <  2:51481/275‑1 (MQ=255)
                                                                                  tCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAACCTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAaca                                                                                                                                                                               >  4:36830/1‑275 (MQ=255)
                                                                                       aaaTAAGGTCATACGCTTCGGAATACTTAACGCCAAAAACCAGCGAAATGGGCGGCTTGCCGACCATCAAAACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATc                                                                                                                                                                          <  2:215473/275‑1 (MQ=255)
                                                                                                               aCTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAAACCGGCGGCAGCGGGGCCTAACACGACCCCCAACAACACTGTGCTGCACGAg                                                                                                                                                  >  1:112123/1‑275 (MQ=255)
                                                                                                                   aaCGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGtgcccgtcccttat                                                                                                                                               >  2:117913/1‑263 (MQ=255)
                                                                                                                       ccAAACACAAGCGAAAGGAGCGGCTTGCCGACGATCACCACCGCGAGCGCAACCAAAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCGAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAGCACTGTGCTGCACGAGTTACgcgc                                                                                                                                          <  2:198989/275‑1 (MQ=255)
                                                                                                                                 gCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTAcc                                                                                                                                                                                                                                                               >  8:448235/1‑148 (MQ=255)
                                                                                                                                                         aTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcgg                                                                                                                                                                                                                                        >  6:550059/1‑148 (MQ=255)
                                                                                                                                                         aTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACAGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTAccc                                                                                                                                                                                                                                                              <  5:110599/126‑1 (MQ=255)
                                                                                                                                                         aTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACAGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTAccc                                                                                                                                                                                                                                                              <  5:522911/126‑1 (MQ=255)
                                                                                                                                                         aTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACAGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTAccc                                                                                                                                                                                                                                                              >  6:110599/1‑126 (MQ=255)
                                                                                                                                                         aTCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACAGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTAccc                                                                                                                                                                                                                                                              >  6:522911/1‑126 (MQ=255)
                                                                                                                                                                 cgcgAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                                                                                                                                                                                                 >  5:504468/1‑115 (MQ=255)
                                                                                                                                                                 cgcgAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                                                                                                                                                                                                 <  6:504468/115‑1 (MQ=255)
                                                                                                                                                                  gcgAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCcggcg                                                                                                                                                                                               >  2:74512/1‑180 (MQ=255)
                                                                                                                                                                   cgAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAACACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGTTCCCCACAAACACTGTGCTGCACGACTTACGCCCACAACCGACGGACTGGGCAACGTTGGGTGACCAGTCAAAACt                                                                                                >  4:74068/1‑274 (MQ=255)
                                                                                                                                                                      gcgcCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaagaa                                                                                                                                                                                                                           <  7:448235/148‑1 (MQ=255)
                                                                                                                                                                      gcgcCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAg                                                                                                                                                                                           <  4:74511/180‑1 (MQ=255)
                                                                                                                                                                          gaccaAAATACCGTACCCGGCCGCAAGTAAACCGGACCTCACACCTAACAACCACGGTCTGGCGGGGCGCGGATCTAAACGCATCACCTCTGGGTAAACACTTTTACCCAGCAACCCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACgcgc                                                                                        <  4:481293/274‑1 (MQ=255)
                                                                                                                                                                           accaaaaccacGATTCCCGCCGCCACTAGACCCGACTTCACACCTTACATCCCCGGTCTGGTGGTGCGCGGATCTAACCGCATCACCTCTGCGGAACAACTTTTACCAGGCAAACCCGCTGGCCTTCCGGCTGCGTCGAAGAATGGCATGGCGATTTTAAGTAACCCGGCGGCATCGGGTCCTAACACGATCCCGAACAACACTGTGCTGCACGAGTTACGCGCCGAGCAGATGGTGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGcc                                                                                       <  4:435326/265‑1 (MQ=255)
                                                                                                                                                                            ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACtttt                                                                                                                                                                                                                                                                  >  6:32900/1‑103 (MQ=255)
                                                                                                                                                                            ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACtttt                                                                                                                                                                                                                                                                  <  7:365927/103‑1 (MQ=255)
                                                                                                                                                                            ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACtttt                                                                                                                                                                                                                                                                  <  5:32900/103‑1 (MQ=255)
                                                                                                                                                                            ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACtttt                                                                                                                                                                                                                                                                  >  8:365927/1‑103 (MQ=255)
                                                                                                                                                                            ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                                                                                                                                                                                                 >  7:628024/1‑104 (MQ=255)
                                                                                                                                                                            ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                                                                                                                                                                                                 <  8:628024/104‑1 (MQ=255)
                                                                                                                                                                            ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTAcc                                                                                                                                                                                                                                                               >  6:124169/1‑106 (MQ=255)
                                                                                                                                                                            ccaGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTAcc                                                                                                                                                                                                                                                               <  5:124169/106‑1 (MQ=255)
                                                                                                                                                                                      gATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                                                                                                                                                                                                 <  8:979352/94‑1 (MQ=255)
                                                                                                                                                                                      gATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTa                                                                                                                                                                                                                                                                 >  7:979352/1‑94 (MQ=255)
                                                                                                                                                                                       aTTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTACAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAAAAAGATCCCAACCAACACTGTGCTGCACGAGTTAAGCGACGACAAGATGGAGTGGGCAATGTTGGTTGCCCAGCCAAAACACAACGCGACTTTAATGTACCg                                                                           >  1:115341/1‑275 (MQ=255)
                                                                                                                                                                                          ccgcccgcCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGc                                                                        >  1:104591/1‑275 (MQ=255)
                                                                                                                                                                                            gcccgccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCCATGTTTGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGCaa                                                                      >  4:343951/1‑274 (MQ=255)
                                                                                                                                                                                               cgccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCAga                                                                                                                                                                                                                                                           >  6:873866/1‑90 (MQ=255)
                                                                                                                                                                                               cgccAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCAga                                                                                                                                                                                                                                                           <  5:873866/91‑2 (MQ=255)
                                                                                                                                                                                                   agaaaaCCGTACTATACCCCTACCACCCCNTGTCTGGTGGGGCCCGGACCTAACCNCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAANTAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGCAGACTCaaa                                                               <  3:104591/272‑1 (MQ=255)
                                                                                                                                                                                                     tAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTACCGGCGGCGTCgaa                                                                                                                                                                                                                              <  6:513107/114‑1 (MQ=255)
                                                                                                                                                                                                     tAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGAGCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaa                                                                                                                                                                                                                              >  5:513107/1‑114 (MQ=255)
                                                                                                                                                                                                        aCCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaa                                                                                                                                                                                                                              >  5:515109/1‑111 (MQ=255)
                                                                                                                                                                                                        aCCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaa                                                                                                                                                                                                                              <  6:515109/111‑1 (MQ=255)
                                                                                                                                                                                                            gACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACACGCTGGCGTTCAGGCGGCTTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCCCAGACCAGAAGGACGTGGCAATGTTGGTTTACCAGACAAAACCACACACCCCTTTAATGTATCGGGCAGACTCAAACAGATTCac                                                      >  4:317989/1‑274 (MQ=255)
                                                                                                                                                                                                                tCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGt                                                                                                                                                                                                                        <  5:337991/109‑1 (MQ=255)
                                                                                                                                                                                                                tCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGt                                                                                                                                                                                                                        >  6:337991/1‑109 (MQ=255)
                                                                                                                                                                                                                 cacaCCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCACCACTGTTCTGCCCGAGCTACGCGCCGACCAGATGGAGTGGGCCATGTTTGTTTACCAGACCAAACTCCACCGGCCTTTAATGTTGCGGGCAGACCCCAACAGAACCCATTTGa                                                 >  4:146237/1‑275 (MQ=255)
                                                                                                                                                                                                                   cacCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcgg                                                                                                                                                                                                                                        <  7:1074093/90‑1 (MQ=255)
                                                                                                                                                                                                                   cacCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCcgg                                                                                                                                                                                                                                        >  8:1074093/1‑90 (MQ=255)
                                                                                                                                                                                                                    acctactcacCTCGGTCGGGGGGTGAGCCGGTCTCAACGGCTCACCTTCGGGTAAAAACTTTTCACCAGCCAAACCGCTGGCCTTCCGGCGGCGTCTAAGAATGTCCTGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTCACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGCAGACTCAAACAGATTCAATTTGAAg                                               <  4:145491/266‑1 (MQ=255)
                                                                                                                                                                                                                    acCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGCAGACTCAAACAGATTCAATTTGAAg                                               <  2:36831/274‑1 (MQ=255)
                                                                                                                                                                                                                       tAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCAga                                                                                                                                                                                                                                                           <  5:478280/67‑2 (MQ=38)
                                                                                                                                                                                                                       tAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAGTTTTTACCCAga                                                                                                                                                                                                                                                           >  6:478280/1‑66 (MQ=38)
                                                                                                                                                                                                                        atggacCACGGTCGGGTGGTGCGCGGATCTACACGACTAACGTATCAGTAAAAACTATAACCCGGCAAACCCGCGTGCGGTCGGGCGGCGTCGATGAAGGTCATGGCGAGTTTAAATAACCCGGCGGCAGCGGGTCATAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGGTGGTTGACCAGACAAAACTCCGCGCGCCTTTAATGTATCGGGCAGACTCAAACAGATTCAATTTGAAGGCGtt                                          <  4:94529/271‑1 (MQ=255)
                                                                                                                                                                                                                          caacCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaaga                                                                                                                                                                                                                            <  6:1078410/95‑1 (MQ=255)
                                                                                                                                                                                                                          caacCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCgaaga                                                                                                                                                                                                                            >  5:1078410/1‑95 (MQ=255)
                                                                                                                                                                                                                          caacCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTaa                                                                                                                                                                                                          <  5:53627/113‑1 (MQ=255)
                                                                                                                                                                                                                          caacCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTaa                                                                                                                                                                                                          >  6:53627/1‑113 (MQ=255)
                                                                                                                                                                                                                               aCGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAt                                                                                                                                                                                                               >  8:928439/1‑103 (MQ=255)
                                                                                                                                                                                                                               aCGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGAt                                                                                                                                                                                                               <  7:928439/103‑1 (MQ=255)
                                                                                                                                                                                                                                             cgcgGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTAcgc                                                                                                                                            >  4:347102/1‑156 (MQ=255)
                                                                                                                                                                                                                                             cgcgGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGCAGACTCAAACAGATTCAAt                                                     >  4:465694/1‑243 (MQ=255)
                                                                                                                                                                                                                                               cgGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGCAGACTCAAACAGATTCAATTTGAAGGCGTTATGGATCTTTCGGCGGCGTAAt                    >  2:493120/1‑274 (MQ=255)
                                                                                                                                                                                                                                                 gATCTAAACGCATCACCTCTGGGTAAAAACTTTTCCCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCAcg                                                                                           <  1:75028/201‑1 (MQ=255)
                                                                                                                                                                                                                                                 gATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCg                                                                                                                                        <  2:347103/156‑1 (MQ=255)
                                                                                                                                                                                                                                                 gATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGCAGACTCAAACAGATTCAATTTGa                                                 <  2:465695/243‑1 (MQ=255)
                                                                                                                                                                                                                                                        aCGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCcca                                                                                                                                                                      >  7:72878/1‑119 (MQ=255)
                                                                                                                                                                                                                                                        aCGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCcca                                                                                                                                                                      <  8:72878/119‑1 (MQ=255)
                                                                                                                                                                                                                                                            aTCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGCCCAGACAAAACTCCACGCGCCTTTAATGTATCGGGCAGACTCAAACAGATTCCATTTGACGGCGTTATGGATATTTCGGCGGCGTAATTCCCGCCCTGCaac      >  4:232022/1‑274 (MQ=255)
                                                                                                                                                                                                                                                            aTCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGCAGACTCAAACAGATTCAAttt                                                   >  2:465140/1‑230 (MQ=255)
                                                                                                                                                                                                                                                            aTCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGCAGACTCAAACAGATTCAATTTGAAGGCGTTATGGATNTTTCGNNGGCGTAATTAGCGGCCCGCaaa      >  3:16895/1‑275 (MQ=255)
                                                                                                                                                                                                                                                              cACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGCAGACTCAAACAGATTCAATTTGAAGGCGTTATGGATATTTCGGCGGCGTAATTCGCGCGCGGCAAac     <  2:278181/274‑1 (MQ=255)
                                                                                                                                                                                                                                                              cACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGCAGACTCAAACAGATTCAATTTGAAGGCGTTATGGATATTTCGGCGGCGTAATTCGCGCGCGGCAAac     >  1:7281/1‑274 (MQ=255)
                                                                                                                                                                                                                                                                ccTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTcc                                                                                                                                                                                    >  5:469583/1‑97 (MQ=255)
                                                                                                                                                                                                                                                                ccTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTcc                                                                                                                                                                                    <  6:469583/97‑1 (MQ=255)
                                                                                                                                                                                                                                                                 ctctGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGCAGACTCAAACAGATTCAATTTGAAGGCGTTATGGATATTTCGGCGGCGTAATTCGCGCGCGGCAAaccac  <  2:343952/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tctGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCggg                                                                                                                                                                                       <  5:403268/92‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tctGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCggg                                                                                                                                                                                       >  6:403268/1‑92 (MQ=255)
                                                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                                     
TACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAAATCAGGTCATACGCTTCGAGATACTTAACGCCAAACACCAGCGAAATGAGCGGCTTGCCGACGATCAACACCGCGAGCGCCACCAGAATACCGATTCCGCCCGCCAGTAAACCGGACTTCACACCTAACAACCACGGTCTGGTGGTGCGCGGATCTAAACGCATCACCTCTGGGTAAAAACTTTTACCCAGCAAACCCGCTGGCGTTCCGGCGGCGTCGAAGAATGTCATGGCGATTTTAAATAACCCGGCGGCAGCGGGTCCTAACACGATCCCCACCAACACTGTGCTGCACGAGTTACGCGCCGACCAGATGGAGTGGGCAATGTTGGTTGACCAGACAAAACTCCACGCGCCTTTAATGTATCGGGCAGACTCAAACAGATTCAATTTGAAGGCGTTATGGATATTTCGGCGGCGTAATTCGCGCGCGGCAAACCAC  >  NC_000913/3131256‑3131785

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