Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,741,977 C→G 7.4% intergenic (‑23/+61) 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,9770CG 7.4% 290.5 / 13.1 95intergenic (‑23/+61)PP_1538/PP_1539hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (46/42);  new base G (3/4);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.78e-01

GTGACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCGAACCTTCCG  >  NC_002947/1741836‑1742119
                                                                                                                                             |                                                                                                                                              
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:1518302/1‑150 (MQ=255)
     tGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATa                                                                                                                                   >  1:371914/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                                                                                                                         <  2:225293/150‑1 (MQ=255)
               aaGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGagag                                                                                                                         <  1:997728/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                                            >  1:655665/1‑150 (MQ=255)
                                                                                            gATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGcc                                            >  2:1739773/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:1243048/1‑126 (MQ=255)
                                                                                                    gTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTgc                                                            >  1:1243031/1‑126 (MQ=255)
                                                                                                    gTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTgc                                                            <  2:1243031/126‑1 (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                                                 >  2:1381763/1‑127 (MQ=255)
                                                                                                              ttGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGt                                                 <  1:1381763/127‑1 (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)
                                                                                                                                            tCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGc                                                      >  1:959093/1‑92 (MQ=255)
                                                                                                                                            tCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGc                                                      <  2:959093/92‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                              
GTGACTGAAGAGGGAAAGCCTGGCTTAACCGGGCGTCGGGTGACTGTCTTCCGTCCTTCCAGGATGGCTCGGACCATCGGCCCGCTGAACAGGATCGGGCGTTCCTTTGCTTGAGACATAGCTTCGCCTTTGCCGCCATATCGCGGCAGTGAATAGAGTGGAGAGGGGTTACAGCGGGGTGGGGTACAGATGTGCTCTGGGGTTAAGAGCTGGCACCGGGGCGTGCGCAAGCGTTGTCAGCCATGGCGCGGCACTCTACCTGGTGCGGTCCGTCGAACCTTCCG  >  NC_002947/1741836‑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: