Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,028,201 2 bp→CG intergenic (‑123/+13) 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,2010GC98.8% 1636.9 / ‑1.0 428intergenic (‑123/+14)PA14_RS22975/PA14_RS22980Nramp family divalent metal transporter/ribbon‑helix‑helix domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base C (193/229);  total (198/230)
Fisher's exact test for biased strand distribution p-value = 1.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0084635,028,2020TG99.1% 1665.2 / 4.2 428intergenic (‑124/+13)PA14_RS22975/PA14_RS22980Nramp family divalent metal transporter/ribbon‑helix‑helix domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base G (195/229);  minor base . (3/1);  total (198/230)
Fisher's exact test for biased strand distribution p-value = 3.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 428 total aligned reads displayed.

GCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGGTGGGTTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTCC  >  NC_008463/5028059‑5028338
                                                                                                                                              ||                                                                                                                                        
gCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTtccctc                                                                                                                                 <  1:9714988/151‑1 (MQ=255)
    cggccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTccctcccct                                                                                                                             >  2:10157945/1‑151 (MQ=255)
     ggccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTccctcccctc                                                                                                                            <  1:5386548/151‑1 (MQ=255)
     ggccggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTccctcccctc                                                                                                                            <  2:1888916/151‑1 (MQ=255)
        cggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAgg                                                                                                                         >  2:7328083/1‑151 (MQ=255)
         ggcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGc                                                                                                                        <  2:6450378/151‑1 (MQ=255)
          gcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGcc                                                                                                                       >  1:610679/1‑151 (MQ=255)
          gcACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGcc                                                                                                                       <  1:3071376/151‑1 (MQ=255)
           acgcGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCGCCCCTCAGGCCt                                                                                                                      <  2:6777224/148‑1 (MQ=255)
            aCCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTg                                                                                                                     >  2:9066536/1‑151 (MQ=255)
            aCCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTg                                                                                                                     <  1:4313033/151‑1 (MQ=255)
            aCCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTg                                                                                                                     >  2:6123993/1‑151 (MQ=255)
              cGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTggg                                                                                                                   >  1:9282354/1‑151 (MQ=255)
                 cACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgc                                                                                                                <  1:4111784/151‑1 (MQ=255)
                 cACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGAGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgc                                                                                                                >  2:2230604/1‑151 (MQ=255)
                  aCGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcg                                                                                                               >  1:5142908/1‑151 (MQ=255)
                   cGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGgcgcgc                                                                                                              >  1:3318954/1‑151 (MQ=255)
                    gCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAACCGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCt                                                                                                             <  2:3323663/151‑1 (MQ=255)
                    gCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCt                                                                                                             <  2:887226/151‑1 (MQ=255)
                       gtgGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGc                                                                                                          >  1:2582061/1‑151 (MQ=255)
                           cTGACCGCCGAGGGACAGAACGGCGCGGTGGCCGGGGTGGCCAGGCAGTAATTCACGAGGCGACGCCGGATCGGGGAGGCACAGCGACCATGCACAACGTTTCACCACTGTCGAGTCGTCCCTCCCCTCAGGCCTGGCCGCGCTGGCCGAc                                                                                                      >  1:2896755/1‑151 (MQ=255)
                            tGACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCGCCCCTCAGGCCTGGGCGCGCTGGCCGAcc                                                                                                     >  2:1520136/1‑151 (MQ=255)
                             gACCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCa                                                                                                    <  1:1597442/151‑1 (MQ=255)
                              aCCGCCGAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAg                                                                                                   <  1:2609056/151‑1 (MQ=255)
                                  ccgAGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTc                                                                                               >  1:9570683/1‑151 (MQ=255)
                                    gagccacacAACGCCGCGGTGGCAGTGGAGGGCAGTGACTATATCACGAGGCGACGCCGGATCGCGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGACGCGTCTACCCTCCCCTCAGGCCTGTGCGCGCTGGCCGACCAGGCTCAg                                                                                             >  2:7037396/6‑151 (MQ=255)
                                     aGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGc                                                                                            <  2:5496949/151‑1 (MQ=255)
                                     aGGGACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGc                                                                                            >  2:6586546/1‑151 (MQ=255)
                                       ggACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTc                                                                                          <  2:9548834/151‑1 (MQ=255)
                                       ggACAGAACGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTc                                                                                          >  2:8535420/1‑151 (MQ=255)
                                             aaCGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGAc                                                                                    >  2:5304194/1‑151 (MQ=255)
                                               cGGCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCt                                                                                  >  1:10023974/1‑151 (MQ=255)
                                               cGGCGCGGTGGCCGGGGTGGGCAGGGAGTCTTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTGTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGGCGACGAGGCTCTGCGCGGCGTCCt                                                                                  >  1:4903648/1‑151 (MQ=255)
                                                ggCGCGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCtt                                                                                 >  1:379215/1‑151 (MQ=255)
                                                 gcgcGGTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTtc                                                                                <  2:2667951/151‑1 (MQ=255)
                                                     ggTGGCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcgc                                                                            >  1:9361605/1‑151 (MQ=255)
                                                       gtggccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcgcg                                                                           <  1:7579349/147‑1 (MQ=255)
                                                       gtggccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcgcg                                                                           <  2:10362327/147‑1 (MQ=255)
                                                       ggggccGTGGGGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGGTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTcgcg                                                                           <  2:10604642/147‑1 (MQ=255)
                                                         ggccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGcgcc                                                                         <  2:3422363/149‑1 (MQ=255)
                                                         ggccGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGcgcc                                                                         >  2:5581971/3‑151 (MQ=255)
                                                            cGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCt                                                                      <  1:8762467/151‑1 (MQ=255)
                                                               ggtggGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTgcg                                                                   <  1:9892970/151‑1 (MQ=255)
                                                                   ggCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAg                                                               <  1:1560801/151‑1 (MQ=255)
                                                                     cAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGt                                                             <  1:649922/151‑1 (MQ=255)
                                                                      aGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTa                                                            >  2:9226755/1‑151 (MQ=255)
                                                                       gggAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAg                                                           >  1:8231577/1‑151 (MQ=255)
                                                                       gggAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAg                                                           <  1:2165478/151‑1 (MQ=255)
                                                                        ggAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGc                                                          >  2:9680324/1‑151 (MQ=255)
                                                                         gAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCg                                                         >  2:4424722/1‑151 (MQ=255)
                                                                         gAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGAAGCg                                                         >  2:8416695/1‑151 (MQ=255)
                                                                          aGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCgg                                                        <  1:3029523/151‑1 (MQ=255)
                                                                          aGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCgg                                                        <  2:5386262/151‑1 (MQ=255)
                                                                           gTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCggg                                                       <  2:4768033/151‑1 (MQ=255)
                                                                           gTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCggg                                                       <  2:6870214/151‑1 (MQ=255)
                                                                            tATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGt                                                      <  1:2474606/151‑1 (MQ=255)
                                                                                 cACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAAc                                                 >  2:8594063/1‑151 (MQ=255)
                                                                                 cACGACGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTTGCGGGTGCAAc                                                 <  1:1428225/151‑1 (MQ=255)
                                                                                     aGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGac                                             <  2:3705892/151‑1 (MQ=255)
                                                                                      ggCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGaca                                            >  2:2929571/1‑151 (MQ=255)
                                                                                        cGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACg                                          <  1:8650758/151‑1 (MQ=255)
                                                                                        cGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACg                                          <  1:7498480/151‑1 (MQ=255)
                                                                                        cGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACg                                          <  2:5874496/151‑1 (MQ=255)
                                                                                            gCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCagg                                      <  1:842003/151‑1 (MQ=255)
                                                                                            gCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCagg                                      <  1:10520717/151‑1 (MQ=255)
                                                                                             ccGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCagga                                     <  2:7864027/151‑1 (MQ=255)
                                                                                              cGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaa                                    <  2:7650681/151‑1 (MQ=255)
                                                                                               ggATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaag                                   >  1:8991826/1‑151 (MQ=255)
                                                                                               ggATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaag                                   >  2:1094341/1‑151 (MQ=255)
                                                                                                gATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaagg                                  >  2:7462507/1‑151 (MQ=255)
                                                                                                 ctCGGGGAGGAAGTGCGACCAGCCGAAACGTTTCACCCCGGTCGCGACTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGANCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaagga                                 <  1:5631461/150‑1 (MQ=255)
                                                                                                 aTCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCaggaagga                                 <  2:1783263/151‑1 (MQ=255)
                                                                                                  tCGGGGAGGAAGTGCGACCAGCCTATACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCTGCTCGGCGACCTTCTCGCGCAGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAaggag                                <  1:10612657/151‑1 (MQ=255)
                                                                                                  tCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAaggag                                <  1:7016814/151‑1 (MQ=255)
                                                                                                   cGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAaggagg                               <  1:10180909/151‑1 (MQ=255)
                                                                                                     gggAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCg                             <  2:4618603/151‑1 (MQ=255)
                                                                                                     gggAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCg                             >  1:10361877/1‑151 (MQ=255)
                                                                                                      ggaggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGa                            <  2:7274229/151‑1 (MQ=255)
                                                                                                      ggaggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGa                            <  2:1486069/151‑1 (MQ=255)
                                                                                                       gaggaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTCGCGGGTGCAACTGACACGCAGGAAGGAGGCGaa                           >  1:8507784/1‑151 (MQ=255)
                                                                                                          gaAGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGtt                        >  2:2353859/1‑151 (MQ=255)
                                                                                                           aaGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTg                       <  2:98206/151‑1 (MQ=255)
                                                                                                            aGTGCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGTAGGAGGCGAAGTTGa                      >  1:4704936/1‑151 (MQ=255)
                                                                                                            aGTGCGACCAGACTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGa                      <  1:968614/151‑1 (MQ=255)
                                                                                                              tGCGACCAGCCTAAACGTTTCACCCCGGACGCGTCTTACCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGc                    <  1:3599571/151‑1 (MQ=255)
                                                                                                               gCGACCAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCa                   <  2:509626/151‑1 (MQ=255)
                                                                                                                   ccAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTc               >  2:4071967/1‑151 (MQ=255)
                                                                                                                   ccAGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTc               <  1:5918152/151‑1 (MQ=255)
                                                                                                                     aGCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgc             >  2:10539374/1‑151 (MQ=255)
                                                                                                                      gCCTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgcc            <  2:8008107/151‑1 (MQ=255)
                                                                                                                       ccTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgccg           >  1:1835485/1‑151 (MQ=255)
                                                                                                                       ccTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgccg           >  2:9269050/1‑151 (MQ=255)
                                                                                                                        cTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTCGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCTCCTcgccgc          >  2:1156974/1‑151 (MQ=255)
                                                                                                                        cTAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTcgccgc          >  2:2037787/1‑151 (MQ=255)
                                                                                                                         tAAACGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCcgcg         <  2:3909197/151‑1 (MQ=255)
                                                                                                                          aaaCGTTTCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGa        <  1:7945259/151‑1 (MQ=255)
                                                                                                                              gTTTCACCCCGGTCGCGTCTTCCCTCCCCTCATGCCTGGGCGCGCGGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGt    >  1:120894/1‑151 (MQ=255)
                                                                                                                               tttCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTc   <  1:7820309/151‑1 (MQ=255)
                                                                                                                               tttCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTc   <  2:4300410/151‑1 (MQ=255)
                                                                                                                                ttCACCCCGGTCGCGTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTcc  <  1:4005204/151‑1 (MQ=255)
                                                                                                                                              ||                                                                                                                                        
GCGTCGGCCGGCACCGCCACGCTGTGGCTGACCGCCGAGGGACAGAACGGCGCGGTGGCCCGTGGTGGGCAGGGAGTATTTCACGAGGCGACGCCGGATCGGGGAGGAAGTGCGACCAGCCTAAACGTTTCACCCCGGGTGGGTTCTTCCCTCCCCTCAGGCCTGGGCGCGCTGGCCGACCAGGCTCAGCTCGGCGACCTTCTCGCGCCGCTGCGCCAGGTAGCGGGTGCAACTGACACGCAGGAAGGAGGCGAAGTTGAGCACCTCGCCGCGATAGTCC  >  NC_008463/5028059‑5028338

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