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% 139.3 / ‑4.8 42intergenic (+2346/+397)ydfD/ynfPlysis protein/protein YnfP
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (20/21);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

TCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAT  >  NC_000913/1652097‑1652573
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
tCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGCATTCATCTTATTAAGCCAGCAATCGACGCCCTCCGGAGTCAGATGAAATTAAAGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGc                                                                                                                                                                                                                                      >  2:173530/1‑249 (MQ=255)
   gTTTATACTGGAGTTCGCGCTGGTGTGCTGTGTGCACTTTGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAACTTTAGCATAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATCAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGAAAAGAGCATATCATTGAGGTTCCCCTCAGAGAGTATGGcaga                                                                                                                                                                                                                                   <  2:539136/249‑1 (MQ=255)
   gTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcag                                                                                                                                                                                                                                    <  1:367531/248‑1 (MQ=255)
       ataCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaacaga                                                                                                                                                                                                                              <  2:433847/250‑1 (MQ=255)
       ataCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGAGCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGcagaacaga                                                                                                                                                                                                                              <  1:411865/250‑1 (MQ=255)
         aCTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGaaaa                                                                                                                                                                                                                           <  1:261872/251‑1 (MQ=255)
            ggAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAAc                                                                                                                                                                                                                         <  1:428070/250‑1 (MQ=255)
                ttCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGGAATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGAATGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaaa                                                                                                                                                                                                                     >  1:516996/1‑250 (MQ=255)
                ttCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATAGGGTCATTCATCTTATTAAACCAGCAAACGACGATCTCCGGGGGCAGATGACATTAACGAGAACGAGCAAAGAGCATAACAATGATGGTCACCCCCGAGCGTATGGCATAACAGAAAAACcaaa                                                                                                                                                                                                                     >  2:332257/1‑250 (MQ=255)
                      cTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATAGACGCTCTCCGGAGTCAGATGACATTAAAGAGACTGAGCAAAGAGCATATCATTGATGTTTCACTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGcccc                                                                                                                                                                                                              >  1:50033/1‑248 (MQ=255)
                             gCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGttt                                                                                                                                                                                                       <  1:466035/251‑1 (MQ=255)
                                    gCACTTGGATGGGAGGACATAGCTCTGACAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACCGGCCATTCATCTTTTTAAGCCAGCAATCGAAGCTCCCCGGGGTCAGATGACCTTAACGAGACTGAGCAAAGAGCATATCCTTGATGTTTACCTCCGAGAGTATGGCAGAACAGAAAAACAAAAATCCCCCTTTTTTTTTCaaac                                                                                                                                                                                                >  2:617327/1‑249 (MQ=255)
                                      aCTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTg                                                                                                                                                                                              <  2:7265/251‑1 (MQ=255)
                                              gggAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCa                                                                                                                                                                                      <  1:645959/251‑1 (MQ=255)
                                                gaggaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTAACGGTACCAAAAACACAGGCGGGAAACAATCGGGTCATTTATCTTATTAAACCAGGAATCTACGATCTTCGGAGTCAGAAGAAATTAAAGAGAATGAGCAAAAAGCATATCATTGAAGTTCATCTCAGAGAGGAAGGCTGAACAGGAAAAAAAAAATGCAACTTTGtttgtacaacggaagggtgagc                                                                                                                                                                                    >  1:234683/1‑232 (MQ=255)
                                                  ggaCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACAGTTTCACGGTACCAAAAACACAGGCGGGAAACAATACGGGCATTCCTCTTATTAAGACAGCAATCGACGCTCTCCGGAGTAAAATGACATTAAAGAGACTGAGCAAAGAGCAAATCATTGAAGTTCACCTCAGAGAGAATGGGAGAACAGAAAAACAAAAAAGCAACCTTGTTTTTCCACATGAAATGTCAGCg                                                                                                                                                                                   >  1:371507/1‑250 (MQ=255)
                                                    aCATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCgaga                                                                                                                                                                                <  1:173530/251‑1 (MQ=255)
                                                                  aaTGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATgg                                                                                                                                                                  >  1:287006/1‑251 (MQ=255)
                                                                               tgatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCCCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGTTGTTCACCTCAGAGAGTATGGCAGAAAAGAAAAACAAAAATGCACCTTTTTTTTTTAACCTGAAGTGTCAGCGAGAGTAAAAAAATATTGGGACCATTTTAcc                                                                                                                                                     >  1:199786/1‑251 (MQ=255)
                                                                                gatCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAAACAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGGCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGAAGTTCACCTCAGAGAGTATGGCAGAACAGAAAAAAAAAAATGAAACTTTGTTTTTCAACCTGAAGCGTCAGCGAGAGGAAAAAAATATGGTGAACATTTTACCg                                                                                                                                                    >  2:622576/1‑251 (MQ=255)
                                                                                 atCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGt                                                                                                                                                   <  2:270669/251‑1 (MQ=255)
                                                                                       agaAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAAACGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAAAGCATATCATTGATGTTTACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAACGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGaatc                                                                                                                                             >  1:15438/1‑248 (MQ=255)
                                                                                          aaTNTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGTGAACCTATCGTGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATTATTGATGTTCACCTCAGAGAAGATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGTAGTGACAGCGAGAGTAAAAAAGTATGGTGGCCATTTGACCTTTGTCTCAAt                                                                                                                                          <  2:596088/251‑1 (MQ=255)
                                                                                            tttAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGGACCAACCGGGGCAAACACATTATTAAGACAGCAATCGCCGCCCTCCGGAGTCAGATGACCATAACCAGACTCAACAAAGAGCATATCCTTTATGTTCAACTCAGAGAGTATGGAAGAACAGAAAAACAAAAATGCACCTTTGTTTTTTAACCTGAAGGGTCAGCGAGAGCAAAAAAATATGGTGGCCATTTTACCGTTGACTCAATa                                                                                                                                         >  1:438844/1‑250 (MQ=255)
                                                                                               aGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACAAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGGAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGc                                                                                                                                     <  1:575515/251‑1 (MQ=255)
                                                                                                  aaaaGACCGTTTCACGGTACCAAAAACACAGGCGGGAAAAAATCGGGTCCTAGAACTTTTTAAACCAGAAATCGACGCCCTCCGGAGGCAGATGACATTAACGAGACTCGGCCAAAAAACTATCATTGATTTTTAACTCAGAGAGTATGGCAGAAAAGAAAAACAAAAAACCACCTTTGGTTTTCAACCTGAAGTGTCAACGAAAGGAAAAAATTATGGTGACCATTTTACCGTTTACCCAATaaagc                                                                                                                                     >  1:380336/1‑245 (MQ=255)
                                                                                                                cGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACaaaa                                                                                                                                                                                                                     <  1:342170/154‑1 (MQ=255)
                                                                                                                cGGTACCAAAAACACAGGCGGGAACCAATCGGGTCAATCATCTTATTAAGCCCGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTTATGTTCACCTCAGAGAGTATGGCAGAACAGaaaaaaaaaa                                                                                                                                                                                                                     >  2:342170/1‑154 (MQ=255)
                                                                                                                    aCCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATaa                                                                                                                <  1:615671/251‑1 (MQ=255)
                                                                                                                                        ccAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCAt                                                                                            >  1:46944/1‑251 (MQ=255)
                                                                                                                                                  tcattcatCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATc                                                                                    >  2:360977/1‑250 (MQ=255)
                                                                                                                                                           ttattaAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCg                                                                          >  1:92200/1‑251 (MQ=255)
                                                                                                                                                                      gCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTa                                                               <  1:144237/251‑1 (MQ=255)
                                                                                                                                                                           cGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGcc                                                           <  2:516996/250‑1 (MQ=255)
                                                                                                                                                                                           cAGATCACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGc                                          >  2:421778/1‑251 (MQ=255)
                                                                                                                                                                                             gATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTg                                        >  1:571393/1‑251 (MQ=255)
                                                                                                                                                                                                     gaaCGAGACTGAGCAAAGAGCATATCATTGATGTTCACCGCAGAGAGTATGGCAGAACAGAAAAACAAAAATGAACCTTTGTGTTGCAACCTGAAGTGTCAGCGAGATTAAAAAATTATGGGGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCCTCGAAAATCATATCAGGCGAGACATACTTATTCCTGCTGGTCGCTGACAGCTGGTGCTAAc                                <  1:228126/250‑1 (MQ=255)
                                                                                                                                                                                                                    aaaGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGc                  >  1:269475/1‑251 (MQ=255)
                                                                                                                                                                                                                     aaGAGCATATCATTGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGTCTCCGCCATCGAAAATCATATCATTCGAGACATACTGATGCCTGCTGGTCGCTGGCAGCTGGTGCTAACCCGGCATTTATAGCa                 <  2:362569/251‑1 (MQ=255)
                                                                                                                                                                                                                                  tGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCg                                                                                                                                                                                   >  2:260866/1‑76 (MQ=255)
                                                                                                                                                                                                                                  tGATGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCg                                                                                                                                                                                   <  1:260866/76‑1 (MQ=255)
                                                                                                                                                                                                                                    aTGTTCACCTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAt  <  1:622576/251‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
TCTGTTTATACTGGAGTTCGCCCTGGTGAGCTGTGTGCACTTGGATGGGAGGACATAGATCTGAAAAATGGAACAATGATGATCAGGAGAAATTTAGCAAAAGACCGTTTCACGGTACCAAAAACACAGGCGGGAACCAATCGGGTCATTCATCTTATTAAGCCAGCAATCGACGCTCTCCGGAGTCAGATGACATTAACGAGACTGAGCAAAGAGCATATCATTGATGTTCACTTCAGAGAGTATGGCAGAACAGAAAAACAAAAATGCACCTTTGTTTTTCAACCTGAAGTGTCAGCGAGAGTAAAAAATTATGGTGACCATTTTACCGTTGACTCAATAAGGCAGATGTGGGACGCAGCGATAAAACGTGCCGGACTCCGCCATCGAAAATCATATCAGTCGAGACATACTTATGCCTGCTGGTCGCTGACAGCTGGTGCTAACCCGGCATTTATAGCAAACCAGATGGGCCAT  >  NC_000913/1652097‑1652573

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