Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,707,947 C→A intergenic (‑242/+669) dppA ← / ← proK dipeptide ABC transporter periplasmic binding protein/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,707,9470CA97.6% 132.3 / ‑6.4 41intergenic (‑242/+669)dppA/proKdipeptide ABC transporter periplasmic binding protein/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (18/22);  total (19/22)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

GCATTATTCTGCTCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGCCAGGGCCAGCCATAAGTAACGCAAAAAAAGATTTAAGCAAATAT  >  NC_000913/3707702‑3708187
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
gCATTATTCTGCTCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGAGCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTCCACCAACGGTGAGAATCCacaaacaac                                                                                                                                                                                                                                             >  1:268134/1‑250 (MQ=255)
              cAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACtt                                                                                                                                                                                                                               >  1:442566/1‑251 (MQ=255)
                   gtgATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAAc                                                                                                                                                                                                                          <  2:139003/251‑1 (MQ=255)
                   gtgATGTTTGTTGTTGTATCCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCGCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCaaa                                                                                                                                                                                                                           <  1:516479/250‑1 (MQ=255)
                    tgATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGAGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAAc                                                                                                                                                                                                                          >  1:90525/1‑250 (MQ=255)
                        gtttgttGTTTTTACCCGTTGCAGTGGGTTGTCGCTTCCTGGTCTTTTTTGTTTTTTTCCCTGTCGGGGCAACGTTATGAGGTGGGGATGACGTACTAATTAACATCAGTGTAACAGTGCGGAGACTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGATGCAAAAAACAAAGATTCAAAAAACTTCaatctgct                                                                                                                                                                                                                      <  2:376633/250‑7 (MQ=255)
                         tttgttGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTAt                                                                                                                                                                                                                    >  1:440258/1‑251 (MQ=255)
                          ttgttgTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCACCGTTATGAGGTGGGGATGCCGTACTAATTCACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAACACTTCAAACAGCTAtt                                                                                                                                                                                                                   >  2:243076/1‑251 (MQ=255)
                                     aCCCTTTGCATTGGGTTGTCGCTGCCTGACCTTTTTTGGTTTTTTCCCGGTCGGGTCAACGTTATGAGTTGGGTATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACTGGGAGAATCCACAAACAAAGAGTAAAAAAACTTCAAACAGCTATTTTCAGCAgaga                                                                                                                                                                                                        <  2:552131/251‑1 (MQ=255)
                                      cccTTTGCAGTTGGTTGTCGCTGCCTTACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGaa                                                                                                                                                                                                       <  1:432340/251‑1 (MQ=255)
                                        cTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCAGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATAAATGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAAAAACAATTACAAAAACTTCAAACAGCTCTGTGCAGCAGCGAAAt                                                                                                                                                                                                     >  1:542514/1‑251 (MQ=255)
                                         tttGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGGCAACGTTATGAGGTGGGGCTGCCGTACTAATTAACATCGGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAACAACTTCACACAGCTATTTGCAGCAGAGAAAtt                                                                                                                                                                                                    >  2:211280/1‑251 (MQ=255)
                                             cAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTtgtg                                                                                                                                                                                                >  2:484081/1‑251 (MQ=255)
                                              aGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAAAATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCCACCGGGATTTACACCAACGGTGAGAATCAACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGc                                                                                                                                                                                               >  1:234853/1‑251 (MQ=255)
                                                    ttGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTcc                                                                                                                                                                                         >  1:27218/1‑251 (MQ=255)
                                                            gtcTGACGTTTGTTGTTTTTTGCCCAGTAGGGTCAACTTTAGGAGGTAGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGGTACCGACAATCTGTCATCAGAATGTGAAAACTTCAATACAGCCAACCGGGATGTACTCCAACGGTGAGAAGCGACAAATTAAGATTAATAAAACTGCAAACAGCTATGGGCAGCAGGGATTTTTGTGCTACTCCAACATGAc                                                                                                                                                                                 <  2:521837/249‑1 (MQ=255)
                                                            gCCTGGCCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATg                                                                                                                                                                                   <  1:155302/249‑1 (MQ=255)
                                                            gCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGGACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGa                                                                                                                                                                                  <  1:85622/250‑1 (MQ=255)
                                                                     tttttgttttttgCCCGGACGGGTCAACGTTATGAGGTGGGGTTGCCGTACTAATTAACATCAGTGAACCAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGAACAGAACAAt                                                                                                                                                                        >  2:556386/1‑251 (MQ=255)
                                                                                   ccGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTa                                                                                                                                                          <  1:386963/251‑1 (MQ=255)
                                                                                   ccGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTa                                                                                                                                                          <  1:225443/251‑1 (MQ=255)
                                                                                    cGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATAttt                                                                                                                                                           >  2:543308/1‑249 (MQ=255)
                                                                                         gggTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACa                                                                                                                                                    <  1:469274/251‑1 (MQ=255)
                                                                                                                         aaTTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTGTTATGCATTTAATGa                                                                                                                      <  2:234853/249‑1 (MQ=255)
                                                                                                                          aTTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGaaa                                                                                                                    >  1:455585/1‑250 (MQ=255)
                                                                                                                                     tgtATCAGTGCGGGTCCTCCACAAAATGCCCCTCGTGACCTACAATCTGACAACATAATGTGAAAACGTCAATACAGCCAACCGGGATTTACTCCAACGGTGAGAATCCACAAACAAATATTAAAAAATCTTCAAACAGCTATTTTCAGCAGAGAAATTTATGCGGCTTCAACATGACGAGAACAATCATCTTAATATTTAGTAGCATTTGGTGATTAAGGTTTAAGCGTTTAATGAAAATAATCTGAgg                                                                                                         <  2:299312/250‑1 (MQ=255)
                                                                                                                                       tAACAATGCGGATCCTCCATAAAATGCACCTCGTGACCTACAATCTGTCAACAGAATGAGAAAACGTCAATACAGCCAAACGGGATTTACACCAACGGTGAGAATCCACAAACAACGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGCAAACAATCTGAggca                                                                                                       >  2:523453/1‑248 (MQ=255)
                                                                                                                                           aGTGCGGATCCTCCATAAAAAGCCGCTCGTGCCCTACAATATGTGCAAAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCCGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTTATGAAAAAAATCTGCGGAAAAggg                                                                                                  >  2:494282/1‑250 (MQ=255)
                                                                                                                                                 gATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTg                                                                                                                                                <  1:442041/199‑1 (MQ=255)
                                                                                                                                                 gATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTg                                                                                                                                                >  2:442041/1‑199 (MQ=255)
                                                                                                                                                                       gTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCa                                                                                                                                                                      >  1:1856/1‑155 (MQ=255)
                                                                                                                                                                       gTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCa                                                                                                                                                                      <  2:1856/155‑1 (MQ=255)
                                                                                                                                                                                              aaTGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAAt                                                <  2:440258/251‑1 (MQ=255)
                                                                                                                                                                                                         gTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGCCAGGgccag                                     <  2:442566/251‑1 (MQ=255)
                                                                                                                                                                                                         gTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCAGTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGCCAGGgccag                                     <  1:556386/251‑1 (MQ=255)
                                                                                                                                                                                                             aTACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACAGTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATTCCAGGGCCAGCCAt                                 <  2:398038/251‑1 (MQ=255)
                                                                                                                                                                                                               aCAGCCAACCGGGATTTACACCAACGGTGAGAATCCACAAACAAATATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCATAACAATCAGCTTAATATTTAGCATCTTTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGATGAAAAGGTTGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGCCAGGGCCAGCCATaa                               <  1:252282/251‑1 (MQ=255)
                                                                                                                                                                                                                       ccGGGATTTACCCCAACGGTGAGAATCCACAAACAAATATTAAAAAAACTTCAAACAGCTATTTGGAGCAGAGAAATTTGTGCTACTCCAACATTACCAGAACAATCAGCTTAATAGTTAGCAACTTTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAAGCTGAGTAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCAAAAACGCACAGCGGAAAATGCCAGGGCCAGCCATAAGTAACGCa                       <  2:475202/251‑1 (MQ=255)
                                                                                                                                                                                                                        cGGGATTTAAACCGACGGCGAGAATCCACAATCAACGATTAAAAGAACTTCAAACAACTATTTGCAGCCGAGAAATTTGTGCTACTCCATCATGACCAGAACAATGAGCTTAATATTTAGAAACATTTGGCGATTAAGTTTTATGCATTTAATGACAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTACAACGCACAGCGGAAAATGCAAGGGCCAGCCATAAGTAACGCa                       >  2:32557/1‑250 (MQ=255)
                                                                                                                                                                                                                                   ccAACGGTGAGAATCCACAAACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGCCAGGGCCAGCAATAAGTAACGCAAAAAAAGATTTa           <  1:257320/251‑1 (MQ=255)
                                                                                                                                                                                                                                          tGAGAATCCACAAACACAGATTAAAAAAACTTCAAACCGCTCTTTGCAGCAGAGAAATTTGTGCTACTCAACCATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGCCAGGGCCAGCCATAAGTAACGCAAAAAAAGATTTAAGCAAat    >  2:576051/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gaATCCACAAACAAAGATTAAAAACACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGCCAGGGCCAGCCATAAGTAACGCAAAAAAAGATTTAAGCAAatat  >  2:541906/1‑250 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
GCATTATTCTGCTCCAATTGTGATGTTTGTTGTTTTAACCCTTTGCAGTGGGTTGTCGCTGCCTGACCTTTTTTGTTTTTTGCCCGGTCGGGTCAACGTTATGAGGTGGGGATGCCGTACTAATTAACATCAGTGTAACAGTGCGGATCCTCCATAAAATGCCCCTCGTGACCTACAATCTGTCAACAGAATGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAAAGGTGGATATCTGAGCGATTAATACCATGAACGCTAAAACGCACAGCGGAAAATGCCAGGGCCAGCCATAAGTAACGCAAAAAAAGATTTAAGCAAATAT  >  NC_000913/3707702‑3708187

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