Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 3,114,336 A→G intergenic (+69/‑63) PA14_RS14235 → / → PA14_RS14240 glucose/quinate/shikimate family membrane‑bound PQQ‑dependent dehydrogenase/TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,114,3360AG97.9% 896.3 / ‑5.8 294intergenic (+69/‑63)PA14_RS14235/PA14_RS14240glucose/quinate/shikimate family membrane‑bound PQQ‑dependent dehydrogenase/TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base G (100/182);  minor base T (6/0);  total (112/182)
Fisher's exact test for biased strand distribution p-value = 2.27e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 298 total aligned reads displayed.

TGATCGTCGCTGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAAGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCACCTGTTGCCCGGCTCCCTCGCGCTCTGCGGCCTGGCCATCGCCTGCGCCG  >  NC_008463/3114188‑3114483
                                                                                                                                                     |                                                                                                                                                   
tGATCGTCGCTGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAggg                                                                                                                                                   <  1:10231510/151‑1 (MQ=255)
tGATCGTCGCTGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAggg                                                                                                                                                   >  2:8426862/1‑151 (MQ=255)
  aTCGTCGCTGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGa                                                                                                                                                 <  2:9668277/151‑1 (MQ=255)
   tcgtctCTGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCAACGCCAACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAg                                                                                                                                                <  1:4255247/151‑1 (MQ=255)
   tcgtcgCTGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAAGGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAg                                                                                                                                                <  1:9838561/151‑1 (MQ=255)
    cgtcgCTGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAgg                                                                                                                                               <  2:1028544/151‑1 (MQ=255)
    cgtcgCTGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAgg                                                                                                                                               <  1:1799477/151‑1 (MQ=255)
     gtcgCTGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCCTCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGgc                                                                                                                                              <  1:2088451/151‑1 (MQ=255)
     gtcgCTGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAATGGGAGgc                                                                                                                                              >  1:3466710/1‑151 (MQ=255)
     gtcgCTGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGgc                                                                                                                                              >  2:10047120/1‑151 (MQ=255)
     gtcgCTGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGgc                                                                                                                                              <  1:9630209/151‑1 (MQ=255)
      tcgCTGGCGGCCACGGCTCCTGCGGCACCCGCATGGGCGACTACATCACCGCCTACGCCCGGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGgcg                                                                                                                                             <  1:7411409/151‑1 (MQ=255)
      tcgCTGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGgcg                                                                                                                                             >  2:8804146/1‑151 (MQ=255)
       cgCTGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGgcgc                                                                                                                                            >  2:2999821/1‑151 (MQ=255)
          tGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGtt                                                                                                                                         <  2:5291283/151‑1 (MQ=255)
           ggCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGGGTTGCCAAGGGGGGGCGCGGTc                                                                                                                                        >  2:4508710/1‑151 (MQ=255)
           ggCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTc                                                                                                                                        <  2:6566392/151‑1 (MQ=255)
           cgCGGACCCGGCTCCTACGGCACCCGCAGGGGCGACTACACCCTCGCCTACGCCCTGCCCAGGCAATGACCCGCCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTc                                                                                                                                        <  1:4779283/150‑1 (MQ=255)
             cGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCCAGGGGCGGCGCGTTccc                                                                                                                                      >  2:5668868/1‑151 (MQ=255)
             cGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGCGGcgcggtacc                                                                                                                                      >  2:2728703/1‑146 (MQ=255)
             cGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTccc                                                                                                                                      <  2:2460282/151‑1 (MQ=255)
              ggACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGTCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCg                                                                                                                                     <  2:9210782/151‑1 (MQ=255)
              ggACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCg                                                                                                                                     <  1:6227351/151‑1 (MQ=255)
               gACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCGTTTGCAAGGGGCGGCGCGGTCCCGc                                                                                                                                    >  1:4331655/1‑151 (MQ=255)
               gACACGCCTCCTTCGGCACCCGCAGGGGCGACTACAACAACGCCCCCGCCCTGCCCCGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCGCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGc                                                                                                                                    <  1:3611407/151‑1 (MQ=255)
                     gCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTgg                                                                                                                              <  2:4432229/151‑1 (MQ=255)
                     gCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTgg                                                                                                                              <  2:2121391/151‑1 (MQ=255)
                       tCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCCCCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCa                                                                                                                            <  1:6043839/151‑1 (MQ=255)
                       tCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGGGTTGGCAAGGGGGGGGGCGTTCCCGcccggggga                                                                                                                            >  1:2434823/1‑145 (MQ=255)
                       tCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCa                                                                                                                            <  1:4940752/151‑1 (MQ=255)
                          gtCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAg                                                                                                                         <  2:7918218/150‑1 (MQ=255)
                           tCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGc                                                                                                                        <  2:6850098/151‑1 (MQ=255)
                            cGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGcc                                                                                                                       <  2:10406463/151‑1 (MQ=255)
                              gCACCCGCATGGGCGACTACATCATCGCCCACGCCCTGCCCAGGCAATGCCCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCATCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGcccc                                                                                                                     <  1:165279/151‑1 (MQ=255)
                               cACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGGGTTTGCCAGGGGGGGCGCGGTGCCCGCCATGGGAAAGcccca                                                                                                                    >  1:10001599/1‑150 (MQ=255)
                                aCCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGGGTTTCCCAGGGGAGGCGCGTTACCCCCCATGGCAAAGccccca                                                                                                                   >  1:1637606/1‑149 (MQ=255)
                                 cccGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACCGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTc                                                                                                                  <  2:4367936/151‑1 (MQ=255)
                                  ccGCATGGGCGACTCCATCCTCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCg                                                                                                                 <  2:7134117/151‑1 (MQ=255)
                                  ccGCAGGGGCGACTACATCATCGCCAACGCCCTGCCCAGGCACTGACCCGTCCACCGGGAGCCCACAGCCGGCACCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGGTCCCGCCCTTGGCAAAGCCCCTTCg                                                                                                                 <  1:9743273/151‑1 (MQ=255)
                                     cATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTTGCAAGGGGGGGCGCGGTGCCGCCCATGGCAAAGCCCCATAGAtt                                                                                                              >  1:9439378/1‑150 (MQ=255)
                                      aTGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGGTTTTCCAAGGGGTGGCGCGGTCCCGCCCATGGCCAAGCCCCGGCGATAc                                                                                                             >  2:59533/1‑151 (MQ=255)
                                      aTGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCGGTGCCCACGGGAGGCGCGTTCCCGCCCATGGCCAAGCCCCGTCgctcc                                                                                                             >  2:8144263/1‑147 (MQ=255)
                                       tGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCCAGGGGAGGCGCGGTCCCGCCCGTGGCACAGCCCCTTCGATAcc                                                                                                            >  2:10600031/1‑151 (MQ=255)
                                       tGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCGTGGCAAAGCCCCTTCGATAcc                                                                                                            >  2:7718024/1‑151 (MQ=255)
                                          gCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACCGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTg                                                                                                         <  2:3130578/151‑1 (MQ=255)
                                           cGACTCCATCAGCGCCTACGCCCTGCCCCGGCAATGACCCGTCCACCGGGAGCCCACCGCCGGCTCCCGGCCTCCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTCGc                                                                                                        <  1:1301467/151‑1 (MQ=255)
                                           cGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTTCCAAGGGGAGGCCCGTTCCCCCCCATGGGCAAGCCCCCTAGATAcaaagc                                                                                                        >  2:4124532/1‑146 (MQ=255)
                                           cGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGc                                                                                                        <  1:5842782/151‑1 (MQ=255)
                                           cGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGGTCCCGCCCTTGGCCAAGCCCCGTCGCTAACATGc                                                                                                        >  2:8312415/1‑151 (MQ=255)
                                            gACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCg                                                                                                       >  2:8021202/1‑151 (MQ=255)
                                            gACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGACTTCCGGGCTTCGGCCCGGAAGTGGATTTGCCAAGGGGAGGCGCGTTACCGCCCCTGTCAAAGACCCATCGATACCATGCg                                                                                                       >  2:612328/1‑151 (MQ=255)
                                            gACTACATCAGCGCCAACGCCCGGCCCAGCCAATGACCCCTCCCCCGGGAGCCCACACCCGGCTCCCGGCCTTCCGGGCTGCGGCCCGGAAGGGGCTTCGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCATCGAGACCTTGCg                                                                                                       <  2:10744754/151‑1 (MQ=255)
                                            gACTACACCATCGCCGACCCCCTGCCCAGGCACTGCCCCGCCCACCGGGACCCCACCGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCGTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCg                                                                                                       <  2:5030258/151‑1 (MQ=255)
                                               tACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGGGTTGGCAAGGGGGGGCGCGGTCCCCCCCATGGCAAAGCCCCGACGATAACGTGCGagc                                                                                                    >  1:5189205/1‑149 (MQ=255)
                                                   tcatcGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCcgcg                                                                                                <  1:8796253/151‑1 (MQ=255)
                                                   tcatcGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCcgcg                                                                                                <  1:7884695/151‑1 (MQ=255)
                                                    catcGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCGTTTCCAAGGGGAGGCGCGGTCCCGCCCTTGGCAAAGACCCTTCGATTACATGCGAGACGCgg                                                                                               >  2:1992897/1‑151 (MQ=255)
                                                          cTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGGCAAGGGGAGGCGCGGGGCCGCCCGTGGGCAAGCCCCTTCGATAACGTGCGATCCGCGGAAcccg                                                                                         >  2:1105312/1‑150 (MQ=255)
                                                           tACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTg                                                                                        <  1:1599778/151‑1 (MQ=255)
                                                           tACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTg                                                                                        <  2:2314646/151‑1 (MQ=255)
                                                              gCCCTGCCCCGGCAATGACCCGTCCCCCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGACGGGGCTTTGCCATGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCAt                                                                                     <  1:4610007/151‑1 (MQ=255)
                                                              gCCCTGCCCAGGCATTGACCCGTCCACCGGGAGCCCACCGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGTGGCGCGTTCCCGCCCTTGGCAATGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCAt                                                                                     <  1:1153721/151‑1 (MQ=255)
                                                              gCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAATGGGAGGCGCGTTCCCGCCCTTGGCAAATCCCCGTAGATACCTTGCGATCCGCGGAACCCAGCAt                                                                                     >  2:6298575/1‑151 (MQ=255)
                                                                ccTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCATTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGc                                                                                   <  2:9508426/151‑1 (MQ=255)
                                                                 cggCCCAGGCAATGCCCCGTCCCCCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGcc                                                                                  <  1:10065168/149‑1 (MQ=255)
                                                                  tGCCCAGGCATTGACCCGTCCACCGGGAGCCCACAGCGGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCTAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCg                                                                                 <  2:9991072/151‑1 (MQ=255)
                                                                   gCCCTGTCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGt                                                                                <  1:10456874/151‑1 (MQ=255)
                                                                     ccAGGCAATGCCCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTcc                                                                              <  1:8493209/151‑1 (MQ=255)
                                                                     ccAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTcc                                                                              <  1:10662379/151‑1 (MQ=255)
                                                                     ccAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCGTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTcc                                                                              >  2:9739772/1‑151 (MQ=255)
                                                                     ccAGGCAATGACCCGCCCACCGGGAGCCCACCGCCGGCTCCCGGCCTTCCGGGCTACGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTcc                                                                              <  1:6789594/151‑1 (MQ=255)
                                                                          cAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCCAGGGGCGGCGCGGTCCCGCCCGTGGCCAAGCCCCCTCGAGACCTTGCGATCCGCGGCACCCAGCAAGCCGTCCCGAcc                                                                         >  2:6009577/1‑151 (MQ=255)
                                                                          cAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAATGGGAGGCGCGTTCCCGCCCTTGGCCAAGCGCCTTAGCTACCTTGCGATCCGCGGAACCCAGCATGCCGTCCCGAcc                                                                         >  2:10534055/1‑151 (MQ=255)
                                                                             tGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCg                                                                      <  1:9502059/151‑1 (MQ=255)
                                                                                  cggccACCGGGAGCCCACAGCCGGCTCCCGCCCTCCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCa                                                                 <  1:2750180/148‑1 (MQ=255)
                                                                                  cGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGGTTTGCCCAGGGGCGGCGGGTTCCCGCCCGGGGGCCAGCCCCGTCGCGTCCAAGCGCTCCGCGGACCCCGGCATGCCGGCCCGACCCCCGCCCa                                                                 >  2:7232786/1‑151 (MQ=255)
                                                                                  cGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCa                                                                 <  2:8644808/151‑1 (MQ=255)
                                                                                   gTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGGTTTTCCAAGGGGGGGCGCGTTCCCGCCCTTGGCAAAACCCCCTCGCTACCATGCGATCCCCGGGACCCAGCAAGCCGTGCCCACCCCCTccccg                                                                >  1:8103696/1‑149 (MQ=255)
                                                                                    tCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGc                                                               <  2:8194500/151‑1 (MQ=255)
                                                                                               cccACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGGTTTTCCAAGGGGGGGGGCGTTCCCCCCCGTGGGAAACCCCCAGCGATACCATGCGATCCGCGGAACCCAGCATGGCGTCCCCACCCCGGCCCCGCGTCGTCACCCa                                                    >  1:7041428/1‑151 (MQ=255)
                                                                                               cccACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTGGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCa                                                    <  1:3094514/151‑1 (MQ=255)
                                                                                                ccACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCATGGGGAGGCGCGTTCCCGCCCTGGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCAc                                                   <  2:4844298/151‑1 (MQ=255)
                                                                                                  acaGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCGTTGCCAAGGGGAGGAGCGGTCCCGCCCGTGGCCAACCCCCGTAGATACCGTGCGATCGGCGGAACCCTGCCTGGCGTCCCCACCCCGGCCCCGCTTCGGCTCCCAcca                                                 >  2:3056942/1‑150 (MQ=255)
                                                                                                         gCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCACGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCATCGATACCTGGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTACGTCTCCCACCAGTTGccc                                          <  2:1649182/151‑1 (MQ=255)
                                                                                                         gCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCACCTGTTGccc                                          <  2:3655999/151‑1 (MQ=255)
                                                                                                                 ccTTCCGGGCTTCGGCCCGGAAGGGGCATTTCCCAAGGGAGGCGCGGTCCCGCCCATTGCAAAAACCCATAGAAACCTTGCGAACCCCGGAACCCAGCATGACCTCCCGACCCCGACCCAGCTGCCTCACCCCCCAGGTTCCCGGCGcccc                                  >  2:1439804/1‑150 (MQ=255)
                                                                                                                  cTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCACCTGTTGCCCGGCTCCCTc                                 <  1:10590718/151‑1 (MQ=255)
                                                                                                                  cTTCCGGGCTTCGGCCCGGAAGGGGCGTTGCCCAGGGGAGGCGCCTTCCCCCCCATGGCAAAGCCCCGACGATACCCTGCCATCCCCGGAACCCAGCAAGCCGTCCCCACCCCCTCCCAGCTGCGTCTCCCACCAGTTGCCCCGCAcccac                                 >  2:4014786/1‑149 (MQ=255)
                                                                                                                      cGGGCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCACCTGTTGCCCGGCTCCCTCgcgc                             <  1:8388548/151‑1 (MQ=255)
                                                                                                                       gggCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCACCTGTTGCCCGGCTCCCTCGCGct                            <  1:3942557/151‑1 (MQ=255)
                                                                                                                         gCTTCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTGCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCACCTGTTGCCCGGCTCCCTCGCGctct                          <  1:10811470/151‑1 (MQ=255)
                                                                                                                           ttCGGCCCGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCACCTGTTGCCCGGCTCCCTCGCGCTCTGc                        <  1:7571952/151‑1 (MQ=255)
                                                                                                                                 ccGGAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCACCTGTTGCCCGGCTCCCTCGCGCTCTGCGGCCTg                  <  1:3357180/151‑1 (MQ=255)
                                                                                                                                    gAAGGGGCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCACCTGTTGCCCGGCTCCCTCGCGCTCTGCGGCCTGGcc               <  1:8247877/151‑1 (MQ=255)
                                                                                                                                    gAAGGGGCTTTGCCAAGGGGAGGCGCGGTCCCGCCCTGGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCACCTGTTGCCCGGCTCCCTCGCGCTCTGCGGCCTGGcc               <  1:9066598/151‑1 (MQ=255)
                                                                                                                                        ggggCTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCACCTGTTGCCCGGCTCCCTCGCGCTCTGCGGCCTGGCCATc            <  2:2269878/151‑1 (MQ=255)
                                                                                                                                          ggCTTTGCCAAGGGGAGGCGCGTGCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCACCTGTGGCCCGGCTCCCTCGCGCTCTGCGGCCTGGCCATCGc          <  1:4432409/151‑1 (MQ=255)
                                                                                                                                          cgCTTTGCCAAGGGGAGGCGCGTACCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCCTGCCGACCCGACCCCGTCACAGCTTCGTCTCCCACCTGATGCCCGGCTCCCTCGCGCTCTGCGGCCTGGCCAACGc          >  2:3402834/2‑151 (MQ=255)
                                                                                                                                            cTTTGCCAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCACCTGTTGCCCGGCTCCCTCGCGCTCTGCGGCCTGGCCATCGCCt        >  1:7356028/1‑151 (MQ=255)
                                                                                                                                                  cAAGGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCACCTGTTGCCCGGCTCCCTCGCGCTCTGCGGCCTGGCCATCGCCTGCGCCg  >  2:546934/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TGATCGTCGCTGGCGGACACGGCTCCTTCGGCACCCGCATGGGCGACTACATCATCGCCTACGCCCTGCCCAGGCAATGACCCGTCCACCGGGAGCCCACAGCCGGCTCCCGGCCTTCCGGGCTTCGGCCCGGAAGGGGCTTTGCCAAAGGGAGGCGCGTTCCCGCCCTTGGCAAAGCCCCTTCGATACCTTGCGATCCGCGGAACCCTGCATGCCGTCCCGACCCCGTCCCAGCTTCGTCTCCCACCTGTTGCCCGGCTCCCTCGCGCTCTGCGGCCTGGCCATCGCCTGCGCCG  >  NC_008463/3114188‑3114483

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