Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,028,197 Δ1 bp intergenic (‑119/+18) PA14_RS22975 ← / ← PA14_RS22980 Nramp family divalent metal transporter/ribbon‑helix‑helix domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084635,028,1950G.98.8% 2061.1 / 15.8 426intergenic (‑117/+20)PA14_RS22975/PA14_RS22980Nramp family divalent metal transporter/ribbon‑helix‑helix domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (5/0);  new base . (207/214);  total (212/214)
Fisher's exact test for biased strand distribution p-value = 2.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 463 total aligned reads displayed.

TGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGGTGGGTTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTCC  >  NC_008463/5028045‑5028338
                                                                                                                                                      |                                                                                                                                               
tGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCg                                                                                                                                                 >  2:2754937/1‑150 (MQ=255)
tGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCg                                                                                                                                                 >  2:9150487/1‑150 (MQ=255)
   gcAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCggtc                                                                                                                                              >  1:1884012/1‑148 (MQ=255)
    cAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCggtcg                                                                                                                                             >  2:7622290/1‑147 (MQ=255)
      gcggcgAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCggtcgcg                                                                                                                                           <  2:8073137/150‑6 (MQ=255)
        ggcgAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCggtcgcgtc                                                                                                                                         <  1:2120111/150‑8 (MQ=255)
                gTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTtccctcc                                                                                                                                >  2:3222788/1‑150 (MQ=255)
                gTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTtccctcc                                                                                                                                <  1:4718193/150‑1 (MQ=255)
                 tCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTtccctccc                                                                                                                               >  2:175437/1‑150 (MQ=255)
                    gccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGACCGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTccctcccctc                                                                                                                            <  2:5563527/150‑1 (MQ=255)
                    gccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAATCGTTTCACCCCGGTCGCGTCTTccctcccctc                                                                                                                            >  2:6858315/1‑150 (MQ=255)
                    gccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTccctcccctc                                                                                                                            >  1:3599918/1‑150 (MQ=255)
                     ccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCa                                                                                                                           >  1:6155/1‑150 (MQ=255)
                       ggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCCGGGAGTATTTCACGCGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTGCCCTCCCCTCAgg                                                                                                                         >  2:2133524/1‑150 (MQ=255)
                       ggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAgg                                                                                                                         <  1:9552807/150‑1 (MQ=255)
                         cACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGcc                                                                                                                       <  2:4898155/150‑1 (MQ=255)
                         cACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGcc                                                                                                                       <  2:6720237/150‑1 (MQ=255)
                         cACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGGGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGcc                                                                                                                       >  2:7915153/1‑150 (MQ=255)
                          aCCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCt                                                                                                                      >  1:4927883/1‑150 (MQ=255)
                            cGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTgg                                                                                                                    <  2:11132576/150‑1 (MQ=255)
                                aCGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgc                                                                                                                >  2:7473966/1‑150 (MQ=255)
                                  gCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcgc                                                                                                              >  1:4049443/1‑150 (MQ=255)
                                  gCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcgc                                                                                                              >  1:7824512/1‑150 (MQ=255)
                                  gCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcgc                                                                                                              <  1:547416/150‑1 (MQ=255)
                                  gCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcgc                                                                                                              >  2:10799520/1‑150 (MQ=255)
                                   cTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCt                                                                                                             <  2:5173347/150‑1 (MQ=255)
                                     gtgGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTgg                                                                                                           <  1:8690621/150‑1 (MQ=255)
                                        gCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCg                                                                                                        >  1:8406693/1‑150 (MQ=255)
                                          tGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGAc                                                                                                      <  2:4090129/150‑1 (MQ=255)
                                           gACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGAcc                                                                                                     >  2:3105686/1‑150 (MQ=255)
                                           gACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGAcc                                                                                                     <  1:9402370/150‑1 (MQ=255)
                                           gACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGAcc                                                                                                     >  2:1963350/1‑150 (MQ=255)
                                              cgccgAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAgg                                                                                                  >  1:3060681/1‑150 (MQ=255)
                                               gccgAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGc                                                                                                 >  1:8107144/1‑150 (MQ=255)
                                                ccgAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCt                                                                                                >  1:3389288/1‑150 (MQ=255)
                                                  gAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCa                                                                                              <  2:5731333/150‑1 (MQ=255)
                                                     ggACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAAGGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCt                                                                                           <  1:2488387/150‑1 (MQ=255)
                                                     ggACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCt                                                                                           <  2:4456831/150‑1 (MQ=255)
                                                           aaCGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGa                                                                                     >  1:6435413/1‑150 (MQ=255)
                                                            aCGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGAc                                                                                    <  1:5617017/150‑1 (MQ=255)
                                                            aCGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGAc                                                                                    <  1:5163878/150‑1 (MQ=255)
                                                               gcgcGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCtt                                                                                 >  2:3351699/1‑150 (MQ=255)
                                                                 gcgGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTtct                                                                               <  1:6861201/150‑1 (MQ=255)
                                                                   ggTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcg                                                                             <  2:3843576/150‑1 (MQ=255)
                                                                      tggccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcgcg                                                                           >  1:5926532/4‑150 (MQ=255)
                                                                      tggccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcgcg                                                                           >  2:1174004/4‑150 (MQ=255)
                                                                        gccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGcgcc                                                                         <  1:10504966/149‑1 (MQ=255)
                                                                         ccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGcgccg                                                                        >  2:6950259/1‑150 (MQ=255)
                                                                           gtggtgGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCGCGGCGACCTTCTCGCGCCGCt                                                                      >  1:1294908/1‑150 (MQ=255)
                                                                            tggtggGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTg                                                                     >  2:10653927/1‑150 (MQ=255)
                                                                               tggGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTgcgc                                                                  <  1:2900488/150‑1 (MQ=255)
                                                                                gggCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGcc                                                                 <  2:4785618/150‑1 (MQ=255)
                                                                                      ggAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAg                                                           >  2:9818498/1‑150 (MQ=255)
                                                                                      ggAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAg                                                           >  2:5092089/1‑150 (MQ=255)
                                                                                        aGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCg                                                         >  1:3836964/1‑150 (MQ=255)
                                                                                         gTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCgg                                                        >  1:11065416/1‑150 (MQ=255)
                                                                                           aTTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGt                                                      >  1:5533252/1‑150 (MQ=255)
                                                                                            tttCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTg                                                     <  2:10373590/150‑1 (MQ=255)
                                                                                              tCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCa                                                   <  2:297342/150‑1 (MQ=255)
                                                                                               cACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCATCTCGCGCCGCTGCGCCAGGTAGCGGGTGCaa                                                  >  1:10416526/1‑150 (MQ=255)
                                                                                                aCGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAAc                                                 >  2:3030631/1‑150 (MQ=255)
                                                                                                aCGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAAc                                                 <  2:2353106/150‑1 (MQ=255)
                                                                                                 cGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACt                                                <  1:9220741/150‑1 (MQ=255)
                                                                                                 cGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACt                                                >  1:3174196/1‑150 (MQ=255)
                                                                                                 cGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACt                                                <  1:9136916/150‑1 (MQ=255)
                                                                                                  gAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTg                                               <  1:7482212/150‑1 (MQ=255)
                                                                                                   aGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTGCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGa                                              >  2:6271669/1‑150 (MQ=255)
                                                                                                   aGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCGGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGa                                              >  1:1418625/1‑150 (MQ=255)
                                                                                                     gCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGaca                                            <  1:8561638/150‑1 (MQ=255)
                                                                                                       gACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACg                                          >  2:4581695/1‑150 (MQ=255)
                                                                                                        aCGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGc                                         >  1:10787790/1‑150 (MQ=255)
                                                                                                        aCGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGc                                         >  2:3605839/1‑150 (MQ=255)
                                                                                                        aCGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGGCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGc                                         >  2:5919475/1‑150 (MQ=255)
                                                                                                          gCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCag                                       >  2:7539077/1‑150 (MQ=255)
                                                                                                          gCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCag                                       <  1:8896521/150‑1 (MQ=255)
                                                                                                           ccGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCagg                                      >  2:8691786/1‑150 (MQ=255)
                                                                                                             ggATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggag                                    >  2:10076469/1‑149 (MQ=255)
                                                                                                             ggATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaa                                    <  1:9596947/150‑1 (MQ=255)
                                                                                                              gATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaag                                   <  1:10008568/150‑1 (MQ=255)
                                                                                                                 cGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAaggag                                <  1:2308034/150‑1 (MQ=255)
                                                                                                                  ggggAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAaggagg                               <  2:7030076/150‑1 (MQ=255)
                                                                                                                   gggAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGc                              <  1:7671712/150‑1 (MQ=255)
                                                                                                                   gggAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGc                              <  1:6625584/150‑1 (MQ=255)
                                                                                                                    ggaggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCg                             <  1:10147421/150‑1 (MQ=255)
                                                                                                                    ggaggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCg                             >  1:7350477/1‑150 (MQ=255)
                                                                                                                     gaggaAGTGCGACCAGCCTAAACGTTTTACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACGGACACGCAGGAAGGAGGCGa                            >  1:2092657/1‑150 (MQ=255)
                                                                                                                     gaggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGa                            >  1:9757056/1‑150 (MQ=255)
                                                                                                                     gaggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGa                            >  1:6357312/1‑150 (MQ=255)
                                                                                                                      aggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGaa                           <  2:2547080/150‑1 (MQ=255)
                                                                                                                       ggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAg                          <  2:1448074/150‑1 (MQ=255)
                                                                                                                        gaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGt                         <  1:1673780/150‑1 (MQ=255)
                                                                                                                          aGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTg                       <  2:2234980/150‑1 (MQ=255)
                                                                                                                          aGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTg                       >  2:1091177/1‑150 (MQ=255)
                                                                                                                             gCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGc                    >  1:9996607/1‑150 (MQ=255)
                                                                                                                                 ccAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCt                >  1:2603609/1‑150 (MQ=255)
                                                                                                                                  cAGCCTACACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTc               <  2:8603615/150‑1 (MQ=255)
                                                                                                                                     ccTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgcc            <  2:9577609/150‑1 (MQ=255)
                                                                                                                                       tAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgccgc          <  1:425624/150‑1 (MQ=255)
                                                                                                                                             tttCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGt    >  2:5350339/1‑150 (MQ=255)
                                                                                                                                               tCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCGGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTcc  >  2:8496792/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                               
TGCGCAGCGGCGAGGCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGGTGGGTTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTCC  >  NC_008463/5028045‑5028338

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