Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,723,048 Δ1 bp coding (587/1527 nt) PA14_RS16920 → inactive transglutaminase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,723,0480T.88.9% 1518.1 / 213.6 386coding (587/1527 nt)PA14_RS16920inactive transglutaminase family protein
Reads supporting (aligned to +/- strand):  ref base T (23/20);  new base . (153/190);  total (176/210)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 390 total aligned reads displayed.

AGAAGCCCTGCTGGCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCAACGTGGTCACCGGCGAA  >  NC_008463/3722899‑3723193
                                                                                                                                                     |                                                                                                                                                 
aGAAGCCCTGCTGGCCCCCAGTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATAtcgt                                                                                                                                                   <  1:8407404/150‑1 (MQ=255)
    gCCCTGCTGGCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCgaga                                                                                                                                              <  2:8299558/150‑1 (MQ=255)
       ctgctgGCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGAccc                                                                                                                                           >  2:2410338/1‑150 (MQ=255)
       ctgctgGCCCCCAGTCGCCAGCACTCGGCCGCCGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGAccc                                                                                                                                           <  2:5966680/150‑1 (MQ=255)
        tgctgGCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCct                                                                                                                                          <  2:6159554/150‑1 (MQ=255)
          ctgGCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCctgc                                                                                                                                        <  1:10176816/150‑1 (MQ=255)
          ctgGCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCctgc                                                                                                                                        <  1:3729094/150‑1 (MQ=255)
            ggCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCctgctg                                                                                                                                      <  1:10387885/150‑1 (MQ=255)
             gCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGt                                                                                                                                     >  1:6904646/1‑150 (MQ=255)
              cccccATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTc                                                                                                                                    >  2:4629770/1‑150 (MQ=255)
              cccccATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTc                                                                                                                                    <  1:819352/150‑1 (MQ=255)
               ccccATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTcc                                                                                                                                   >  1:8906638/1‑150 (MQ=255)
               ccccATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTcc                                                                                                                                   >  1:1405819/1‑150 (MQ=255)
                 ccATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCAt                                                                                                                                 <  1:4292013/150‑1 (MQ=255)
                      cGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGccc                                                                                                                            <  2:7458397/150‑1 (MQ=255)
                       gCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGACAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCa                                                                                                                           >  2:10269439/1‑150 (MQ=255)
                          aGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACCCCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGt                                                                                                                        >  2:2483206/1‑150 (MQ=255)
                          aGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGt                                                                                                                        >  2:5503859/1‑150 (MQ=255)
                             aCTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGATAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGcc                                                                                                                     <  2:3173839/150‑1 (MQ=255)
                                  gCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATCGCCCACGTGCCGATg                                                                                                                 >  2:10600618/1‑150 (MQ=255)
                                    cGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAg                                                                                                              <  1:4569646/150‑1 (MQ=255)
                                      aCGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCg                                                                                                            <  1:1468577/150‑1 (MQ=255)
                                        gTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCggg                                                                                                          >  1:10021042/1‑150 (MQ=255)
                                        gTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCggg                                                                                                          <  2:6708513/150‑1 (MQ=255)
                                            agaCCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTcc                                                                                                       <  1:2223990/150‑1 (MQ=255)
                                            agaCCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCca                                                                                                      <  1:5244455/150‑1 (MQ=255)
                                             gaCCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCcac                                                                                                     >  2:1645764/1‑150 (MQ=255)
                                             gaCCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCcac                                                                                                     >  2:6995810/1‑150 (MQ=255)
                                                cTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACAcc                                                                                                  >  1:5127830/1‑150 (MQ=255)
                                                   cGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCAt                                                                                                <  1:7918879/150‑1 (MQ=255)
                                                    gTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATcc                                                                                              <  1:8797609/150‑1 (MQ=255)
                                                     tCAGCGAAGCCATCAAGCGCGCCAACACCGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCg                                                                                             <  1:10697839/150‑1 (MQ=255)
                                                       aGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGcc                                                                                           <  1:1989872/150‑1 (MQ=255)
                                                       aGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGcc                                                                                           >  1:90098/1‑150 (MQ=255)
                                                        gCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCt                                                                                          <  2:4903396/150‑1 (MQ=255)
                                                           caGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGt                                                                                       <  2:9472725/149‑1 (MQ=255)
                                                           aaGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGt                                                                                       <  1:1219753/150‑1 (MQ=255)
                                                            aGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTc                                                                                      >  2:3286723/1‑150 (MQ=255)
                                                            aGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTc                                                                                      <  1:7959813/150‑1 (MQ=255)
                                                                aTCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCg                                                                                  <  1:5438654/150‑1 (MQ=255)
                                                                   aaGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGAcc                                                                               >  2:177520/1‑150 (MQ=255)
                                                                    aGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCt                                                                              <  2:6378644/150‑1 (MQ=255)
                                                                    aGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCt                                                                              >  2:2201247/1‑150 (MQ=255)
                                                                      cgcgcCAACACCGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGACAGCGCATATCGCGAGACCCTGCGGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCCTCCGCCTGGTCGCCGACCTGc                                                                            <  2:8861710/150‑1 (MQ=255)
                                                                      cgcgcCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGc                                                                            >  1:9962287/1‑150 (MQ=255)
                                                                          ccAACACCGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACa                                                                        <  2:7676386/150‑1 (MQ=255)
                                                                           caacaaCGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACag                                                                       >  1:7114473/1‑150 (MQ=255)
                                                                           caacaaCGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACag                                                                       >  1:7820812/1‑150 (MQ=255)
                                                                            aacaacGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACaga                                                                      <  1:4111511/150‑1 (MQ=255)
                                                                              caacGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGc                                                                    >  1:3553078/1‑150 (MQ=255)
                                                                                acGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGccc                                                                  <  1:1052229/150‑1 (MQ=255)
                                                                                 cGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGccc                                                                  <  1:9371917/150‑1 (MQ=255)
                                                                                   gCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGa                                                               >  2:7333907/1‑150 (MQ=255)
                                                                                   gCAACGACGACAACGTCAAGCAGCTGATGCCCGGCGACACCAGACTGGAAAAGAAAGCGCATATCGCGAGACCCTGAAGTCCATCGCCAACGAGCCGAAGGATCGGGTCCACACCATCACCCTGCTCGCCGACCTGCAACATAGCACGGa                                                               >  1:7379529/1‑150 (MQ=255)
                                                                                        gacgacAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACtgt                                                          <  2:4393740/150‑1 (MQ=255)
                                                                                        gacgacAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACtgt                                                          <  1:662500/150‑1 (MQ=255)
                                                                                          cgacAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTgg                                                        <  1:3247188/150‑1 (MQ=255)
                                                                                           gacAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGc                                                       <  1:8544385/150‑1 (MQ=255)
                                                                                              aaCGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTgc                                                    <  1:6392919/150‑1 (MQ=255)
                                                                                               aCGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGACCTGTGGCTgcg                                                   >  2:8126672/1‑150 (MQ=255)
                                                                                               aCGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTgcg                                                   <  2:1558905/150‑1 (MQ=255)
                                                                                                cGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCCTCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTgcgc                                                  <  2:7914264/150‑1 (MQ=255)
                                                                                                    aaGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCt                                              <  1:4700957/150‑1 (MQ=255)
                                                                                                     aGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCtt                                             >  1:5726593/1‑150 (MQ=255)
                                                                                                      gctgctGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTc                                            >  2:36460/1‑150 (MQ=255)
                                                                                                       ctgctgATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTc                                            >  2:5994632/1‑150 (MQ=255)
                                                                                                           tgATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACgg                                       <  2:8049968/150‑1 (MQ=255)
                                                                                                           tgATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACgg                                       <  2:2596046/150‑1 (MQ=255)
                                                                                                            gATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACgg                                       >  2:9052609/1‑150 (MQ=255)
                                                                                                                 ccGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGaa                                 >  2:5316037/1‑150 (MQ=255)
                                                                                                                 ccGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGaa                                 >  1:9653926/1‑150 (MQ=255)
                                                                                                                  cGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAAc                                <  2:2978107/150‑1 (MQ=255)
                                                                                                                  cGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAAc                                >  2:5164356/1‑150 (MQ=255)
                                                                                                                    gCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACt                               >  2:9680653/1‑150 (MQ=255)
                                                                                                                    gCGACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTg                              <  1:7608384/150‑1 (MQ=255)
                                                                                                                      gACACCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGc                            <  1:10293631/150‑1 (MQ=255)
                                                                                                                       acacCAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCt                           <  1:1705596/150‑1 (MQ=255)
                                                                                                                          ccAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTa                        <  1:3068213/150‑1 (MQ=255)
                                                                                                                           cAGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTAc                       >  1:6233522/1‑150 (MQ=255)
                                                                                                                            aGCCTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCAACCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACt                      <  2:6605893/150‑1 (MQ=255)
                                                                                                                              ccTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACtt                     >  1:2857834/1‑150 (MQ=255)
                                                                                                                              ccTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGACCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACtt                     >  2:9276375/1‑150 (MQ=255)
                                                                                                                              ccTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTc                    >  2:10378038/1‑150 (MQ=255)
                                                                                                                              ccTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTc                    <  1:6023550/150‑1 (MQ=255)
                                                                                                                              ccTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTc                    <  2:364015/150‑1 (MQ=255)
                                                                                                                               cTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCa                   <  2:4735767/150‑1 (MQ=255)
                                                                                                                               cTGGAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCa                   >  1:9157095/1‑150 (MQ=255)
                                                                                                                                 ggAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCAAc                 <  2:1943243/150‑1 (MQ=255)
                                                                                                                                  gAAAAGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCAACg                >  2:6452304/1‑150 (MQ=255)
                                                                                                                                   aaaaGAAAGCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCAACGt               <  1:477889/150‑1 (MQ=255)
                                                                                                                                    aaagaaagCGCATATCGTCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCAACGt               >  1:9814399/1‑150 (MQ=255)
                                                                                                                                     aagaaagCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCAACGTgg             <  1:2696710/150‑1 (MQ=255)
                                                                                                                                      agaaagCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCAACGTGGt            <  2:5747870/150‑1 (MQ=255)
                                                                                                                                       gaaagCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCAACGTGGTc           <  1:2474102/150‑1 (MQ=255)
                                                                                                                                       gaaagCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGACGAGAACTGGCTTTACTTCAACGTGGTc           <  2:1184867/150‑1 (MQ=255)
                                                                                                                                       gaaagCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGGGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCAACGTGGTc           <  2:1719976/150‑1 (MQ=255)
                                                                                                                                         aagCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCAACGTGGTCAc         >  2:10099524/1‑150 (MQ=255)
                                                                                                                                          agCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCAACGTGGTCAcc        <  1:4925652/150‑1 (MQ=255)
                                                                                                                                          agCGCATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCAACGTGGTCAcc        <  1:4051942/150‑1 (MQ=255)
                                                                                                                                           gcgcATATCGCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCAACGTGGTCACCg       <  2:2854553/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
AGAAGCCCTGCTGGCCCCCATTCGCCAGCACTCGGCCGACGTCGAGACCTTCGTCAGCGAAGCCATCAAGCGCGCCAACAACGGCAACGACGACAACGTCAAGCTGCTGATGGCCGGCGACACCAGCCTGGAAAAGAAAGCGCATATCGTCGAGACCCTGCTGTCCATCGCCCACGTGCCGATGGAGCGGGTCCACACCATCCGCCTGGTCGCCGACCTGCAACAGAGCCCGGAACTGTGGCTGCGCAGCTTCAACGGCGAGAACTGGCTTTACTTCAACGTGGTCACCGGCGAA  >  NC_008463/3722899‑3723193

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: