Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,652,331 T→C intergenic (+2346/+397) ydfD → / ← ynfP lysis protein/protein YnfP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,652,3310TC98.4% 236.6 / ‑6.2 62intergenic (+2346/+397)ydfD/ynfPlysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (29/32);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGC  >  NC_000913/1652092‑1652575
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
gtgtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTa                                                                                                                                                                                                                                            <  2:486770/250‑1 (MQ=255)
  gtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATgg                                                                                                                                                                                                                                         <  2:366225/251‑1 (MQ=255)
  gtTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATg                                                                                                                                                                                                                                          <  2:1122051/250‑1 (MQ=255)
   tttCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACC‑TTAACGAGACTGAGCAAAGAGCATATCATTTTTGTTCACCTCCGAGAGTATGGc                                                                                                                                                                                                                                        >  2:136538/1‑251 (MQ=255)
     tCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGca                                                                                                                                                                                                                                       <  1:6718/250‑1 (MQ=255)
            ataCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaacaga                                                                                                                                                                                                                                >  1:513633/1‑250 (MQ=255)
            ataCTGGAGCTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTGATTCATCTTATTAAGCCACCAATCG‑ACGCTCTCCGGAGGAAGATGACA‑TTTACGAGACTGAGCAAAAAGCATATCATTTATGTTCACCTCAGAGAGTATGGcagaacaga                                                                                                                                                                                                                                >  1:567867/1‑250 (MQ=255)
              aCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATTTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGGCAGATGACA‑TTAACGAGACTGAGCAAAAAGCATATCATTGATGTTCAACTCAGAGAGGATGGCAGAACAGaaa                                                                                                                                                                                                                              >  1:812336/1‑250 (MQ=255)
              aCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTACGAGAGTATGGCagaaaagaaa                                                                                                                                                                                                                              >  1:443285/1‑250 (MQ=255)
                        gctcTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAGCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGc                                                                                                                                                                                                                   <  1:980023/248‑1 (MQ=255)
                                 aGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGGCATTCATCTTATTAAGCCAGCAATCG‑AAGCTCTCCGGAGTCAGATGACA‑TTAACGGGACTGAGCAAAGAGCATATCATTGATGTTCACCTCTGAGAGGATAGCAGAACAGACACCAAAAAATTCACCTTTGtt                                                                                                                                                                                                          >  1:987952/1‑251 (MQ=255)
                                   cTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTTATCTTATTAAGCCAGCAATCG‑ACGCTCTCCCGAGTCAGATGACA‑TTAAAGAGGATGAGCAAAGAGCATAACATTGATGTTCACCTCAGAGAGTGAGGCAGAACAGAAAAACAAAAAAGCAACttttgttt                                                                                                                                                                                                        >  1:268471/1‑246 (MQ=255)
                                          cACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCCGGGCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAAAGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCt                                                                                                                                                                                                 >  1:220129/1‑251 (MQ=255)
                                            cTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGGCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCaaactt                                                                                                                                                                                                >  2:948842/1‑246 (MQ=255)
                                              tGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGa                                                                                                                                                                                               <  2:1046728/249‑1 (MQ=255)
                                                gATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAgtg                                                                                                                                                                                           <  2:262904/251‑1 (MQ=255)
                                                 aTGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAgtgt                                                                                                                                                                                          <  1:136538/251‑1 (MQ=255)
                                                             aGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaa                                                                                                                                                                              <  2:1024175/251‑1 (MQ=255)
                                                             aGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCCGAGTCAGATGACA‑TTAACAAGACTGAGCAAAGAGCCTATCATTGATGTTCACCCCAGAGAGTATAGCAGAAAAGAAAAACAAAAATTCACCTTTGTTTTTTAACCCGAAAGGTCAACGAGAGTa                                                                                                                                                                               >  1:223934/1‑250 (MQ=255)
                                                              gATCTGAACAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGACAACAAAAGGGGCATTCATCTTTTTAAGGCAACAATCGCACCCT‑TCCCGGGTCAGGTGAAAATTAA‑GGGACTTAGCAAAGTGCCAATCCATGTTGTTCTCCCCCGAGAGTATGGCAGCACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGGAGTGTCAGCgaggggaa                                                                                                                                                                              >  2:769779/1‑245 (MQ=255)
                                                                          ggAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCCATCT‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGAATGAGCCAAGAGCATATCATTGGTGTTCACCTCAGAGAGGATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGAGAGAGGAAAAAATTATGGTg                                                                                                                                                                  >  1:1150896/1‑250 (MQ=255)
                                                                          gaaaCAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTg                                                                                                                                                                  <  1:412950/248‑1 (MQ=255)
                                                                                           gagaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCAACTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGGATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTAACGGTGAcc                                                                                                                                                >  1:791134/1‑250 (MQ=255)
                                                                                             gaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACTCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCa                                                                                                                                              <  2:220129/251‑1 (MQ=255)
                                                                                                      cAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCCATCGGGTCGGTCATCTTATTAAAACAGCAAACG‑ACGCCCACCGGGGTCAGATGACA‑TTAAAGAGACTGAGCAAAGAGCATATCATTTATGTTCACCTCAGAGAGGATGGCAGAAAAAAAAAAAAAAAATGACCCTTTGTTTTTCAACCCGAAGTGTCAGCGAGAGGAAAAAATTATGGTGACCATTTTACCGTTTAATCAATaaagc                                                                                                                                       >  2:236218/1‑246 (MQ=255)
                                                                                                               gTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGAAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGAc                                                                                                                            >  1:547484/1‑251 (MQ=255)
                                                                                                                             aaaaCACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACg                                                                                                              <  1:948842/251‑1 (MQ=255)
                                                                                                                               aaCACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTg                                                                                                            <  2:1165210/251‑1 (MQ=255)
                                                                                                                                   caGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTa                                                                                                                                                                               >  2:355456/1‑180 (MQ=255)
                                                                                                                                   caGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTa                                                                                                                                                                               <  1:355456/180‑1 (MQ=255)
                                                                                                                                              cAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATc                                                                                             <  2:443285/251‑1 (MQ=255)
                                                                                                                                              cAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATc                                                                                             <  2:547484/251‑1 (MQ=255)
                                                                                                                                                  cGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGaaa                                                                                         <  1:961271/251‑1 (MQ=255)
                                                                                                                                                     gTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATc                                                                                      >  1:753914/1‑251 (MQ=255)
                                                                                                                                                                tattaAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAgg                                                                              <  2:1150896/248‑2 (MQ=255)
                                                                                                                                                                         aGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACtt                                                                  <  1:731665/251‑1 (MQ=255)
                                                                                                                                                                             aTCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGc                                                              >  1:682409/1‑251 (MQ=255)
                                                                                                                                                                              tCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCAtt                                                                                                                                                            <  2:482132/156‑1 (MQ=255)
                                                                                                                                                                               tcgccGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCAtt                                                                                                                                                            >  1:482132/5‑156 (MQ=255)
                                                                                                                                                                               cG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAc                                                                                                                                                                <  1:685063/151‑1 (MQ=255)
                                                                                                                                                                               cG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGAc                                                                                                                                                                >  2:685063/1‑151 (MQ=255)
                                                                                                                                                                                g‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctgc                                                          >  2:961769/1‑251 (MQ=255)
                                                                                                                                                                                         ccGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTc                                                     <  2:655901/249‑1 (MQ=255)
                                                                                                                                                                                          cGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAAGCATATCAGTCGAGACATACTTATGCCGGCTGGTCGCt                                                  <  2:139987/251‑1 (MQ=255)
                                                                                                                                                                                                cAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGc                                            <  2:268471/251‑1 (MQ=255)
                                                                                                                                                                                                cAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAA‑CAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCt                                           <  1:961769/251‑1 (MQ=255)
                                                                                                                                                                                                 aGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCt                                           >  1:1025485/1‑251 (MQ=255)
                                                                                                                                                                                                    tGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTggg                                        >  1:265210/1‑250 (MQ=255)
                                                                                                                                                                                                          ttAACGAGCCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACCGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGCCGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGATGGTCGCTGCCAGCTGGTGCTa                                    >  2:420191/1‑250 (MQ=255)
                                                                                                                                                                                                            aaCGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCTTCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAAcc                                 <  1:668143/251‑1 (MQ=255)
                                                                                                                                                                                                               gagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACa                                              >  2:542574/1‑235 (MQ=255)
                                                                                                                                                                                                               gagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACa                                              <  1:542574/235‑1 (MQ=255)
                                                                                                                                                                                                                 aacggAGCAAAGAGCATATCATTTATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAAGTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCGAACCCGGCa                            <  2:503920/247‑1 (MQ=255)
                                                                                                                                                                                                                      aGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTtat                       >  1:380277/1‑251 (MQ=255)
                                                                                                                                                                                                                            gagCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCaaa                 <  2:380277/251‑1 (MQ=255)
                                                                                                                                                                                                                                  atCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGAt           <  1:1074522/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATgg         <  2:682409/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aTTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATggg        >  2:348374/1‑251 (MQ=255)
                                                                                                                                                                                                                                      ttGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCgg                                                                                                        <  1:671740/154‑1 (MQ=255)
                                                                                                                                                                                                                                      ttGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCgg                                                                                                        >  2:671740/1‑154 (MQ=255)
                                                                                                                                                                                                                                       tggtcatCACTTCAGAGAGTATCGCAGAACAGACAAATAACAATGCACCTTTGTTTTTCAACCTGAAGAGAAAGCGAGAGTAAAAAATTATGGTGACTATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCTAAAATCATAACAGCGGTGACATACTTATGCTTGCTGGTCGCTGACAGCTGGTGCTAATCCGGCCTTTATAGCACAACGGATGGGcc      >  2:1182111/7‑251 (MQ=255)
                                                                                                                                                                                                                                             tCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGc  >  2:196793/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                    
GTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCG‑ACGCTCTCCGGAGTCAGATGACA‑TTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCATGC  >  NC_000913/1652092‑1652575

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