Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,741,972 C→G 7.8% intergenic (‑18/+66) PP_1538 ← / ← PP_1539 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,741,9720CG 7.8% 275.2 / 13.2 90intergenic (‑18/+66)PP_1538/PP_1539hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (44/39);  new base G (3/4);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.71e-01

GGAACACCTCGGTGACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCGAACCTTCCG  >  NC_002947/1741825‑1742119
                                                                                                                                                   |                                                                                                                                                   
ggAACACCTCGGTGACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCat                                                                                                                                                   >  2:63240/1‑150 (MQ=255)
           gTGACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCGATATGGCGGCaaa                                                                                                                                        <  2:1311864/150‑3 (MQ=255)
           gTGACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCGATATGGCGGCaaa                                                                                                                                        <  2:1692706/150‑3 (MQ=255)
           gTGACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGt                                                                                                                                        <  1:461508/150‑1 (MQ=255)
             gACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCCTATCGCGGCAGTGa                                                                                                                                      >  2:429578/1‑150 (MQ=255)
              aCTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGaa                                                                                                                                     <  1:1785367/150‑1 (MQ=255)
               cTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCACGGCAGTGAAt                                                                                                                                    >  1:1717725/1‑150 (MQ=255)
                tGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGGCATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATa                                                                                                                                   >  1:1231432/1‑150 (MQ=255)
                tGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATa                                                                                                                                   >  1:371914/1‑150 (MQ=255)
                tGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATa                                                                                                                                   >  1:1518302/1‑150 (MQ=255)
                   agagGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTGCGTCCTGCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTCCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATagag                                                                                                                                >  2:1657215/1‑150 (MQ=255)
                   agagGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATagag                                                                                                                                >  2:1251260/1‑150 (MQ=255)
                    gagGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGt                                                                                                                               <  2:837204/150‑1 (MQ=255)
                       ggAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGga                                                                                                                            <  2:1565738/150‑1 (MQ=255)
                         aaaGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGgagg                                                                                                                          >  2:930588/1‑149 (MQ=255)
                         aaaGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGgaga                                                                                                                          >  1:358673/1‑150 (MQ=255)
                          aaGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGagag                                                                                                                         <  1:997728/150‑1 (MQ=255)
                          aaGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGagag                                                                                                                         <  2:225293/150‑1 (MQ=255)
                            gCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCGTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAggg                                                                                                                       <  1:1251260/150‑1 (MQ=255)
                             ccTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAgggg                                                                                                                      <  1:1760431/150‑1 (MQ=255)
                              cTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGGCATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGt                                                                                                                     <  2:1231432/150‑1 (MQ=255)
                                 gCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTAc                                                                                                                  >  2:240817/1‑150 (MQ=255)
                                   ttAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGCGGAGAGGGGTTACAg                                                                                                                >  2:1520100/1‑150 (MQ=255)
                                     aaCCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCg                                                                                                              >  1:1584659/1‑150 (MQ=255)
                                          ggCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGGCACAGCggggtg                                                                                                         >  2:586778/1‑150 (MQ=255)
                                            cGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCggggtggg                                                                                                       <  1:163453/150‑1 (MQ=255)
                                            cGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCggggtggg                                                                                                       <  1:240817/150‑1 (MQ=255)
                                            cGTCGGGTGACGGTATCGCGTCCTTCTAGGATGGCTTGGCGCATCGGCGCGCGGGACAGGCGCGGGCGTTCCCGTTCTTGACACATAGCTTCGCCTTTTTCGCCACATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCgaggtggg                                                                                                       <  2:799410/150‑1 (MQ=255)
                                             gTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCggggtgggg                                                                                                      >  2:1302857/1‑150 (MQ=255)
                                             gTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCggggtgggg                                                                                                      <  2:1765390/150‑1 (MQ=255)
                                                 ggTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACa                                                                                                  <  1:1850595/150‑1 (MQ=255)
                                                   tGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTAc                                                                                                                  <  2:1300599/132‑1 (MQ=255)
                                                   tGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTAc                                                                                                                  >  1:1300599/1‑132 (MQ=255)
                                                   tGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGa                                                                                                >  1:650254/1‑150 (MQ=255)
                                                     aCTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGAtg                                                                                              >  1:777785/1‑150 (MQ=255)
                                                           ttCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGctct                                                                                        >  1:1047036/1‑150 (MQ=255)
                                                             ccGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACCGCCGGGGGGGGTACAGATGTGCTCTgg                                                                                      >  2:904626/1‑150 (MQ=255)
                                                              cGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTggg                                                                                     <  2:1584659/150‑1 (MQ=255)
                                                              cGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTggg                                                                                     <  2:650254/150‑1 (MQ=255)
                                                                   ttCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCGATATGGCGGCaaag                                                                                                                                       >  1:300174/1‑92 (MQ=255)
                                                                   ttCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCGATATGGCGGCaaag                                                                                                                                       <  2:300174/95‑4 (MQ=255)
                                                                   ttCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGAt                                                                                               <  1:695445/135‑1 (MQ=255)
                                                                   ttCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGAt                                                                                               >  2:695445/1‑135 (MQ=255)
                                                                   ttCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTaa                                                                                >  2:1596066/1‑150 (MQ=255)
                                                                     ccAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTAGAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAaga                                                                              <  1:1838683/150‑1 (MQ=255)
                                                                     ccAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAaga                                                                              >  1:1183853/1‑150 (MQ=255)
                                                                       aGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGc                                                                            >  1:11368/1‑150 (MQ=255)
                                                                          aTGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGCGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTgg                                                                         <  2:1215666/150‑1 (MQ=255)
                                                                          aTGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTgg                                                                         <  1:1302857/150‑1 (MQ=255)
                                                                          aTGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTgg                                                                         <  1:1657215/150‑1 (MQ=255)
                                                                               tCGGACCATCGGCCCGCTGAACGGGATCGGGCGTTCCTTTGATTGAGACATATCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGGGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGGACCg                                                                    >  2:1709145/1‑150 (MQ=38)
                                                                                      aTCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTg                                                             >  2:1383280/1‑150 (MQ=255)
                                                                                       tCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTgc                                                            >  1:760220/1‑150 (MQ=255)
                                                                                                       gATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGcc                                            >  2:1739773/1‑150 (MQ=255)
                                                                                                       gATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGcc                                            >  1:655665/1‑150 (MQ=255)
                                                                                                          cGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATg                                         <  2:835291/150‑1 (MQ=255)
                                                                                                           gggCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATgg                                        >  1:571876/1‑150 (MQ=255)
                                                                                                              cGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGGGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGgcgc                                     <  1:904626/150‑1 (MQ=255)
                                                                                                              cGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGgcgc                                     <  2:358673/150‑1 (MQ=255)
                                                                                                               gTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTgc                                                            <  2:1243048/126‑1 (MQ=255)
                                                                                                               gTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTgc                                                            >  1:1243031/1‑126 (MQ=255)
                                                                                                               gTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTgc                                                            <  2:1243031/126‑1 (MQ=255)
                                                                                                               gTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTgc                                                            >  1:1243048/1‑126 (MQ=255)
                                                                                                                 tCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAggg                                                                                                                       >  1:1713289/1‑65 (MQ=255)
                                                                                                                 tCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAggg                                                                                                                       <  2:1713289/65‑1 (MQ=255)
                                                                                                                  ccTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGGGGAGAGGGGTTACAGCGGGGGGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCa                                 >  2:1231989/1‑150 (MQ=255)
                                                                                                                   cTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCAc                                <  1:1596066/150‑1 (MQ=255)
                                                                                                                   cTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGCTGTGAGCCATGGCGCGGCAc                                >  1:186521/1‑150 (MQ=255)
                                                                                                                    tttGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCAct                               >  2:538177/1‑150 (MQ=255)
                                                                                                                         ttGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGt                                                 <  1:1381763/127‑1 (MQ=255)
                                                                                                                         ttGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGt                                                 >  2:1381763/1‑127 (MQ=255)
                                                                                                                         ttGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTAcc                          >  1:318938/1‑150 (MQ=255)
                                                                                                                         ttGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTAcc                          >  2:721877/1‑150 (MQ=255)
                                                                                                                            agaCATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTgg                       <  1:538177/150‑1 (MQ=255)
                                                                                                                                 tAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCgg                  >  1:367367/1‑150 (MQ=255)
                                                                                                                                  aGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGt                 <  2:1806785/150‑1 (MQ=255)
                                                                                                                                   gCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTc                >  1:392512/1‑150 (MQ=255)
                                                                                                                                      tCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGt             <  2:392512/150‑1 (MQ=255)
                                                                                                                                       cGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTc            <  1:63240/150‑1 (MQ=255)
                                                                                                                                       cGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTc            <  1:429578/150‑1 (MQ=255)
                                                                                                                                       cGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGAGCAATAACACTGAATTCGGGCGTTCAGAGCTGGCACCGGGGAGTGCGCAAACGTCGTCAGCCACTGCGCGGCACTCTACCTGGTGCGGTCCGTc            <  1:1383280/150‑1 (MQ=19)
                                                                                                                                        gCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCg           >  1:1916231/1‑150 (MQ=255)
                                                                                                                                         ccTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGt                                                 <  2:1564225/111‑1 (MQ=255)
                                                                                                                                         ccTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGGTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGt                                                 >  1:1564225/1‑111 (MQ=255)
                                                                                                                                          cactGCCGCGATATGGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCGaa         >  2:181493/4‑150 (MQ=255)
                                                                                                                                          cTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCGaa         <  2:318938/150‑1 (MQ=255)
                                                                                                                                          cTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCGaa         >  2:174917/1‑150 (MQ=255)
                                                                                                                                          cTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCGaa         <  2:1518302/150‑1 (MQ=255)
                                                                                                                                           actGCCGCGATATGGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGt                 <  1:128344/139‑1 (MQ=255)
                                                                                                                                           actGCCGCGATATGGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGt                 >  2:128344/3‑141 (MQ=255)
                                                                                                                                                cgccATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGTGTTAATAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAg                                              <  2:1163060/107‑1 (MQ=255)
                                                                                                                                                cgccATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAg                                              >  1:1163060/1‑107 (MQ=255)
                                                                                                                                                cgccATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCGAACCTTcc   <  2:367367/150‑1 (MQ=255)
                                                                                                                                                 gccATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTGGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCGAACCTTCCg  >  2:895348/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGAACACCTCGGTGACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCGAACCTTCCG  >  NC_002947/1741825‑1742119

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: