Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 4,997,786 T→C K173E (AAG→GAG)  PA14_RS22840 ← DUF853 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084634,997,7860TC99.8% 1824.0 / NA 469K173E (AAG→GAG) PA14_RS22840DUF853 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (213/255);  minor base A (0/1);  total (213/256)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 469 total aligned reads displayed.

GCCCCCGCCCGTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTTGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTGGC  >  NC_008463/4997639‑4997933
                                                                                                                                                   |                                                                                                                                                   
  ccccGCCCGTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGa                                                                                                                                                <  2:10513108/151‑1 (MQ=255)
   cccgcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAg                                                                                                                                               >  1:6101794/1‑151 (MQ=255)
   cccgcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAg                                                                                                                                               >  1:496253/1‑151 (MQ=255)
    ccgcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCTGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAgg                                                                                                                                              <  1:7369631/151‑1 (MQ=255)
      gcccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGAt                                                                                                                                            <  1:3059468/151‑1 (MQ=255)
        ccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGAtct                                                                                                                                          <  1:10162421/151‑1 (MQ=255)
        ccgTCCGGCTCGGGGCGGAGGATGTCTTCCAGGTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGCGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGAtct                                                                                                                                          <  1:4126574/151‑1 (MQ=255)
           tCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCg                                                                                                                                       <  1:4473104/151‑1 (MQ=255)
           tCCGGCACGGGGCGGAGGATGTCTTCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCgg                                                                                                                                      <  2:209709/151‑1 (MQ=255)
             cGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGa                                                                                                                                     >  2:10064904/1‑151 (MQ=255)
             cGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGa                                                                                                                                     >  2:2792131/1‑151 (MQ=255)
              ggCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCTGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGAt                                                                                                                                    <  1:7847977/151‑1 (MQ=255)
                  cGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTc                                                                                                                                >  1:8744287/1‑151 (MQ=255)
                  cGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTc                                                                                                                                >  1:6825638/1‑151 (MQ=255)
                     ggCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCtt                                                                                                                             >  2:8396570/1‑151 (MQ=255)
                       cGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGAGTGTCCGTCa                                                                                                                           >  2:9218357/1‑151 (MQ=255)
                        ggaggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAg                                                                                                                          <  2:5405857/151‑1 (MQ=255)
                        ggaggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAg                                                                                                                          <  2:3552380/151‑1 (MQ=255)
                        ggaggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAg                                                                                                                          <  2:941694/151‑1 (MQ=255)
                          aggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAggt                                                                                                                        <  2:5916247/151‑1 (MQ=255)
                          aggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAggt                                                                                                                        <  2:10300831/151‑1 (MQ=255)
                          aggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAggt                                                                                                                        <  1:3285587/151‑1 (MQ=255)
                           ggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGTGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAggtg                                                                                                                       <  2:2881645/151‑1 (MQ=255)
                           ggaTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAggtg                                                                                                                       >  2:107299/1‑151 (MQ=255)
                             aTGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAggtggt                                                                                                                     >  1:5138005/1‑151 (MQ=255)
                              tGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGtt                                                                                                                    >  1:10584535/1‑151 (MQ=255)
                              tGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGtt                                                                                                                    <  2:3374676/151‑1 (MQ=255)
                                 cTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGag                                                                                                                 <  1:1808813/151‑1 (MQ=255)
                                          gcagggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGcctt                                                                                                        <  2:1502514/151‑1 (MQ=255)
                                          gcagggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGcctt                                                                                                        <  1:6383315/151‑1 (MQ=255)
                                           cagggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttc                                                                                                       <  1:5549719/151‑1 (MQ=255)
                                             gggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcag                                                                                                     <  1:10656869/151‑1 (MQ=255)
                                             gggcaggCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcag                                                                                                     <  2:2708343/151‑1 (MQ=255)
                                                   gCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtcctt                                                                                               >  1:143831/1‑151 (MQ=255)
                                                     tCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttca                                                                                             >  2:564317/1‑151 (MQ=255)
                                                       gCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGccttcaggtccttcagg                                                                                           >  2:1208331/1‑151 (MQ=255)
                                                          gaagaaGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCg                                                                                        >  1:9593066/1‑151 (MQ=255)
                                                           aagaagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCGTCAGGTCCTTCAGGTCGa                                                                                       >  1:1236529/1‑151 (MQ=255)
                                                              aagGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagca                                                                                    <  1:4635420/151‑1 (MQ=255)
                                                               agGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcag                                                                                   <  2:8871816/151‑1 (MQ=255)
                                                                ggCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcagc                                                                                  <  2:7851658/151‑1 (MQ=255)
                                                                  ccTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcagcag                                                                                <  1:7671228/151‑1 (MQ=255)
                                                                  ccTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGagcagcag                                                                                <  2:1676211/151‑1 (MQ=255)
                                                                       gcgcCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGc                                                                           >  2:9741686/1‑151 (MQ=255)
                                                                       gcgcCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGc                                                                           <  1:6699941/151‑1 (MQ=255)
                                                                        cgcCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGcc                                                                          >  2:6102189/1‑151 (MQ=255)
                                                                        cgcCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGcc                                                                          >  1:4661371/1‑151 (MQ=255)
                                                                         gcCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCt                                                                         >  2:8107580/1‑151 (MQ=255)
                                                                         gcCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGAGTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCt                                                                         >  1:2770112/1‑151 (MQ=255)
                                                                             tgctgcTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTccc                                                                     <  1:3180892/151‑1 (MQ=255)
                                                                              gctgctCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCg                                                                    >  2:9030347/1‑151 (MQ=255)
                                                                                tgctCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGt                                                                  <  1:9846362/151‑1 (MQ=255)
                                                                                 gctCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTc                                                                 >  1:5899500/1‑151 (MQ=255)
                                                                                 gctCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTc                                                                 <  2:9474691/151‑1 (MQ=255)
                                                                                    ccAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTGCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGc                                                              >  1:693909/1‑151 (MQ=255)
                                                                                      agagTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCa                                                            <  1:7127705/151‑1 (MQ=255)
                                                                                           ggCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTgg                                                       <  2:3158600/151‑1 (MQ=255)
                                                                                             ccAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGaa                                                     <  1:4726335/151‑1 (MQ=255)
                                                                                               aGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAgg                                                   <  1:9238606/151‑1 (MQ=255)
                                                                                                   ggcgCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAggcggc                                               <  2:10565271/151‑1 (MQ=255)
                                                                                                   ggcgCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCGCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGCAggcggc                                               <  1:8520209/151‑1 (MQ=255)
                                                                                                   ggcgCAGCAGCGCCTGCCCGGAAGCGCCGGAGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAggcggc                                               <  1:9067657/151‑1 (MQ=255)
                                                                                                    gcgcAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGgcggcg                                              <  2:6852682/151‑1 (MQ=255)
                                                                                                    gcgcAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGgcggcg                                              <  1:852767/151‑1 (MQ=255)
                                                                                                       cagcagCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAAc                                           <  1:891243/151‑1 (MQ=255)
                                                                                                         gccgcgcCTGCCCGGCAGCGCCGGAGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAg                                         <  2:5848354/148‑1 (MQ=255)
                                                                                                          cagcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAgc                                        <  1:2536447/151‑1 (MQ=255)
                                                                                                           agcTCCTGCCCGGAAGCGCCGGTGACCAGCGCCCGGTCCTCGCCGATGCTCTCCGGAGTGTCCTACAGGTGGATGAGCAGCGCCTTCAGGTCCTTCACGACGAGCAGCAGCAGGCCTTCCCCGTCGGCCACCTGGAAGGCGGCGAACAgcg                                       <  2:5168110/151‑1 (MQ=255)
                                                                                                           agcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAgcg                                       >  2:4882432/1‑151 (MQ=255)
                                                                                                           agcGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTACAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAgag                                       >  2:8544179/1‑149 (MQ=255)
                                                                                                            gcgcCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAgcgc                                      <  2:1782590/151‑1 (MQ=255)
                                                                                                              gcCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGACCAGCGCgg                                    <  1:3594702/151‑1 (MQ=255)
                                                                                                                cTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCGTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGcc                                  >  2:4141188/1‑151 (MQ=255)
                                                                                                                 tGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCct                                 >  1:3907262/1‑151 (MQ=255)
                                                                                                                  gCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCctg                                >  1:9469862/1‑151 (MQ=255)
                                                                                                                  gCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCctg                                <  1:10368141/151‑1 (MQ=255)
                                                                                                                   cccGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCctgc                               <  1:1001300/151‑1 (MQ=255)
                                                                                                                    ccGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCctgct                              >  2:7487915/1‑151 (MQ=255)
                                                                                                                     cGGAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCctgctg                             >  2:98206/1‑151 (MQ=255)
                                                                                                                      ggAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTgg                            >  2:4596282/1‑151 (MQ=255)
                                                                                                                       gAAGCGCCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGc                           <  2:2935667/151‑1 (MQ=255)
                                                                                                                          gcgcCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTAt                        >  2:1345449/1‑151 (MQ=255)
                                                                                                                           cgcCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATc                       <  1:441712/151‑1 (MQ=255)
                                                                                                                            gcCGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCg                      <  2:239887/151‑1 (MQ=255)
                                                                                                                              cGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGt                    <  1:5627176/151‑1 (MQ=255)
                                                                                                                              cGGTGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGAGTGTCCTTCAGGGGGTTGAGCAGCGCCTTCAGCTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGt                    <  2:6436798/151‑1 (MQ=255)
                                                                                                                                 tGAACAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCGTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCaa                 <  2:7281546/151‑1 (MQ=255)
                                                                                                                                     cAGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTcc             >  1:6197484/1‑151 (MQ=255)
                                                                                                                                      aGCGCCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCGTCAGGTCCTGCAGGTCGAGCAGCAGCTGGCCTTCGCGGTCGGCCACCTGGCAGGCGGCGAAGCGCGCGGCCTGGTGGCTATCGGTCTACTccg            >  2:42689/1‑150 (MQ=255)
                                                                                                                                         gcCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCagga         >  2:548225/1‑151 (MQ=255)
                                                                                                                                         gcCCGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCagga         >  1:9930158/1‑151 (MQ=255)
                                                                                                                                           ccGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCGaggagg       <  1:2065181/151‑1 (MQ=255)
                                                                                                                                           ccGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCaggagg       <  1:7474306/151‑1 (MQ=255)
                                                                                                                                           ccGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCaggagg       <  2:7027633/151‑1 (MQ=255)
                                                                                                                                           ccGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCaggagg       >  2:3803711/1‑151 (MQ=255)
                                                                                                                                           acGGTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGATGGCTATCGGTCAACTCCaggagg       >  2:9121562/2‑151 (MQ=255)
                                                                                                                                             ggTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGtt     >  1:8161777/1‑151 (MQ=255)
                                                                                                                                              gTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTg    <  2:4483734/151‑1 (MQ=255)
                                                                                                                                              gTCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTg    >  1:2466/1‑151 (MQ=255)
                                                                                                                                               tCCTCGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTgg   >  1:8241836/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GCCCCCGCCCGTCCGGCTCGGGGCGGAGGATGTCTTCCAGTTGCAGGGCAGGCTCGCCGAAGAAGGCCTCGGCGCCCTGCTGCTCCAGAGTGGCCAGGCGGCGCAGCAGCGCCTGCCCGGAAGCGCCGGTGAACAGCGCCCGGTCCTTGCCGAGGATCTCCGGATTGTCCTTCAGGTGGTTGAGCAGCGCCTTCAGGTCCTTCAGGTCGAGCAGCAGCAGGCCTTCCCGGTCGGCCACCTGGAAGGCGGCGAACAGCGCGGCCTGCTGGCTATCGGTCAACTCCAGGAGGTTGGC  >  NC_008463/4997639‑4997933

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: