Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,368,894 C→T 18.3% S77L (TCG→TTG)  fxsA → suppressor of F exclusion of phage T7

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,368,8940CT18.3% 122.1 / 4.6 71S77L (TCG→TTG) fxsAsuppressor of F exclusion of phage T7
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base T (12/2);  total (39/32)
Fisher's exact test for biased strand distribution p-value = 1.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCG  >  NC_000913/4368663‑4369131
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         
cTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGAt                                                                                                                                                                                                                                                    <  2:889221/238‑1 (MQ=255)
    gcgcTGGTTACCTTTTATTGCCATTTTCCTTTATTGCTTTAATGAGAGTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGGCCCTCGTGCTGGTTATATTCACGGCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTg                                                                                                                                                                                                                                              <  1:58405/239‑1 (MQ=255)
             aCCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGGGATGGTTAAAAGTGTTTCGCTGGTCATTGCTGGgtt                                                                                                                                                                                                                                     >  2:813332/1‑237 (MQ=255)
              ccTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTtgct                                                                                                                                                                                                                                  <  2:841909/239‑1 (MQ=255)
                            tttCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTTGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCgg                                                                                                                                                                                                                 <  2:845669/239‑1 (MQ=255)
                             ttCCTTTATTTATATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCggg                                                                                                                                                                                                                <  2:752239/239‑1 (MQ=255)
                                  ttATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGGTTAAAAGGGTTTCGCTGGTCATTGCTGGTTTGCGGCCTTTTTTTACGGGCttttt                                                                                                                                                                                                          >  2:661471/1‑239 (MQ=255)
                                         tataTTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGGTGATTAAAAGGGTTTCGCTGATCATTGCTGGTTTTCTGCTTTTATTACCGGGCTTTTTTACCGa                                                                                                                                                                                                    >  2:652871/1‑239 (MQ=255)
                                          atatTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGTCCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGGAGTGGGTGAAAAACACACGGCTGAGATGGTTAACAGGGATTCCGTGATTATTGATGGTTTTCTGCTTGTATTATCGCGCTTTTTTAcctgg                                                                                                                                                                                                   >  2:344850/1‑236 (MQ=255)
                                          atatTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTGGTGCTGATGCAGCAAAAAATGGCGGCGGGTGGAAAACCAGCGGCGGAGATGATTTAAAGTGTTTTGCTTCTCCTTGCTGGTTTGCTCCCTTTATTACCGGGCTTTTTTTCCg                                                                                                                                                                                                     >  1:262126/1‑237 (MQ=255)
                                          atatTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGGTTTGCTGATCACTGCTGGTTTGCTGGCTTTTTTACCCGGCGTTTTTTCCg                                                                                                                                                                                                     >  2:874328/1‑237 (MQ=255)
                                           tatTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCACCGCCGGAGATGGTTAAAAGTGGTTTGCCGATCCTTGCTGGGTTGTTGCCTTTATTAACGGGCTTTTTTACCGCCt                                                                                                                                                                                                 >  2:852915/1‑239 (MQ=255)
                                                agaTTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGGGTTTGTTTCAAGTGTTTTGTTGGTATTTGCTGGTTTGCTCTCTTTTTTTACCGGTGttttttcctactttcc                                                                                                                                                                                               >  1:413062/1‑228 (MQ=255)
                                                agaTTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCCTTTATTAACGGGCTTTTTTACCGACTTCCTc                                                                                                                                                                                            >  2:393352/1‑239 (MQ=255)
                                                 gaTTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGAACCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGGAATGGTTTAAAGGGGTTCGCTGGTCATTGCGGGTTTGCTGCTTTTTTTTCCGGGCTTTTTTTCCCCCTTccctg                                                                                                                                                                                           >  2:257103/1‑236 (MQ=255)
                                                       aaTTTTTTTTTAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGtcttc                                                                                                                                                                                    <  1:202147/239‑1 (MQ=255)
                                                            ttATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGGGTTTCGCTGATCATTGCTTGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACtt                                                                                                                                                                              >  2:811683/1‑239 (MQ=255)
                                                             tATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGCGGAAACCCAGCGGGGGAGGGGAGTTAAAAGTGTTTTCCTTCACTATTCGGTTTTTGTGGTTTTTTTAACGGGCTTTTTTTCCGatcctcccgggtttctcttt                                                                                                                                                                                 >  1:825308/1‑220 (MQ=255)
                                                                       gCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGGGATGGTTTAAAAGGTTTCTCCGATTATTTGCTGGTTTGCTGTTTTTTTTACGGGCTTTTTTTTCCCCCTCCTCGGGTTTCTACCTTTATTtccgcc                                                                                                                                                                   >  1:326844/1‑239 (MQ=255)
                                                                            tttATTGGGGGTTTTGCTGACCCCCGGTCTGGTTATATTCACGTCAGTTTGCGGTATGTCCCTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGc                                                                                                                                                             <  1:696229/237‑1 (MQ=255)
                                                                                  gggggTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTTGAAACCCAGCGGCGGAGATGGTTAAAAAGGTTTCGATGGTCATTGCTGGTTTGCTGCTTTTATTTCCGGGCTTTTTTTCCGACTTCCCCGGGTTTCTACTTTTTTTTCCCCCGGGTGCaaaac                                                                                                                                                        >  2:859486/1‑238 (MQ=255)
                                                                                   ggggTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAACGGCGGGGATGATTAAAAGGGGTTCGCTGATCCTTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGGCTTCCTCGGTCTTCTACTTTTATTTTCCCCGGGTCaaaaaaa                                                                                                                                                      >  2:128343/1‑236 (MQ=255)
                                                                                     ggTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGGGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGGCTTCTACTTTTTTTTCCCCCGGGGCCAAAGCa                                                                                                                                                      >  2:185141/1‑237 (MQ=255)
                                                                                       tATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTTCTGCTTTTATTTACGGGCTTTTTTACCGACTTCCCTGGGCTTCTACTTTTATTGCCGCCGGTGCCAAAACATCtt                                                                                                                                                  >  2:262414/1‑238 (MQ=255)
                                                                                       tATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGGAAACCCAGCGGCGGGGAAGATTAAAAAGTTTTTGCCGGCAATTGCTGGTTTGCTGGTTTTATTTCCGGGCGTTTTTCCCCACTTTCTCGGGCTTTTACTTTTTTTGCCGCCCGTGCaaaaaaaactg                                                                                                                                                  >  1:583286/1‑232 (MQ=255)
                                                                                        tttggcgccccttgtttttgTTATTTTCACGTCAGTTATCGGTATGTCACTGGTACGGAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGa                                                                                                                                                 <  1:185141/220‑1 (MQ=255)
                                                                                        aTTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTGAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAAACCAGCCGCGGGGGAGATTAAAAATGGTTTGCTGGTCATTGTTTGGTTGTTTGCTTTTTTTACCGGGCTTTTTTCCGACTTTCCCGGTCTTTTACCTTTTTTTCCCCCCGGTGCaaaaaaactt                                                                                                                                                   >  1:345390/1‑232 (MQ=255)
                                                                                          tGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCCGCGGGGGGGGTGGATTAAAAGGGTTCGCTGGACCATGCGAGGTTTTTGCCTTTTTTAACGGGCTTTTTTTACGCCTTTCCCCGGTTTTTTTCTTTTTTTCCCCCCGGGGAAAAAACCTCTCAc                                                                                                                                                >  1:665436/1‑239 (MQ=255)
                                                                                          tGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGGTCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCCACTTCCTCGGTCTTCTACTTTTATTTCCGCCCGTGCAAAAGCATCttgca                                                                                                                                               >  2:364658/1‑235 (MQ=255)
                                                                                          tGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGa                                                                                                                                                 <  2:700114/237‑1 (MQ=255)
                                                                                            cttccccTCTTGCTGTTTATTTTTACTTCAGTTATGGGGATGTCACGGGTACGAAACCAGGGTTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGt                                                                                                                                             <  1:849258/235‑1 (MQ=255)
                                                                                              gACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCCGCGGGTGAAAACCCAACGGCGGGAAATCATAAAAGGGGTTCGCTGCACATATCCGGGTTTCTACGTTTATTATCAGGGTATTTTTACGTCCTTACTCGTCTTCTACATTTTTTTGCGCCAGTGCTAAAGCATcacaacgaga                                                                                                                                           >  1:266155/1‑230 (MQ=255)
                                                                                               aCCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGCGGCGGAGATGATTAAAAAGGGTTTGCTGATCATTTCTGGGTTGCTGCCTTTTTTACCGGGCTTTTTTTCCGACTTCCCCGGTGTTTTACTTTTTTTTCCGCCCGGGCCAAAGCAACTTACAGTGaa                                                                                                                                         >  2:215177/1‑239 (MQ=255)
                                                                                                    cGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGCGGCGGGGTTGATTTAAAAAGTTTTGCTGAGATATTCTGGGTTGGTGCTTTTTTTATCCGGGTTTTTTTCCACCATTCTCGGGTTTTCTCTTTTTTTTGCCCCCGGGGAAAAAACATCttaaaggaagtgtg                                                                                                                                      >  1:489241/1‑226 (MQ=255)
                                                                                                    cGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGCGGAAACCCACCGGCGGGGATGAATAAAAAGGTTTTGCCTGATATTGCCGGGTTTTTGCTTTTTTTAACAGGGTTTTTTTCCCGATTCCCCGGCCTTCTACTTTTTTTGCCCCCCGGGGCCAAGCATCtctaccttgaatt                                                                                                                                       >  1:700114/1‑225 (MQ=255)
                                                                                                       tcTGGTTTTTTTTCCGTCAGTTTTCGGTATGTCCCTGGTACGGAACCAGGGCTTTAAGAATTTCGGGCTGATGCAGCAAAAAATGGGGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGc                                                                                                                                 <  1:215177/238‑1 (MQ=255)
                                                                                                       gttgtttATATTCATGTTATTTTTCGGTTTGTCACTGGTACGGAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGc                                                                                                                                 <  1:257103/234‑1 (MQ=255)
                                                                                                       gcgggTTATATTCACGTCAGTTAACGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGAt                                                                                                                                                                                                                                                    <  2:334486/132‑1 (MQ=255)
                                                                                                             tataTTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCAt                                                                                                                            <  2:489241/238‑1 (MQ=255)
                                                                                                                     cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGGGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGGCTTCCTCGGGCTTCTACTTTTATTGCCGCCGGTGCCAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCttt                                                                                                                   >  2:252272/1‑239 (MQ=255)
                                                                                                                     cGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGCGGCGGAGATGATTTAAAGTGTTTCGCTGGTCCTTGCTGGTTTTCGGCTTTTTTTACCGGGGTTTTTTTCCGCCTTCCTCGGTCTTTTTCTTTTTTTTCCCCCCGGGCAAAAACATCTTACAGGGAAGTTGGTTCCGCATTTTCGGttt                                                                                                                   >  1:708149/1‑239 (MQ=255)
                                                                                                                      gTCAGTTATCGGTATGTTAGTGGTTCGTAACCAGGGCTTTAAGACTTTCGTGCGTATTCAGAAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCtttt                                                                                                                  <  1:96926/239‑1 (MQ=255)
                                                                                                                         aGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGGTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTTCCGCCCGTGCCAAAGCATCTGCACGTGAAGTTGTTTTCGCAATTGCGGTTTtct                                                                                                                >  2:447297/1‑238 (MQ=255)
                                                                                                                           ttATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAACCCAGCGGCGGAGAGGATTAAAAGTGTTTCGCTGATCATTGCTGGGTTGCTGCTTTTTTTACCGGGCTTTTTTACCGACTTCCTCGGTCTTTTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACCGGGAAGTTGTTGCCCCATTTGCGCTTTtttccc                                                                                                             >  2:328609/1‑237 (MQ=255)
                                                                                                                            tATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCa                                                                                                            <  2:708149/239‑1 (MQ=255)
                                                                                                                             aTCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGGGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTTTTTCCCCCGGTGCAAAAGCACCTTACCGTGGAGTTGATTCCGCCTTTTGGCtttttctgcc                                                                                                            >  2:96926/1‑232 (MQ=255)
                                                                                                                                     gTCACTGGAACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGGGAAAAACCAGCGGCGGAGATGCTTTAAAAAGTTTTGCCGAACCTTGCTCGGTTGCCTGCTTTCTTACCGGGCTTTTTTTCCCCGTTTCCCGGCTTTTTACTTTTTTTGCACCCCGTGCAAAAACAATTGCGAGTGTAGGTTCATGCCGAATTGGCGCTTTTTCCGCAGGCCgggg                                                                                                    >  1:505322/1‑238 (MQ=255)
                                                                                                                                            gTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGGGATGGTTAAAAGGGTGTCGCTGGTCCTTTGTGGTTTGATGCTTTTTTTTCCTGGCCTTTTTTACGCGTTTCTCGGTCTTCTACATTTTTTTCCGCCGATGCAAAACCATCTGACGGTGACATTCTTACCGCATATGCGCGTTTTTCGCacaccgcggggagggga                                                                                            >  1:370880/1‑223 (MQ=255)
                                                                                                                                             tgcGTAACCAGGGCTTTACACTTTTGGTGGTGATGCACCAAAAAATGGGGGCGGGTTAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCCTGTATGACGGTGTTATTCCCGGGCTTTTTTACCGACTTCCTCGGTCTGCTACTGTTATTGCCGCCGGTGCAAAATCATCTGACAGTGAAGTTGATGCCGCATTTGCGCGTTTCTCGCATGCCTGGCGGCGGtttt                                                                                           <  1:500616/237‑1 (MQ=255)
                                                                                                                                               cGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGGGTTTCGCTGGTCATTGCTGGTTTGCCGCTTTTATTACCGGGCTTTTTTACCGCCTTCCTCGGTCTTCTACTTTTATTTCCGCCGGTGCAAAAGCATCTGACAGGGAAGTTGCTTCCGCATTTTCGCTTTTTTCGCATGCCTGGCGGCGGGTTTat                                                                                         >  2:782451/1‑238 (MQ=255)
                                                                                                                                                   aCCAGGGCTTTAAGAATTTGGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCAGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGcc                                                                                     <  2:262126/239‑1 (MQ=255)
                                                                                                                                                       gggCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGGGCTGGTTAAAAAGGGTTTGCTGATCATTGCTGGTTTTCTGCTTTTTTTTCCCGGCCTTTTTTACGCCTTCCTCGGTCTTCTACTTTTTTTTCCGCCCGTTGAAAAGCATCTGACAGGGGAGTTGTTGCCGCATTTGCGCTTTTCTCGCATGCCCGGGGGCGGGTTTTGAGCCgcg                                                                                  >  1:176333/1‑236 (MQ=255)
                                                                                                                                                             tAAGAATTTCGTGGTGATGCAGCAAAAAATGGCTGCGGGTGAAACCCCATCGGCGGAGATTATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTGCCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTg                                                                           <  2:825308/239‑1 (MQ=255)
                                                                                                                                                                     tCGTGCTTGTGCAGCAAAAAAGGGCGGCGGGTTAAACCCCGGGGGCGGAGATGATTAAAAGTGTTTCGCTGATCATGTCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAAt                                                                   <  1:852915/239‑1 (MQ=255)
                                                                                                                                                                     tCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAAt                                                                   <  2:527117/239‑1 (MQ=255)
                                                                                                                                                                     tCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAAt                                                                   <  2:277314/239‑1 (MQ=255)
                                                                                                                                                                                       aaaTGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTAc                                                 <  2:326844/239‑1 (MQ=255)
                                                                                                                                                                                             ggggtggTGAAACCCCAGGGGCGGAGGTGATTAAAAGTGTTTCGCTGGTCATTGCTGGTTTGGTGCTTTTATTACCGGGTTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCTCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                                            <  2:266155/233‑1 (MQ=255)
                                                                                                                                                                                             cggcggGTGAAAACCCAGCGGGGGAGATGATTAAAAGTTTTTCCGTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCCCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCg                                            <  1:811683/238‑1 (MQ=255)
                                                                                                                                                                                               ggggggGAAAACCCAGCGGGGGTGATTATTAAAGGGGTTTCCGGGATCTTTGCTGTTTTTCTGCTTTTATTACCGGGTTTTTTTCCCGAGTTCCTCGTTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGaa                                          <  1:447297/236‑1 (MQ=255)
                                                                                                                                                                                                      taacacccggtggttggtttgatTAAAAGTTTTTTGCTGATCTTTGCTGTGTTGCCTTTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGgatg                                    <  1:93639/219‑1 (MQ=255)
                                                                                                                                                                                                      aaaaCCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGGTGGTACCACCGAAAggggtat                                  >  2:45348/1‑234 (MQ=255)
                                                                                                                                                                                                      aaaaCCCAGCGAGGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGgatga                                   <  2:505322/238‑1 (MQ=255)
                                                                                                                                                                                                        gaaacccccggggtgggtttttttAAATGTTTTTTGCTTTTCATTGCTGGTTTGCTGCTTTTATTTCCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATGGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGgatg                                    <  1:262414/216‑1 (MQ=255)
                                                                                                                                                                                                            cAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGGCCGGTGGCGGGAATACTTTTGATGGTGGGGACCAGCGAAAggaggagtggcgc                            >  2:811843/1‑232 (MQ=255)
                                                                                                                                                                                                              ggggcggAAATGGTAAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGTTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGa                          <  2:665436/237‑1 (MQ=255)
                                                                                                                                                                                                                       ggtgatTAAAAGTTTTACGCAGATATTTGCTGGTTTGCTGCTTTTATTGCCGGGCTTTTTTACCGACTGCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGGGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCtt                   <  2:312176/236‑1 (MQ=255)
                                                                                                                                                                                                                        atgatTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTg                  <  2:851937/238‑1 (MQ=255)
                                                                                                                                                                                                                               tAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGCAAGGGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGGTTTAGCGCCGGGGCCGGGGGGGGCAATACTTTTTATTGCGGGGAACAGCGACAGGGAGATGGGCGAGAGCGCCtcggaccaa            >  1:312176/1‑231 (MQ=255)
                                                                                                                                                                                                                                       tttCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGATCCCGACCGCCTTGATCTTAAAGACGa     >  2:256678/1‑238 (MQ=255)
                                                                                                                                                                                                                                       tttCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGGAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACgc     >  2:555116/1‑237 (MQ=255)
                                                                                                                                                                                                                                          cGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATTGTGAGTACCCGCGGAAGGATGATGAGCGCGACCGCCTTGATCCTAAAGACGACCg  >  2:833240/1‑238 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                         
CTTTGCGCTGGTTACCTTTTATTGCCATTTTCCTTTATGTCTATATTGAGATTTCAATCTTTATTCAGGTTGCCCATGTATTGGGGGTATTGCTGACCCTCGTGCTGGTTATATTCACGTCAGTTATCGGTATGTCACTGGTACGTAACCAGGGCTTTAAGAATTTCGTGCTGATGCAGCAAAAAATGGCGGCGGGTGAAAACCCAGCGGCGGAGATGATTAAAAGTGTTTCGCTGATCATTGCTGGTTTGCTGCTTTTATTACCGGGCTTTTTTACCGACTTCCTCGGTCTTCTACTTTTATTGCCGCCGGTGCAAAAGCATCTGACAGTGAAGTTGATGCCGCATTTGCGCTTTTCTCGCATGCCTGGCGGCGGTTTTAGCGCCGGGACCGGTGGCGGTAATACTTTTGATGGTGAGTACCAGCGAAAGGATGATGAGCGCGACCGCCTTGATCATAAAGACGATCG  >  NC_000913/4368663‑4369131

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