Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,566,196 A→C 5.5% intergenic (‑313/‑167) PP_1374 ← / → pcaR hypothetical protein/pca regulon transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,566,1960AC 5.5% 385.9 / 8.9 118intergenic (‑313/‑167)PP_1374/pcaRhypothetical protein/pca regulon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (51/60);  new base C (2/5);  total (53/65)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

AATTTACGCAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAACCGGTGT  >  NC_002947/1566048‑1566345
                                                                                                                                                    |                                                                                                                                                     
  tttACGCAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGc                                                                                                                                                    <  2:301325/150‑1 (MQ=255)
  tttACGCAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGc                                                                                                                                                    <  2:1098680/150‑1 (MQ=255)
    tACGCAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGa                                                                                                                                                  <  2:1703029/150‑1 (MQ=255)
      cGCAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAgg                                                                                                                                                <  2:5402859/150‑1 (MQ=255)
       gCAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCTTTCGGCGAGGc                                                                                                                                               <  1:502030/150‑1 (MQ=255)
        cAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCTTTCGGCGAGGcc                                                                                                                                              <  1:4364876/150‑1 (MQ=255)
             cTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGc                                                                                                                                         <  1:1779131/150‑1 (MQ=255)
                aCATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTg                                                                                                                                      <  2:2816823/150‑1 (MQ=255)
                aCATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTg                                                                                                                                      <  1:2202091/150‑1 (MQ=255)
                  atcatcTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTggg                                                                                                                                    >  2:946773/1‑150 (MQ=255)
                  atcatcTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTggg                                                                                                                                    >  1:2420171/1‑150 (MQ=255)
                  atcatcTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTggg                                                                                                                                    >  1:2197831/1‑150 (MQ=255)
                     atcTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGAc                                                                                                                                 <  1:1714846/150‑1 (MQ=255)
                     atcTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGAc                                                                                                                                 >  2:4557504/1‑150 (MQ=255)
                     atcTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGAc                                                                                                                                 >  2:4669351/1‑150 (MQ=255)
                      tcTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACt                                                                                                                                >  2:2079908/1‑150 (MQ=255)
                       cTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACtt                                                                                                                               >  2:4794740/1‑150 (MQ=255)
                       cTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACtt                                                                                                                               >  2:5423320/1‑150 (MQ=255)
                         ggCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCt                                                                                                                             <  2:2788193/150‑1 (MQ=255)
                           cGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTAc                                                                                                                           <  2:3197304/150‑1 (MQ=255)
                            gTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACg                                                                                                                          >  2:1671367/1‑150 (MQ=255)
                             tAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACgg                                                                                                                         <  1:2150236/150‑1 (MQ=255)
                              aaGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGt                                                                                                                        <  1:976429/150‑1 (MQ=255)
                                 tAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTggg                                                                                                                     >  2:2234237/1‑150 (MQ=255)
                                 tAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTggg                                                                                                                     >  2:2234352/1‑150 (MQ=255)
                                   gCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTa                                                                                                                   >  1:3936975/1‑150 (MQ=255)
                                    cAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTAc                                                                                                                  <  1:1329228/150‑1 (MQ=255)
                                       aTTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACg                                                                                                                          <  2:1605059/139‑1 (MQ=255)
                                       aTTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACg                                                                                                                          >  1:1605059/1‑139 (MQ=255)
                                            gCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAg                                                                                                          >  2:2506240/1‑150 (MQ=255)
                                             cccGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGt                                                                                                         <  1:3421458/150‑1 (MQ=255)
                                              ccGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGtc                                                                                                        >  1:2802528/1‑150 (MQ=255)
                                              ccGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGtc                                                                                                        >  2:1099485/1‑150 (MQ=255)
                                               cgcgAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGtct                                                                                                       >  1:1404216/1‑150 (MQ=255)
                                                   aaCAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTcc                                                                                                   <  1:4677196/150‑1 (MQ=255)
                                                    aCAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCa                                                                                                  >  1:4761978/1‑150 (MQ=255)
                                                      aGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCaaa                                                                                                <  1:3241815/150‑1 (MQ=255)
                                                         tGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAaagaa                                                                                             <  2:2023506/150‑1 (MQ=255)
                                                           ggACCTTTCGTACACAGGTTGCCAGGACTGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCACAGAAGc                                                                                           <  2:1298438/150‑1 (MQ=255)
                                                                tttCGTACAAAGTTGGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACaa                                                                                      <  1:4686590/150‑1 (MQ=255)
                                                                  tCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGt                                                                                    >  1:2428775/1‑150 (MQ=255)
                                                                    gTACAAAGTTTGAAAGGACTGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTATGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTgg                                                                                  <  2:3371827/150‑1 (MQ=255)
                                                                     tACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGc                                                                                 >  1:484816/1‑150 (MQ=255)
                                                                      aCAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCa                                                                                <  1:1916361/150‑1 (MQ=255)
                                                                        aaaGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCAt                                                                               <  1:2478792/149‑1 (MQ=255)
                                                                        aaaGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCAt                                                                               >  2:2478792/1‑149 (MQ=255)
                                                                        aaaGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATg                                                                              >  2:1801393/1‑150 (MQ=255)
                                                                        aaaGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATg                                                                              >  2:1010400/1‑150 (MQ=255)
                                                                             ttGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGaaa                                                                         <  2:1975313/150‑1 (MQ=255)
                                                                              tGAAAGGAATGGCGGATAGGAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAg                                                                        <  2:921918/150‑1 (MQ=255)
                                                                              tGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGTAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCTGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAg                                                                        <  2:1765409/150‑1 (MQ=255)
                                                                              tGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAg                                                                        <  1:4146049/150‑1 (MQ=255)
                                                                              tGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAg                                                                        <  1:4277118/150‑1 (MQ=255)
                                                                                aaaGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGt                                                                      >  2:909627/1‑150 (MQ=255)
                                                                                aaaGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGt                                                                      >  2:909557/1‑150 (MQ=255)
                                                                                    gAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATaa                                                                  >  2:4720314/1‑150 (MQ=255)
                                                                                        ggCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCaa                                                              <  2:1965792/150‑1 (MQ=255)
                                                                                            gATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATaa                                                                  <  1:1939668/142‑1 (MQ=255)
                                                                                            gATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATaa                                                                  >  2:1939668/1‑142 (MQ=255)
                                                                                                     tCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGt                                                 <  1:2504596/150‑1 (MQ=255)
                                                                                                       gTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCa                                               >  1:4078950/1‑150 (MQ=255)
                                                                                                       gTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCa                                               >  2:1620752/1‑150 (MQ=255)
                                                                                                        ttGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCaa                                              <  2:3428688/150‑1 (MQ=255)
                                                                                                        ttGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCaa                                              <  2:693426/150‑1 (MQ=255)
                                                                                                         tGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAAt                                             >  1:5249145/1‑150 (MQ=255)
                                                                                                              tGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGa                                        <  1:818347/150‑1 (MQ=255)
                                                                                                                tCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCaaaa                                                            <  1:1291896/128‑1 (MQ=255)
                                                                                                                tCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCaaaa                                                            >  2:1291896/1‑128 (MQ=255)
                                                                                                                 cAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTa                                     <  2:1404216/150‑1 (MQ=255)
                                                                                                                   gggcgggcCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATc                                   >  1:4247940/1‑150 (MQ=255)
                                                                                                                        ggcCTTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCg                              <  1:4901024/150‑1 (MQ=255)
                                                                                                                         tcttggttGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGt                             <  1:1416635/143‑1 (MQ=255)
                                                                                                                          ccGGAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACGTCTACGGTGGGTACCTGACCAGTCTCTCCAACGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGtt                            <  1:3222732/150‑1 (MQ=255)
                                                                                                                               aTGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTc                       <  2:3135307/150‑1 (MQ=255)
                                                                                                                                tGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTACAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCg                      >  2:5175947/1‑150 (MQ=255)
                                                                                                                                 gAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGa                     >  1:2928794/1‑150 (MQ=255)
                                                                                                                                 gAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGa                     <  1:4132938/150‑1 (MQ=255)
                                                                                                                                 gAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGa                     <  1:4132942/150‑1 (MQ=255)
                                                                                                                                   aaaaCGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATa                   <  1:749971/150‑1 (MQ=255)
                                                                                                                                      aCGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATc                <  1:4492547/150‑1 (MQ=255)
                                                                                                                                        ggCCTCGCCTTTCGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGc                                          <  1:2571679/122‑1 (MQ=255)
                                                                                                                                        ggCCTCGCCTTTCGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGc                                          >  2:2571679/1‑122 (MQ=255)
                                                                                                                                        ggCCTCGCCTTTCGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGa                                        <  2:1308665/124‑1 (MQ=255)
                                                                                                                                        ggCCTCGCCTTTCGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGa                                        >  1:1308665/1‑124 (MQ=255)
                                                                                                                                        ggCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAaagaa                                                                                             >  1:2294931/1‑71 (MQ=255)
                                                                                                                                        ggCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAaagaa                                                                                             <  2:2294931/71‑1 (MQ=255)
                                                                                                                                          ccTCGCCGAACGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCAc            <  2:4959490/150‑1 (MQ=255)
                                                                                                                                            tCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGa          >  1:3384020/1‑150 (MQ=255)
                                                                                                                                              gCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGacc        >  2:3630662/1‑148 (MQ=255)
                                                                                                                                              gCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAAc        >  2:3267466/1‑150 (MQ=255)
                                                                                                                                              gCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAAc        >  2:4712744/1‑150 (MQ=255)
                                                                                                                                              gCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAAc        <  1:5423320/150‑1 (MQ=255)
                                                                                                                                                 gAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAACCgg     >  1:807725/1‑150 (MQ=255)
                                                                                                                                                 gAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAACCgg     >  1:3343573/1‑150 (MQ=255)
                                                                                                                                                 gAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAACCgg     <  2:3936975/150‑1 (MQ=255)
                                                                                                                                                 gAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAACCgg     >  2:3068593/1‑150 (MQ=255)
                                                                                                                                                 gAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAACCgg     <  1:4931120/150‑1 (MQ=255)
                                                                                                                                                  aaaGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTGCGATTATCGAACAAAcggc     >  2:2037469/1‑146 (MQ=255)
                                                                                                                                                  aaaGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTGCGATTATCGAACAAAcggc     <  1:2037469/149‑4 (MQ=255)
                                                                                                                                                    aGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAACCGgtgt  <  2:2928794/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AATTTACGCAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAACCGGTGT  >  NC_002947/1566048‑1566345

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: