Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,298,945 G→A 7.3% H141H (CAC→CAT fdhF ← formate dehydrogenase‑H, selenopolypeptide subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,298,9450GA 7.3% 159.8 / 4.4 55H141H (CAC→CATfdhFformate dehydrogenase‑H, selenopolypeptide subunit
Reads supporting (aligned to +/- strand):  ref base G (25/25);  new base A (2/2);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

TGACATGGCCCATCGCATTCAACAGCGCGATGTTCGAGCCGTTTTTCAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATA  >  NC_000913/4298668‑4299232
                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                               
ggtataagggacaggaaTTAAACAGCGGGGTGTTCGAGCCGTTTTTAAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTAAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTGAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCATGTCattattg                                                                                                                                                                                                                                                                                      <  1:1320/273‑7 (MQ=255)
    aTGGCCCATCGCATTCAACAGAGCGATGTTCGAGCCGTTTTGCAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAgcg                                                                                                                                                                                                                                                                                  <  2:230517/289‑1 (MQ=255)
       gCCCATCGCATTCAACAGCGCGATGTTCGAGCCGTTTTTCAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTCCGTTTAGCGTTAATTACGTTATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTCTTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGCTTGGTGCAGACCTTAACCCGTTGGGCCGTGTCAGACACGAgcgcg                                                                                                                                                                                                                                                                                >  1:310630/1‑287 (MQ=255)
           gtCGCATTCAACAGGGCGTGGTTCGAGCCGTTTTTCACTGAAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATAGCAGACGATAATTTTCGCCCCGTTACGTTGAGCGTTAATTACGTGATTCGCCACGATTGGGTGTGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGt                                                                                                                                                                                                                                                                           <  1:299219/288‑1 (MQ=255)
           aTCGCATTCAACAGCGCGATGTTCGAGCCGTTTTTCAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGt                                                                                                                                                                                                                                                                           <  2:492408/289‑1 (MQ=255)
            tCGCATTCAACAGCGCGATGTTCGAGCCGTTTTTCAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGt                                                                                                                                                                                                                                                                           >  2:217229/1‑288 (MQ=255)
             cGCATTCAACAGCGCGATGTTCGAGCCGTTTTTCAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCa                                                                                                                                                                                                                                                                         >  1:325006/1‑289 (MQ=255)
                  tCAACAGCGCGATGTTCGAGCCGTTTTTCAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTCTCGCCCCGTTACGTTTAGCCTTAATTACGTGTGCCGCCACGCTTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTAGTATCAATTTCGTTAATAGCATTGCGCATTCCGCCATTACCGACCGATTGGTGCAGACCCGCAACCGATGGGCAGTGTCAGGCACGAGCGCAGCCGTCACCGtt                                                                                                                                                                                                                                                                    >  1:118139/1‑289 (MQ=255)
                        gcgcGATGTTCGAGCCGTTTTTCAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGt                                                                                                                                                                                                                                                              >  2:466600/1‑289 (MQ=255)
                         cgcgATGTTCGAGCCGTTTTTCAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGt                                                                                                                                                                                                                                                              <  1:48888/288‑1 (MQ=255)
                         cgcgATGTTCGAGCCGTTTTTCAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCATGTCattattgaagaaaatct                                                                                                                                                                                                                                                                            >  2:1320/1‑258 (MQ=255)
                          gcgATGTTCGAGCCGTTTTTCAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGt                                                                                                                                                                                                                                                              <  2:325006/287‑1 (MQ=255)
                          gcgATGTTCGAGCCGTTTTTCAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCATGTCattat                                                                                                                                                                                                                                                                                        <  2:262196/261‑5 (MQ=255)
                          gcgATGTTCGAGCCGTTTTTCAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCATGTCattat                                                                                                                                                                                                                                                                                        >  1:262196/1‑257 (MQ=255)
                                   gAGCCGTTTTTCAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACg                                                                                                                                                                                                                                                    >  2:450964/1‑288 (MQ=255)
                                           tttCAGTGCAATGTTCATGTCAGCAATGCGCGCCGTTTCAATTTTGCGCGGATCCCAGACGATAATTTTCGCCCCGTTACATTTAGCGTTAAGTACGTGATTGGCCACGATTGGGGGGGAATCCGCCGGATTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCTATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCa                                                                                                                                                                                                                                           <  2:318422/289‑1 (MQ=255)
                                                   cAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAACGCAGCAGTCAACGTTATTGGTACCAATAACCGCCCGACCAAAttttg                                                                                                                                                                                                                                    >  1:494057/1‑287 (MQ=255)
                                                   cAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTCCGCGGATCGTAGACGATACTCTTCGCCCCGTTACGTTTAGCGCTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAATTCGAACACTAAATCGGTATTATCAATTACGTTAATAGCATTGCTCATTGCGCCATTACCGACCGACTCGTGCAGACCTGCAACCGATGGCCCGTGTCAGACACTAGCGCCACCGTCCACGATATTGGTACCCATAACGGCGCGCACaatccttt                                                                                                                                                                                                                                    >  1:137540/1‑282 (MQ=255)
                                                     aggtgCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCa                                                                                                                                                                                                                                 <  2:193628/287‑1 (MQ=255)
                                                     aTGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCa                                                                                                                                                                                                                                 <  2:310630/289‑1 (MQ=255)
                                                      tgtgCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCAt                                                                                                                                                                                                                                >  2:35219/1‑289 (MQ=255)
                                                           aTGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACAt                                                                                                                                                                                                                           >  2:367791/1‑289 (MQ=255)
                                                               cAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGt                                                                                                                                                                                                                        <  1:466600/288‑1 (MQ=255)
                                                                     tGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTc                                                                                                                                                                                                                 >  2:89278/1‑289 (MQ=255)
                                                                       cgcgcgGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTc                                                                                                                                                                                                                                                                          <  2:125564/230‑1 (MQ=255)
                                                                       cgcgcgGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTc                                                                                                                                                                                                                                                                          >  1:125564/1‑230 (MQ=255)
                                                                                 cAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTAcc                                                                                                                                                                                                     <  1:217229/289‑1 (MQ=255)
                                                                                        gcgcGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCgag                                                                                                                                                                                              >  2:64916/1‑289 (MQ=255)
                                                                                                cGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGgt                                                                                                                                                                                      <  1:89278/289‑1 (MQ=255)
                                                                                                                                 cGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGt                                                                                                                                                                                                               <  1:328903/231‑1 (MQ=255)
                                                                                                                                 cGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGt                                                                                                                                                                                                               >  2:328903/1‑231 (MQ=255)
                                                                                                                                 cGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCttt                                                                                                                                                     >  2:189818/1‑289 (MQ=255)
                                                                                                                                           gTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATTGCCTCCGGACCGTACTTCTCTTTGATGGCGCt                                                                                                                                            >  2:450497/1‑288 (MQ=255)
                                                                                                                                           gTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGa                                                                                                                                                   <  1:55643/281‑1 (MQ=255)
                                                                                                                                           gTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGa                                                                                                                                                   >  2:55643/1‑281 (MQ=255)
                                                                                                                                               ttCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGgc                                                                                                                                       >  2:175806/1‑289 (MQ=255)
                                                                                                                                                                      ggccgGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGcc                                                                                                                <  2:31440/288‑1 (MQ=255)
                                                                                                                                                                      cgccgGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCCCTCGGCAACGTAATTCAGttc                                                                                                                 >  1:5895/1‑286 (MQ=255)
                                                                                                                                                                       gccgGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCt                                                                                                               <  1:189818/289‑1 (MQ=255)
                                                                                                                                                                                                         tattatCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCAc                                                                              >  2:28819/1‑288 (MQ=255)
                                                                                                                                                                                                         tattatCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACg                                                                             <  1:450964/289‑1 (MQ=255)
                                                                                                                                                                                                                   ttCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACgg                                                                   >  2:19922/1‑289 (MQ=255)
                                                                                                                                                                                                                          aTAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATgg                                                            <  1:64916/289‑1 (MQ=255)
                                                                                                                                                                                                                            aGCATTGCTCGTGGGGCCATTACCGACCGATTGGTGCAGACCTGCAACCGGTGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATgggg                                                          <  1:28819/289‑1 (MQ=255)
                                                                                                                                                                                                                                ttGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGCTATTGGTACCAATAACGGCGAGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCCCCGCCACGTACTTCAGTGGCTACTCACAGGAAACAGGCTCGAGTTTGCCGCACCGCTGGCTACAGATCAGGGGGGtgt                                                      >  1:317405/1‑287 (MQ=255)
                                                                                                                                                                                                                                                tACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAg                                      <  1:175806/289‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccgaccgaTTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGa                                    >  2:131601/1‑289 (MQ=255)
                                                                                                                                                                                                                                                  ccgaccgaTTGGTGCAGACCTGCAACCGAAGGGCCTGGTCAGACACGAGCGCGGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGa                                    <  1:35219/289‑1 (MQ=255)
                                                                                                                                                                                                                                                    gaccgaTTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGAGCATGGGGGTTTTCAGGCGCGGGGTCAGGaca                                  >  2:389521/1‑287 (MQ=255)
                                                                                                                                                                                                                                                         aTTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGt                             <  1:19922/289‑1 (MQ=255)
                                                                                                                                                                                                                                                         aTTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCCCAGGCGCTCGGCAACGTAATTCAGTGCCCCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGTTGGCGACGGATCATGGGGGTTTTCATGCGCGGGGTCAGGATCTGGGt                             >  2:430795/1‑289 (MQ=255)
                                                                                                                                                                                                                                                                         aaCCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTc              <  2:78878/288‑1 (MQ=255)
                                                                                                                                                                                                                                                                            caaTGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGGGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGGGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAg          <  1:474722/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                             gATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTGTTGGGACCTATTACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCCTCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGGTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGc         <  2:118139/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    cGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCataata  <  1:450497/289‑1 (MQ=255)
                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                               
TGACATGGCCCATCGCATTCAACAGCGCGATGTTCGAGCCGTTTTTCAGTGCAATGTGCATGTCAGCAATGCGCGCGGTTTCAATTTTGCGCGGATCGCAGACGATAATTTTCGCCCCGTTACGTTTAGCGTTAATTACGTGATTCGCCACGATTGGGTGGGAATCCGCCGGGTTGTACCCGAAAACGAACACTAAATCGGTATTATCAATTTCGTTAATAGCATTGCTCATTGCGCCATTACCGACCGATTGGTGCAGACCTGCAACCGATGGGCCGTGTCAGACACGAGCGCAGCAGTCAACGTTATTGGTACCAATAACGGCGCGCGCAAATTTTTGCATTACATAGTTGGTTTCGTTACCCGTACCACGCGAGGAGCCGGTCGTCTGGATGGCATCCGGACCGTACTTCTCTTTGATGGCGCTCAGGCGCTCGGCAACGTAATTCAGTGCCTCATCCCAGGAAACAGGTTCGAGTTTGCCGCCACGCTGGCGACGGATCATGGGGGTTTTCAGGCGCGGGGTCAGGATCTGGGTATCGTTAATGAAGTCCCAGCCATAATA  >  NC_000913/4298668‑4299232

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: