Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,566,194 A→T 5.6% intergenic (‑311/‑169) 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,1940AT 5.6% 399.6 / 8.3 118intergenic (‑311/‑169)PP_1374/pcaRhypothetical protein/pca regulon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (51/60);  new base T (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 124 total aligned reads displayed.

GATAATTTACGCAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAACCGG  >  NC_002947/1566045‑1566342
                                                                                                                                                     |                                                                                                                                                    
gATAATTTACGCAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAATACGGCCTCGCCGa                                                                                                                                                      >  2:4132942/1‑150 (MQ=255)
gATAATTTACGCAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGa                                                                                                                                                      >  1:498827/1‑150 (MQ=255)
gATAATTTACGCAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGa                                                                                                                                                      >  2:689259/1‑150 (MQ=255)
   aaTTTACGCAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAg                                                                                                                                                   <  1:1969899/150‑1 (MQ=255)
     tttACGCAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGc                                                                                                                                                 <  2:301325/150‑1 (MQ=255)
         cGCAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAgg                                                                                                                                             <  2:5402859/150‑1 (MQ=255)
                cTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGc                                                                                                                                      <  1:1779131/150‑1 (MQ=255)
                   aCATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTg                                                                                                                                   <  2:2816823/150‑1 (MQ=255)
                     atcatcTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTggg                                                                                                                                 >  1:2197831/1‑150 (MQ=255)
                     atcatcTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTggg                                                                                                                                 >  1:2420171/1‑150 (MQ=255)
                     atcatcTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTggg                                                                                                                                 >  2:946773/1‑150 (MQ=255)
                        atcTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGAc                                                                                                                              <  1:1714846/150‑1 (MQ=255)
                        atcTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGAc                                                                                                                              >  2:4669351/1‑150 (MQ=255)
                        atcTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGAc                                                                                                                              >  2:4557504/1‑150 (MQ=255)
                          cTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACtt                                                                                                                            >  2:5423320/1‑150 (MQ=255)
                          cTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACtt                                                                                                                            >  2:4794740/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)
                               gTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACg                                                                                                                       >  2:1671288/1‑150 (MQ=255)
                                tAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACgg                                                                                                                      <  1:2150236/150‑1 (MQ=255)
                                 aaGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGt                                                                                                                     <  1:976429/150‑1 (MQ=255)
                                  aGTAGCAGATTGTGCCCGCGAACATTTGGGAACTTTCGTACAAAGTTTGAAAGGAATGGCGGATACAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACtgtg                                                                                                                    <  2:2058437/150‑1 (MQ=255)
                                    tAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTggg                                                                                                                  >  2:2234352/1‑150 (MQ=255)
                                    tAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTggg                                                                                                                  >  2:2234237/1‑150 (MQ=255)
                                      gCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTa                                                                                                                >  2:4931120/1‑150 (MQ=255)
                                      gCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTa                                                                                                                >  1:3936975/1‑150 (MQ=255)
                                       cAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTAc                                                                                                               <  1:1329228/150‑1 (MQ=255)
                                       cAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTAc                                                                                                               >  1:927360/1‑150 (MQ=255)
                                          aTTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACg                                                                                                                       >  1:1605059/1‑139 (MQ=255)
                                          aTTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACg                                                                                                                       <  2:1605059/139‑1 (MQ=255)
                                               gCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAg                                                                                                       >  2:2506240/1‑150 (MQ=255)
                                                cccGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGt                                                                                                      <  1:3421458/150‑1 (MQ=255)
                                                 ccGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGtc                                                                                                     >  2:1099485/1‑150 (MQ=255)
                                                 ccGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGtc                                                                                                     >  1:2802528/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)
                                                           ttGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAaaga                                                                                           <  1:4082590/150‑1 (MQ=255)
                                                            tGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAaagaa                                                                                          <  2:2023506/150‑1 (MQ=255)
                                                              ggACCTTTCGTACACAGGTTGCCAGGACTGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCACAGAAGc                                                                                        <  2:1298438/150‑1 (MQ=255)
                                                              ggACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGc                                                                                        <  2:1298437/150‑1 (MQ=255)
                                                                   tttCGTACAAAGTTGGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACaa                                                                                   <  1:4686590/150‑1 (MQ=255)
                                                                    ttCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAg                                                                                  >  2:1916361/1‑150 (MQ=255)
                                                                     tCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGt                                                                                 >  1:2428775/1‑150 (MQ=255)
                                                                       gTACAAAGTTTGAAAGGACTGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTATGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTgg                                                                               <  2:3371827/150‑1 (MQ=255)
                                                                       gTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTgg                                                                               >  2:3248357/1‑150 (MQ=255)
                                                                       gTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTgg                                                                               >  1:4355111/1‑150 (MQ=255)
                                                                        tACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGc                                                                              >  1:484816/1‑150 (MQ=255)
                                                                         aCAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCa                                                                             <  1:1916361/150‑1 (MQ=255)
                                                                           aaaGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCAt                                                                            >  2:2478792/1‑149 (MQ=255)
                                                                           aaaGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATg                                                                           >  2:1010400/1‑150 (MQ=255)
                                                                           aaaGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATg                                                                           >  2:1801393/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:4277118/150‑1 (MQ=255)
                                                                                 tGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAg                                                                     <  1:4146049/150‑1 (MQ=255)
                                                                                   aaaGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGt                                                                   >  2:909627/1‑150 (MQ=255)
                                                                                   aaaGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGt                                                                   >  2:909557/1‑150 (MQ=255)
                                                                                     aGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGtat                                                                 <  2:4414211/150‑1 (MQ=255)
                                                                                       gAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATaa                                                               >  2:4720314/1‑150 (MQ=255)
                                                                                               gATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATaa                                                               >  2:1939668/1‑142 (MQ=255)
                                                                                                     aGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTcggc                                                 <  1:5048237/150‑1 (MQ=255)
                                                                                                        tCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGt                                              <  1:2504596/150‑1 (MQ=255)
                                                                                                          gTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCa                                            >  2:1620752/1‑150 (MQ=255)
                                                                                                          gTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCa                                            >  1:4078950/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)
                                                                                                                aTGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCg                                      <  1:2624907/150‑1 (MQ=255)
                                                                                                                 tGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGa                                     <  1:818347/150‑1 (MQ=255)
                                                                                                                   tCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCaaaa                                                         >  2:1291896/1‑128 (MQ=255)
                                                                                                                   tCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCaaaa                                                         <  1:1291896/128‑1 (MQ=255)
                                                                                                                   tCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGAtt                                   >  1:3304801/1‑150 (MQ=255)
                                                                                                                      gggcgggcCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATc                                >  1:4247940/1‑150 (MQ=255)
                                                                                                                      gggcgggcCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATc                                >  1:2352230/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)
                                                                                                                                    gAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGa                  <  1:4132938/150‑1 (MQ=255)
                                                                                                                                      aaaaCGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATa                <  1:749971/150‑1 (MQ=255)
                                                                                                                                         aCGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATc             <  1:4492547/150‑1 (MQ=255)
                                                                                                                                           ggCCTCGCCTTTCGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGc                                       >  2:2571679/1‑122 (MQ=255)
                                                                                                                                           ggCCTCGCCTTTCGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGc                                       <  1:2571679/122‑1 (MQ=255)
                                                                                                                                           ggCCTCGCCTTTCGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGa                                     >  1:1308665/1‑124 (MQ=255)
                                                                                                                                           ggCCTCGCCTTTCGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGa                                     <  2:1308665/124‑1 (MQ=255)
                                                                                                                                           ggCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAaagaa                                                                                          >  1:2294931/1‑71 (MQ=255)
                                                                                                                                           ggCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAaagaa                                                                                          <  2:2294931/71‑1 (MQ=255)
                                                                                                                                             ccTCGCCGAACGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCAc         <  2:4959490/150‑1 (MQ=255)
                                                                                                                                                 gCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGacc     >  2:3630662/1‑148 (MQ=255)
                                                                                                                                                 gCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAAc     <  1:5423320/150‑1 (MQ=255)
                                                                                                                                                    gAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAACCgg  <  2:3936975/150‑1 (MQ=255)
                                                                                                                                                    gAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAACCgg  >  2:3068593/1‑150 (MQ=255)
                                                                                                                                                    gAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAACCgg  >  1:3343573/1‑150 (MQ=255)
                                                                                                                                                    gAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAACCgg  <  1:4931120/150‑1 (MQ=255)
                                                                                                                                                     aaaGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTGCGATTATCGAACAAAcggc  <  1:2037469/149‑4 (MQ=255)
                                                                                                                                                     aaaGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTGCGATTATCGAACAAAcggc  >  2:2037469/1‑146 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATAATTTACGCAACCCTGACATCATCTGGCGTAAGTAGCAGATTGTGCCCGCGAACAGTTGGGACCTTTCGTACAAAGTTTGAAAGGAATGGCGGATAGAAGTTCGTTGCAATGTCAGGGCGGGCCGTAATGAAAAACGGCCTCGCCGAAAGGCGAGGCCGTTGCCTGGGGACTTCTACGGTGGGTACCTGACCAGTCTCTCCAAAGAAGCCACAAGTGGCATGGGAAAGGTATAAGCAAAAGTGTCGGCGGTCAATTGCGATTATCGGCCGTTTGTTCGATAATCGCACGAACCGG  >  NC_002947/1566045‑1566342

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: