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,3310TC97.6% 134.4 / ‑6.1 41intergenic (+2346/+397)ydfD/ynfPlysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (15/25);  total (16/25)
Fisher's exact test for biased strand distribution p-value = 3.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAT  >  NC_000913/1652087‑1652573
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
aTGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagag                                                                                                                                                                                                                                              <  1:646693/251‑1 (MQ=255)
 tggggtgtTTCTGTTTATACTGGAGTTCGCCCTGGGGAGCTGTGTGCACTGGGATGGGAGGACATAGATCAAAAAAATGGAACAATGATGATCAGGAGAAATGTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTAGTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCagag                                                                                                                                                                                                                                              <  1:449897/246‑1 (MQ=255)
                      ggAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGGCATTCATCTTATTTAGCCAGCAATCGACGCTCCCCGGGGTAAGATGACAATAAAGAGAATGAGCAAAAAGCATTTCATTGATGTTCCCCTCAGAGAGGATGGGAAAACAGAAAAAc                                                                                                                                                                                                                         >  1:286709/1‑250 (MQ=255)
                             gCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGCAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATCCATCTTATTAAGACAGCAATACACCCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAAAGCATATCCTTTATGATCCACTCAGAGAGAATGGCAGAACAGAAAAACAAAAATGc                                                                                                                                                                                                                 >  1:587975/1‑251 (MQ=255)
                                   gTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCttt                                                                                                                                                                                                           <  1:156072/251‑1 (MQ=255)
                                   gTGAGCTGTGTACACTTGTATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCATCAATCGACGCTCTCCGGATTCAGATGACATTAACGAGACTGAGCAAAAAGCATCTCATTGATGTTCACCTCAGATATTATGGCAGAACAGAAAAACAAAAATGCACCtgt                                                                                                                                                                                                           <  1:156077/251‑3 (MQ=255)
                                       gcagtgtgCACTTGGATGGGATGACATAGATCTGAAAAATGGAACAATGATGATCAGGATAAATTTAGCAAAAGACTGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttt                                                                                                                                                                                                       <  1:494818/248‑1 (MQ=255)
                                         tgtgtgCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGTACCAAACGGGTCATTCATCTTATTAAGCCAGCAAACGACGCCCTCCGGAGTCAGATGACATTAAAAAGACTGAGCCAAACGCATATCATTGATGTTTACCTCAGAGAGTATTGCAAAACAAAAAAAAAAAAATGCACCCtttttttt                                                                                                                                                                                                     >  2:133076/1‑251 (MQ=255)
                                                aCTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAAAGCATATCATTGATGTTCACCTCCAAGAGGATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTg                                                                                                                                                                                              >  1:602242/1‑251 (MQ=255)
                                                        gggAGGACATAGATCTGAAAAATGTAACAAAGATGATCAGGAGAAATTTAGCAAAAGTCCGTTTCACGGTATCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGTCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAAACTGAAGTGTc                                                                                                                                                                                       <  2:438376/250‑1 (MQ=255)
                                                         ggaggaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCCGGTCCTTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGGGTCAGATGACATTAAAGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGGAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTTAAGTGTc                                                                                                                                                                                       >  1:566214/1‑249 (MQ=255)
                                                                                gAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAAACGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTGCAACCTGAAGTGTCAGCGATAGTAAAAAATTATGGTGAc                                                                                                                                                              <  1:593755/251‑1 (MQ=255)
                                                                                       gatgatCAGGAGAAATATAGCAAAAGACCGTTTCACGGTATCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCCCCCCGGAGTCAGCTGACATTACCGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAAt                                                                                                                                                                                                                   <  1:635506/191‑1 (MQ=255)
                                                                                            tCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaacaga                                                                                                                                                                                                                              <  1:379975/175‑1 (MQ=255)
                                                                                            tCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGGCATTCATCTTTTTAAGCCCACAATCGACGCTCTCCCGAGGCAGATGACCATAACGAGACTGAGCAAAGAGCCTATCATTGATGTTCACCTCAGAGAGTATGGgagaacaga                                                                                                                                                                                                                              >  2:379975/1‑175 (MQ=255)
                                                                                                    aaTTTATCAAAAGTCCCTATCACTGTACCAAACATACAAGCGTGAACCATTTGGGTCATTTATCTTATTAAGCCAGTAATCGACGCTCTCCGGTGTCAGATGACATTAACGAGACTGAGCAAAGAGCATAACATTGATGTTCACCTCAGAGAGTATGTCAGAACAGAAAAACAAAAATGAACCTTTTTTTTTCAACCTGAATTTTCAGCGTGAGTAAAAAATTATGGTGTCTATGTTATCGGTGACTCAAt                                                                                                                                          <  2:608028/251‑1 (MQ=255)
                                                                                                         aGCAAAAGACCTTTTCACTGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAAGCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTTATGTTCTCCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAgg                                                                                                                                      <  2:183933/250‑1 (MQ=255)
                                                                                                                 aCCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGc                                                                                                                                                                                    <  1:692714/196‑1 (MQ=255)
                                                                                                                            gtatcAAAAACACAGGCGGGAATCAATCGGGTCATTCATCTTATGAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGGTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGTAGATTTGGGACGCAGCGAt                                                                                                                  <  2:370774/247‑1 (MQ=255)
                                                                                                                                cAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCGGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTaaa                                                                                                                                                                           <  1:119581/190‑1 (MQ=255)
                                                                                                                                                                        taAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCgaga                                                                       <  1:595590/250‑1 (MQ=255)
                                                                                                                                                                               gCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTa                                                               <  1:487138/251‑1 (MQ=255)
                                                                                                                                                                                 aagcGACGCTCTCCGGAGTCAGAGGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTTAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGTCGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATg                                                             <  2:203277/248‑1 (MQ=255)
                                                                                                                                                                                    cGACGCTCTCCGGCGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGCCCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAATTCATATCACTCGAGCCATACTTATGctt                                                          >  2:26655/1‑249 (MQ=255)
                                                                                                                                                                                    cGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCct                                                          <  1:75656/251‑1 (MQ=255)
                                                                                                                                                                                     gACGATCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATAATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGATTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACAATTTTACCGTTGCCTCAATAAGACAGATGTGGGACGCAGCGATAAAACGTGCAGGACTCCGCCATCGAAAATCATATAAGTCGAGACATAATTATGCctg                                                         >  1:640800/1‑251 (MQ=255)
                                                                                                                                                                                        gCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCctgctg                                                      <  1:452497/251‑1 (MQ=255)
                                                                                                                                                                                         ctctCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAACAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCAGGACTCCGCCATCGAACATCATATCAGTCGAGACATCCTTATGCCTGCTgg                                                     >  2:696955/1‑251 (MQ=255)
                                                                                                                                                                                             ccGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTc                                                   <  1:156611/249‑1 (MQ=255)
                                                                                                                                                                                              cGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCt                                                <  1:354481/251‑1 (MQ=255)
                                                                                                                                                                                                     aGATGACATTAACGAGACTGAGCAAAGAGCATATCATGGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGCAGTGTCAGCGAGAGCAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACATGCCGGACTCCGCCATCGAAAATCATATCAGGCGAGACATACGTATGCCTGCTGGTCggtgcaagct                                         >  2:151739/1‑242 (MQ=255)
                                                                                                                                                                                                       aTGACATTAACGAGACTGAGCCAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGAAGATGTGGGACGCAGCGATAAAACGTGCCGGCCTCCGCCATAGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGaaacctgg                                       >  2:14546/1‑244 (MQ=255)
                                                                                                                                                                                                          aCATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGc                                    <  1:418397/251‑1 (MQ=255)
                                                                                                                                                                                                                  gagaCTGAGCAAAGAGCATATCCTTGCTGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCt                                         >  2:262159/1‑238 (MQ=255)
                                                                                                                                                                                                                  gagaCTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCt                                         <  1:262159/238‑1 (MQ=255)
                                                                                                                                                                                                                         aGCAAAGAGCATCTCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGAAAAACATTATGGTGACCATTTTACCGTTGACTCAATAAAGCAGATTTGGGAAGCAGCGATAAAAAGTGGCGGACTCCGACCTCGAAAAACATATCAGTAGAGACATACTTATGCCTGATGGTCACGTACAGCAGGAGCAAACCCAGCATTtat                     >  2:479830/1‑251 (MQ=255)
                                                                                                                                                                                                                            aacgagCATATCATTGATGTTCCCCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCCCCTTTGTTTTTCAACCTAAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGAATCACTAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGCAAATCATATCAGACTAGACATACTTATGCCTGCTGGTAGCTAACATCTGGTGCTTAACCGGCCTTtatacc                  >  2:311562/4‑249 (MQ=255)
                                                                                                                                                                                                                                   atatCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGCCTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAAcacg           >  2:80506/1‑248 (MQ=255)
                                                                                                                                                                                                                                    tatCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTCCCGTTGACTCAATAAGGCAGATGTAGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGGAAATCATCTCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAgg          >  2:546674/1‑250 (MQ=255)
                                                                                                                                                                                                                                     atCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATTCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGa          <  1:352687/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aTGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAt  <  1:26655/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
ATGGTGTGTTTCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAT  >  NC_000913/1652087‑1652573

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