Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,014,019 A→G 14.1% E254G (GAA→GGA)  pqiB → paraquat‑inducible, SoxRS‑regulated MCE domain protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,014,0190AG14.1% 217.3 / 2.1 76E254G (GAA→GGA) pqiBparaquat‑inducible, SoxRS‑regulated MCE domain protein
Reads supporting (aligned to +/- strand):  ref base A (6/57);  new base G (9/3);  total (15/61)
Fisher's exact test for biased strand distribution p-value = 6.76e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGTTCCGTGGCTATCGGGTAGGTTCGGTTGAAACCAGCACCTTCGATACACAAAAACGCAATATCAGCTATCAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCAAAGTGCCATTCTTTGCGCCGAATATGCGTCAGACATTTAACGATGATTA  >  NC_000913/1013786‑1014256
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                
ttttttCCGTGGCTTCCGGGTAGGTTCGGTTGAAACCAGCCCCTTCGATACACAAAAACGCAATATCAGCTATCAACTGTTCATCAACGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGt                                                                                                                                                                                                                                              <  2:1352244/234‑1 (MQ=255)
gttGTTCCGGGGTTTTCGGGTAGGTTCGGTTGAAACCCGCCCCTTCGATCCACAAAAACGCAATATAAGCTATCAACTGTTCATCAATGCACCTTATGGCTGGCTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTc                                                                                                                                                                                                                                             <  2:164715/237‑1 (MQ=255)
atgttttaaaaaagaaggggTTGGTTCGGTTGAAACCAGCCCCTTCGATCCCAAAAAACGCAATATCAGCTATCAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTc                                                                                                                                                                                                                                             <  2:1300911/222‑1 (MQ=255)
 ttGTTCCCGGGCTTTCGGGGTGGTTTGGTTTAAACCCGCACCTTCGATCCCCAAAAAAGCAATATCAGCCATCAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTc                                                                                                                                                                                                                                             <  2:826794/237‑1 (MQ=255)
 cttttCCGTGGCTAGCGGGTAGGTTCGGTTGAAACCAGCACCTTCGATACACAAAAACGCAATATCAGCTATCAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCt                                                                                                                                                                                                                                            <  2:466523/236‑1 (MQ=255)
 cTGTTCCGTGGCTATCGGGTAGGTTCGGTTGAAACCAGCACCTTCGATACACAAAAACGCAATATCAGCTATCAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTc                                                                                                                                                                                                                                             <  1:152265/238‑1 (MQ=255)
  tGTTCCGTGGCTATCGGGTAGGTTCGGTTGAAACCAGCACCTTCGATACCCAAAAACGCAATATCAGCTATCAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCt                                                                                                                                                                                                                                            <  1:943574/238‑1 (MQ=255)
    ttCCGTGGCTATCGGGTAGGTTCGGTTGAAACCAGCACCTTCGATACACAAAAACGCAATATCAGCTATCAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGGGTTAGGTTTTTTGTGGCGGGAggggtggg                                                                                                                                                                                                                                         >  1:599638/1‑233 (MQ=255)
                            gTTGAAACCAGCACCTTCGATACACAAAAACGCAATATCAGCTATCAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCTGCTTTTATGTTCCGGGAGGGCTGTGTTTTAgggtgtctcttttacaca                                                                                                                                                                                                                  >  1:1300911/1‑223 (MQ=255)
                             tttaaaccgaccccttggtaccacaAAAACGCAAGTACAGCTATCAACTTTTCATCAATTCACCTTATGACCGACTGGTGACCAACAATGGTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGaaa                                                                                                                                                                                                              <  2:884387/218‑1 (MQ=255)
                                   ccAGCCCCTTCGATACACAAAAACGCAATTCCAGCTTCCAACTGTTCATCAATGCACCTTATGGCCGACTGGTGACCAACAAGGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCGGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCt                                                                                                                                                                                                        <  2:464614/239‑1 (MQ=255)
                                    cAGCACCTTCGATACACAAAAACGCAATATCAGCTATCAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCtt                                                                                                                                                                                                       <  1:921667/239‑1 (MQ=255)
                                         ccTTCGATACACAAAAACGCAATATCAGCTATCAACTGTTCATCAATGCACCTTATGAACGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGtt                                                                                                                                                                                                  <  1:193400/239‑1 (MQ=255)
                                           ttCGATACACAAAAACGCAATATCAGCTATCAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCGGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGtttt                                                                                                                                                                                                <  1:1213101/239‑1 (MQ=255)
                                              gATTCACAAAAACGCAATATCAGCTATCAACTGGTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTa                                                                                                                                                                                             <  1:758973/239‑1 (MQ=255)
                                                tACACAAAAACCAAATTTCAGCATACGACTGTTAAACAATGGCCCTGATGACCGGCTGGGGACCAACAATGATCGCTGCTGGAAAGAGAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCGGAGTGGGGGTGGCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCGTTCGTTTTGGatg                                                                                                                                                                                           <  2:170989/239‑1 (MQ=255)
                                                       aaaCGCAATATCAGCTATCAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGGGTCAGGTTTTATGTGCCGGAGGGGCGTGGTTTTGGGCGGCCAGTGGGCCCGAAAACAGCTTTCTGTTTGTTTGATGATCa                                                                                                                                                                                    >  1:619585/1‑239 (MQ=255)
                                                           gCAATATCAGCTATCAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGtt                                                                                                                                                                                                  <  1:526652/221‑1 (MQ=255)
                                                                  ccgCTATCACCCTTTCATCAATGCACCTTTTGGCCGTTTGGTTACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTc                                                                                                                                                                         <  2:619585/237‑1 (MQ=255)
                                                                      tATCAACGGTTCACCAATGCACCGTATGGCCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAgg                                                                                                                                                                      <  1:700824/238‑1 (MQ=255)
                                                                       ttaaaCTGTTCAACAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGAt                                                                                                                                                                    <  1:285084/236‑1 (MQ=255)
                                                                       ttCAACTGTTCATCAATGCACCTTAAGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGAt                                                                                                                                                                    <  1:1134083/238‑1 (MQ=255)
                                                                       aTCAACTGTTCATCAATGCACCTTATGACCGACTTGTGACCAACAATGTTCGCTTCTGGAAAGATGGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGAt                                                                                                                                                                    <  1:1325657/239‑1 (MQ=255)
                                                                       aTCAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGAt                                                                                                                                                                    <  1:415017/239‑1 (MQ=255)
                                                                       aTCAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGAt                                                                                                                                                                    <  1:1203093/239‑1 (MQ=255)
                                                                         cAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTATGTGCCGGGAGGGCTGGGTTTAGGGGAGCCCGTGGCACCGAAAAAAAACTTTCGTTTTTATATTGATTCGAAAAGGATTCacggttt                                                                                                                                                                  >  1:229544/1‑233 (MQ=255)
                                                                           acttgtcatcaATGGCCCGTATGACCGACTTGTGACCAACAATGTTCGCTTCTGGAAAGGTCGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGt                                                                                                                                                                <  1:1339266/234‑1 (MQ=255)
                                                                                           ccTTATGACCGACTGGTGACCAACAATGTTCGCGTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACAt                                                                                                                                                <  1:454109/239‑1 (MQ=255)
                                                                                            cTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCt                                                                                                                                                                                                                                            <  1:205923/148‑1 (MQ=255)
                                                                                            cTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAggg                                                                                                                                                                                                                                 <  1:534264/158‑1 (MQ=255)
                                                                                             ttattgcctattGGTTACCAAATCTGGTCGCTTCTGGAAAGGTAGTGGCATTGCGGTTGAGCTGGCGGCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTg                                                                                                                                              <  2:867566/228‑1 (MQ=255)
                                                                                             ttATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTg                                                                                                                                              <  1:801576/239‑1 (MQ=255)
                                                                                              tATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTTTGTGCCGGGAGGGCTGGGTTTAGGGGAGCCCGAGGGACCGAAAAAAGCTTTTGTTTTTTGTGATTGTCAGAAGAGCATTCAGAGTTTGTTGTAAAcccctaccattta                                                                                                                                             >  1:281175/1‑228 (MQ=255)
                                                                                                ggACCGATTGGTGCCCAAAAATGTTGGTTTCTGGAAAGATGGGGGCATTGGGGTTGATCTGACGTCAGGAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGAtt                                                                                                                                           <  2:1263233/238‑1 (MQ=255)
                                                                                                      aCTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTTGTTGTGCGGGAGGGCGGGGTTTCGGGCGGGCCGCGGGACCGAAAACAACATTCGGTTTTTTTGATGGATAGAAGGGCAGTAATAGTTCGTTGTAAAcccgaccaaatttgtatttac                                                                                                                                     >  2:996113/1‑220 (MQ=255)
                                                                                                                         cGCTTCTGTAAAGTGAGTGGCATTGCGGTTGATCTGACATCAGCAGGGATGCGTGTGGAGAGGGGCTCAGTGACAACGCTGCTGATTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGt                                                                                                                 <  2:1035627/239‑1 (MQ=255)
                                                                                                                            ttCTGGAAAGGTTGTGGGTTTGCGGTTGATCGGACGTCAGGAGGGATGCGTGGGGGGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTAcg                                                                                                              <  2:599638/239‑1 (MQ=255)
                                                                                                                            ttCTGGAAAGATAGTGGCTTTGCGGTTGATCTGACGTAAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTAcg                                                                                                              <  1:796207/239‑1 (MQ=255)
                                                                                                                             gcgggAAAGGTAGTGGCAGTGCGGTTGAGCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATGCGGTAcg                                                                                                              <  1:243316/235‑1 (MQ=255)
                                                                                                                                 gAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGAt                                                                                                                                                     <  1:1122526/196‑1 (MQ=255)
                                                                                                                                      aTAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGGCGGGAGGGCTGGATTTTGGGGAGCCAACGGGACCGCAAAAAGCCTTTGGTTTGTTTTATGTTTAGAAGAGGATTCTGGGGTCGTTGTTCACCCCTCACATTGTTTTTATTTTGTTCTTTTAATATTTTGTATACGGGCTGcc                                                                                                     >  2:169712/1‑237 (MQ=255)
                                                                                                                                      aTAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGAt                                                                                                                                                                    <  1:952167/176‑1 (MQ=255)
                                                                                                                                        aGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTTCGCTTTTATGTGCCGGGAGGGCTGGGTTTAGGGCAGCCCGGGGCGCCGAAAAAAGCTTTCGTTTTTTTTTATGATCAGAAGAGGAATCCGGATTCGTTGTTCAACCATCACATTGTTTTTTTTATGTTCTTTTAAAGATTGGTACTCGGTCTGCAAcc                                                                                                  >  1:968914/1‑239 (MQ=255)
                                                                                                                                           ggCATTTCGGTTTAACTGACGTCAGCAGGGATGCGTGTGGAAGTGGGCTCAGTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCggg                                                                                               <  1:508354/239‑1 (MQ=255)
                                                                                                                                            gaaTTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCggg                                                                                               <  1:153284/236‑1 (MQ=255)
                                                                                                                                            gCATTGGGGTTTATCTGGCGGCAGCGGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCggg                                                                                               <  2:968914/238‑1 (MQ=255)
                                                                                                                                                      tGATCTGGCGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGCGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTa                                                                                     <  1:976847/238‑1 (MQ=255)
                                                                                                                                                          cTGACGCCAGCAGGGATGCGTGGGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGt                                                                                 <  1:1104410/238‑1 (MQ=255)
                                                                                                                                                                 cagcagGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGTAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGt                                                                         <  1:51742/239‑1 (MQ=255)
                                                                                                                                                                  agcggGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGGTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTa                                                                        <  1:515737/239‑1 (MQ=255)
                                                                                                                                                                  agcagGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCTGCTTTGTTGTGCCGGGAGGTCTGGGTTTTGGGGAGCCAGGGGGCACGCAAACAGACTTCGTTTTTTATGATGATCAGAAGAGCATTCAGGATTCGTTGTTCACCGATCCCATTTATTATCTTTCGGTCTTTTAAAATTTGGCAAGCGGGCGTCCACCGGGGGGTCCGGGAGAGTTTCggggga                                                                        >  1:1188431/1‑237 (MQ=255)
                                                                                                                                                                   gcagGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGGCCGCTTTGATGTGCCGGGAGGGCTGGGTTTTGGGCAGCCCGTGGGACCCGAAACAGACTTCCTTTTGTTTGATGATCAGAAGAGGATTCAGGGTTCGTTGTACAACGATCACATTTATTTTTTTATTTTTTTTTAAGATTTGGGTCGGGGTCTTCAACCCGGGGGTCCGGGAGAGGTCCggggggt                                                                       >  2:1073638/1‑236 (MQ=255)
                                                                                                                                                                       ggATGGGTGTGGAAGTGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTTAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCg                                                                    <  1:865981/238‑1 (MQ=255)
                                                                                                                                                                         attcGTGTGGAGATGGGCTCATTGACGACGCTGCTGAGTGGCGGTGTCAGCTTGGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGcc                                                                  <  2:922586/235‑1 (MQ=255)
                                                                                                                                                                            cGTGTGGAGATGGGCTCATTGACAACGCTGGTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTgg                                                               <  1:123695/238‑1 (MQ=255)
                                                                                                                                                                             gtgtGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGAtt                                                                                                                                           <  1:456791/162‑1 (MQ=255)
                                                                                                                                                                              tttgGGGATGGGCTCATTGACAACGCTGTTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGt                                                             <  2:376061/236‑1 (MQ=255)
                                                                                                                                                                                   atggtgaaaagagacgggcggtgctgAGTGGCGGTGTCAGCTTTGATGGGCCGGGAGGTTTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGt                                                        <  1:565998/217‑1 (MQ=255)
                                                                                                                                                                                       gggTTCATTGAAAACGGTGCTGAGTGGCGGTGTCAGCTTTGATGGGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTATGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCa                                                   <  2:1115786/239‑1 (MQ=255)
                                                                                                                                                                                       gggCCCATTGACAACGGGGGGGAGGGGGGGTGTCAGCTTTGAGGTGCCGGAAGGTCTGGATTTAGGGCAGGCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCa                                                   <  2:281175/239‑1 (MQ=255)
                                                                                                                                                                                        tgCTCATTTACAACGCTGGTGAGTGGCGGTGGCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCa                                                   <  1:394238/237‑1 (MQ=255)
                                                                                                                                                                                         gCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGACGTGCCGGGAGGTCTGGATTTAGGGCAACCAGTGGCAACGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCCACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCa                                                   <  1:886613/237‑1 (MQ=255)
                                                                                                                                                                                                 aCAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCa                                                   <  1:443217/229‑1 (MQ=255)
                                                                                                                                                                                                         ctgctggtggcgtggttggtttttATGTGCCGGAAGTGCTGGATTTAGGGCAGCCAGTGGGACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCCGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCAAAGTGCCATTCTTTgcg                                  <  2:523148/219‑1 (MQ=255)
                                                                                                                                                                                                               tGGCGGTTTCAGCTTTGATTTGCCGGAAGGTTTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCAAAGTGCCATTCTTTGCGCCGAat                            <  2:229544/238‑1 (MQ=255)
                                                                                                                                                                                                                ggCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCAAAGTGCCATTCTTTGCGCCGAata                           <  1:697276/238‑1 (MQ=255)
                                                                                                                                                                                                                       tCAGCTTGGATGTGCCGGGAGGTCTGGATTTAGGGCAGCAAGGGGCACCGAAAACCGCTGTCGTTTTGTATGATGATTAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCAGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCAAAGTGCCATTCTTTGCGCCGAATATGCGTCAg                   <  1:604960/239‑1 (MQ=255)
                                                                                                                                                                                                                        cAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCAAAGTGCCATTCTTTGCGCCGAATATGCGTCAg                   <  1:503724/238‑1 (MQ=255)
                                                                                                                                                                                                                        cAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCAAAGTGCCATTCTTTGCGCCGAATATGCGTCAg                   <  1:1073638/238‑1 (MQ=255)
                                                                                                                                                                                                                         aGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTCGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCAAAGTGCCATTCTTTGTGCCGAATATGCGTCAg                   >  1:784698/1‑237 (MQ=255)
                                                                                                                                                                                                                         aGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGACAGTGGCACCGAAAACAGCTTTCGTTTTGTATGGTGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGGGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGGACCGCAAGCAAAGGGCCATTCTTTGGGCCGAAAATGCGTCGGaa                 >  2:526053/1‑238 (MQ=255)
                                                                                                                                                                                                                                gATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGTTCCGGAAGAGTTCCGGGGGATTGGCCTGGGAACGGTAAGCAAAGTGCCATTCTTTCCCCCGAACTAGCGTAAGAAATTaaa           >  2:1350873/1‑239 (MQ=255)
                                                                                                                                                                                                                                 aTGTGCCGGAAGGTTTGGATTTAGGGCAGCCAGTGGCACCGAAACCAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCAAAGTGCCATTCTTTGCGCCGAATATGCGTCAGACATTTAAc          <  1:566720/239‑1 (MQ=255)
                                                                                                                                                                                                                                  tgtgCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTacac                                                                                                                                                         <  1:1163518/96‑1 (MQ=255)
                                                                                                                                                                                                                                  tgtgCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTacac                                                                                                                                                         >  2:1163518/1‑96 (MQ=255)
                                                                                                                                                                                                                                      ccGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCAAAGGGCCATTCTTTGCGCCGAATATGCGTCAGCCATTTAAGgataa     >  1:240075/1‑237 (MQ=255)
                                                                                                                                                                                                                                       ggtaaGGTTTTGATTTAGGGCAGTCGGTGGGCCCGTAAACAGTTTTCGTTTTAGATGATGATCAGGAGATCATTCAGTATTGGTTGTACACCTATCGCGTTGATTAGCTGATGTTCTGTAAAGGTTCGGTACGCGGTATGACACCGGGAGCGCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCAAAGTGCCATTCTTTGCGCCGAATATGCGTCAGACATTTAACgatga     <  2:417816/235‑1 (MQ=255)
                                                                                                                                                                                                                                          aaGGTCTGGATTTAGGGCACCCAGTGGCACCGAAAACAGTTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCAAAGTGCCATTCTTTGCGCCGAATATGCGTCAGACATTTAACGATGATTa  <  1:1052372/238‑1 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                
GCTGTTCCGTGGCTATCGGGTAGGTTCGGTTGAAACCAGCACCTTCGATACACAAAAACGCAATATCAGCTATCAACTGTTCATCAATGCACCTTATGACCGACTGGTGACCAACAATGTTCGCTTCTGGAAAGATAGTGGCATTGCGGTTGATCTGACGTCAGCAGGGATGCGTGTGGAGATGGGCTCATTGACAACGCTGCTGAGTGGCGGTGTCAGCTTTGATGTGCCGGAAGGTCTGGATTTAGGGCAGCCAGTGGCACCGAAAACAGCTTTCGTTTTGTATGATGATCAGAAGAGCATTCAGGATTCGTTGTACACCGATCACATTGATTATCTGATGTTCTTTAAAGATTCGGTACGCGGTCTGCAACCGGGAGCTCCGGTAGAGTTCCGGGGTATTCGCCTGGGTACCGTAAGCAAAGTGCCATTCTTTGCGCCGAATATGCGTCAGACATTTAACGATGATTA  >  NC_000913/1013786‑1014256

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