Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,041,297 T→A 8.7% intergenic (+227/‑435) PP_5680 → / → PP_4447 Cro/CI family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,041,2970TA 8.7% 269.9 / 14.0 92intergenic (+227/‑435)PP_5680/PP_4447Cro/CI family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (42/42);  new base A (3/5);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

GACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGAAAAGTCATT  >  NC_002947/5041148‑5041444
                                                                                                                                                     |                                                                                                                                                   
gACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCt                                                                                                                                                     <  2:1146702/150‑1 (MQ=255)
   ggTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTAt                                                                                                                                                  <  2:370481/150‑1 (MQ=255)
       ttGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCg                                                                                                                                              <  2:495885/150‑1 (MQ=255)
         gATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCATAAGGCGGt                                                                                                                                            >  2:1588447/1‑150 (MQ=255)
         gATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCATAAGGCGGt                                                                                                                                            <  1:1588447/150‑1 (MQ=255)
          aTCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGtt                                                                                                                                           <  2:1445848/150‑1 (MQ=255)
           tCCACAGAGGGTTAGGGATCGCCTCAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGTGAAGAGCCAGAGTCCCCTAGCAAAACCCCGGTCCGCCTTAGGGCGGttt                                                                                                                                          <  2:1207544/150‑1 (MQ=255)
           tCCACAGAGGGTTAGGGATAGTCTTAAGGCGTATTTGTCTGATCCCAACACCTCTGCGATGAGGCCATGCCCTCTTCATCCGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCAGAAGGCGGtta                                                                                                                                          <  1:203346/150‑2 (MQ=255)
           tCCACAGAGGGTTAGGGATAGGCTCAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGCACATGACCTCTTCATCGGCGCGCCTTACGCCCGCCGAGAAGAGCCAGATACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGttt                                                                                                                                          <  2:1207399/150‑1 (MQ=255)
              acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtctc                                                                                                                                       >  1:59629/1‑150 (MQ=255)
              acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCATAAGGCGGttaggg                                                                                                                                       <  2:860797/150‑5 (MQ=255)
                 gagGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGc                                                                                                                                    <  2:59629/150‑1 (MQ=255)
                        aGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAg                                                                                                                             >  1:248066/1‑150 (MQ=255)
                             tAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTctgc                                                                                                                        <  1:1780877/150‑1 (MQ=255)
                              aGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTctgct                                                                                                                       <  1:105385/150‑1 (MQ=255)
                                    aaGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAg                                                                                                                 <  1:107708/150‑1 (MQ=255)
                                    aaGGCGTAATTGTCTGATACCAACACATCTGCGATGACGACATGACCTCTTCATCGGCGCGCCTTCAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGCCCCCTAACCGCCTTATGGCTGTTTCTCTGCCTCCCAGTCTGCTGCGAAg                                                                                                                 <  1:1251614/150‑1 (MQ=255)
                                      ggCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGc                                                                                                               >  2:709353/1‑150 (MQ=255)
                                        cGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCg                                                                                                             <  2:1983600/150‑1 (MQ=255)
                                           aaTTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTc                                                                                                                           <  1:1925466/133‑1 (MQ=255)
                                           aaTTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTc                                                                                                                           >  2:1925466/1‑133 (MQ=255)
                                                acTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGTCGGTTTCTCTGCCTCCCAGTCTGCTGCGCCGGCCGATAGCCGa                                                                                                     >  2:970283/2‑150 (MQ=255)
                                                      aCCCACGCCTCGGCGGGGAGGACATGACCTCTTCATCGGCGCGCCGTAGGCCCGCCGGGAAGGGCCAGAGACCCCGAGCAAGGCCCCTAACCGCCTTATGGCGGTTTCTCTGCCCCCCAGGCTGCTGCGAAGGCCGATAGCCGAGTccaa                                                                                               <  2:912199/150‑1 (MQ=255)
                                                        cAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTccaacc                                                                                             >  1:656209/1‑150 (MQ=255)
                                                        cAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTccaacc                                                                                             >  2:1885690/1‑150 (MQ=255)
                                                          acacCCCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTACGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTccaaccaa                                                                                           <  2:2447169/150‑1 (MQ=255)
                                                            acCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCATCCGATGCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAAtt                                                                                         >  1:1981938/1‑150 (MQ=255)
                                                                      tGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGc                                                                                                               >  2:583138/1‑118 (MQ=255)
                                                                      tGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGc                                                                                                               <  1:583138/118‑1 (MQ=255)
                                                                       gAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGa                                                                              >  1:3784/1‑150 (MQ=255)
                                                                        aGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGAc                                                                             <  2:766438/150‑1 (MQ=255)
                                                                         ggACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTccaac                                                                                              >  2:1727276/1‑132 (MQ=255)
                                                                         ggACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTccaac                                                                                              <  1:1727276/132‑1 (MQ=255)
                                                                          gACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCt                                                                           >  2:2243579/1‑150 (MQ=255)
                                                                           aCATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTg                                                                          <  2:656209/150‑1 (MQ=255)
                                                                               gACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCt                                                                      >  1:1620804/1‑150 (MQ=255)
                                                                                 ccTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATgg                                                                                                                                                >  2:392771/1‑74 (MQ=255)
                                                                                 ccTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATgg                                                                                                                                                <  1:392771/74‑1 (MQ=255)
                                                                                 ccTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGa                                                                    <  2:1620804/150‑1 (MQ=255)
                                                                                  ctctTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCCCTTCTCTTTATAAACGTGTCCTGAc                                                                   >  2:1319039/1‑150 (MQ=255)
                                                                                  ctctTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGAc                                                                   >  2:2397376/1‑150 (MQ=255)
                                                                                   tctTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTATGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGAcc                                                                  >  1:2109848/1‑150 (MQ=255)
                                                                                     ttCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGc                                                                >  2:1786274/1‑150 (MQ=255)
                                                                                      tCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCa                                                               <  1:2243579/150‑1 (MQ=255)
                                                                                        aTCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCa                                                                                                                              <  1:1084074/85‑1 (MQ=255)
                                                                                        aTCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCa                                                                                                                              >  2:1084074/1‑85 (MQ=255)
                                                                                        aTCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGt                                                             >  2:2171832/1‑150 (MQ=255)
                                                                                        aTCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGt                                                             >  1:112096/1‑150 (MQ=255)
                                                                                        aTCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGt                                                             >  1:1479514/1‑150 (MQ=255)
                                                                                        aTCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGt                                                             >  1:1106511/1‑150 (MQ=255)
                                                                                             cgcgcCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGc                                                        <  2:1703988/150‑1 (MQ=255)
                                                                                               cgcCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCgt                                                      <  2:890032/150‑1 (MQ=255)
                                                                                               cgcCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCgt                                                      <  1:1319039/150‑1 (MQ=255)
                                                                                                   ttAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCt                                                                           >  2:2383885/1‑125 (MQ=255)
                                                                                                   ttAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCt                                                                           <  1:2383885/125‑1 (MQ=255)
                                                                                                     aaGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTa                                                <  1:578237/150‑1 (MQ=255)
                                                                                                       gCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAAt                                                                                          <  1:575037/106‑1 (MQ=255)
                                                                                                       gCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAAt                                                                                          >  2:575037/1‑106 (MQ=255)
                                                                                                       gCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATc                                              >  1:1159507/1‑150 (MQ=255)
                                                                                                        cccGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCa                                             <  1:950867/150‑1 (MQ=255)
                                                                                                        cccGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTCATGGCGGCGTCGCAGCGTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCa                                             >  2:2039795/1‑150 (MQ=255)
                                                                                                           gccgAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCagag                                          >  1:687325/1‑150 (MQ=255)
                                                                                                              gagaAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGaaa                                       <  1:2237939/150‑1 (MQ=255)
                                                                                                              gagaAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGaaa                                       <  1:994619/150‑1 (MQ=255)
                                                                                                               agaagaGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGaaaa                                      <  1:1786274/150‑1 (MQ=255)
                                                                                                               agaagaGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGaaaa                                      <  2:248066/150‑1 (MQ=255)
                                                                                                                 aagaGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGt                                    <  2:687325/150‑1 (MQ=255)
                                                                                                                  accgCCAGAGACCCCGAGCAAGACCCCTACCCGCCTTCCAACGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTg                                   <  1:22036/147‑1 (MQ=255)
                                                                                                                       cAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATc                              >  1:1788590/1‑150 (MQ=255)
                                                                                                                       cAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATc                              >  2:240883/1‑150 (MQ=255)
                                                                                                                            aCCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAatat                         <  1:1885690/150‑1 (MQ=255)
                                                                                                                              cccGAGCAAGACCCCTCACCGCCTTATGGCGGTTTTTCTGCCTCCCAGTCTGGTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGCCCTGGCCTGACCGCAGTAGAGCTTGTCCTATCAGAGAAAAGTGGTATCAACATcc                       >  2:144554/1‑149 (MQ=255)
                                                                                                                              cccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCa                                                               >  2:596649/1‑110 (MQ=255)
                                                                                                                              cccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCa                                                               <  1:596649/110‑1 (MQ=255)
                                                                                                                              cccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCa                       >  1:9514/1‑150 (MQ=255)
                                                                                                                               ccGAGCAAGCCCCCTACCCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATGCTATTTTTCGACCTGGCCTGATCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAg                      <  2:880633/150‑1 (MQ=255)
                                                                                                                               ccGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAg                      >  2:1952329/1‑150 (MQ=255)
                                                                                                                                          gagaaaCCGCCATAAGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAATTGGTATCAATATCAGCCCCCCACCGa           >  1:2342002/5‑150 (MQ=255)
                                                                                                                                          cccTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCCGAAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGa           >  1:1730319/1‑150 (MQ=255)
                                                                                                                                          cccTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGCGAAAAGTGGTATCAATATCAGCCCCCCACCGa           >  2:149183/1‑150 (MQ=255)
                                                                                                                                          cccTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGa           >  2:712571/1‑150 (MQ=255)
                                                                                                                                          cccTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTAGCCGCTGGCCTGACCGCAGTAGAGCGTGTGCTATGCGAGAAGGGAGCTATCTGCGTCAGCCCCCCACCGa           >  1:42919/1‑150 (MQ=255)
                                                                                                                                           ccTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGaa          >  1:858356/1‑150 (MQ=255)
                                                                                                                                           agaaaCCGCCATAAGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGc                                                  >  1:2154243/4‑110 (MQ=255)
                                                                                                                                           agaaaCCGCCATAAGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGc                                                  <  2:2154243/107‑1 (MQ=255)
                                                                                                                                             aaaCCGCCATAAGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGaaaa        <  2:2342002/149‑1 (MQ=255)
                                                                                                                                              aaCCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGc                     >  2:1899634/1‑136 (MQ=255)
                                                                                                                                              aaCCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGc                     <  1:1899634/136‑1 (MQ=255)
                                                                                                                                                  gCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCa                                                               >  1:433507/1‑90 (MQ=255)
                                                                                                                                                  gCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCa                                                               <  2:433507/90‑1 (MQ=255)
                                                                                                                                                   ccTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGc                                                        <  2:1661992/96‑1 (MQ=255)
                                                                                                                                                   ccTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGc                                                        >  1:1661992/1‑96 (MQ=255)
                                                                                                                                                   ccTTATGGCGGTTTCTATGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGAAAAGTCAtt  <  2:2109848/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTGCTATCAGAGAAAAGTGGTATCAATATCAGCCCCCCACCGAAAAGTCATT  >  NC_002947/5041148‑5041444

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