Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,583,420 T→A 12.1% intergenic (+105/‑219) PP_3158 → / → PP_5589 hypothetical protein/rhs repeat‑associated core domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,583,4200TA12.1% 810.6 / 48.7 322intergenic (+105/‑219)PP_3158/PP_5589hypothetical protein/rhs repeat‑associated core domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (130/144);  new base A (36/9);  total (169/153)
Fisher's exact test for biased strand distribution p-value = 4.28e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 340 total aligned reads displayed.

ACGCGCAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGT  >  NC_002947/3583272‑3583569
                                                                                                                                                     |                                                                                                                                                     
   cgcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAttttt                                                                                                                                                   <  1:5467181/150‑1 (MQ=255)
    gcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTg                                                                                                                                                  <  2:3553334/150‑1 (MQ=255)
    gcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTg                                                                                                                                                  >  2:8756062/1‑150 (MQ=255)
     ctgagCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGACATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGc                                                                                                                                                 <  1:500421/148‑1 (MQ=255)
     cAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGc                                                                                                                                                 <  2:6291178/150‑1 (MQ=255)
       gagCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCa                                                                                                                                               <  1:8460872/150‑1 (MQ=255)
          ccTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGtt                                                                                                                                            <  2:7000934/150‑1 (MQ=255)
          ccTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGtt                                                                                                                                            <  2:8257457/150‑1 (MQ=255)
          ccTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAAAAATGCCAGtt                                                                                                                                            <  2:218259/150‑1 (MQ=255)
          ccTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCGGCCACAGGCCttttttacaggt                                                                                                                                            >  1:904708/1‑143 (MQ=255)
           cTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGtt                                                                                                                                            <  1:5707861/149‑1 (MQ=255)
           cTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGtt                                                                                                                                            >  2:5707861/1‑149 (MQ=255)
             gTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAAAAATGCCAGTTGCa                                                                                                                                         <  1:229472/150‑1 (MQ=255)
             gTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAAAAATGCCAGTTGCa                                                                                                                                         <  2:4802354/150‑1 (MQ=255)
              ttGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAATGGCATTTTTGCCAGTTGCAg                                                                                                                                        >  1:8296525/1‑150 (MQ=255)
               tGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGc                                                                                                                                       <  1:3032713/150‑1 (MQ=255)
                 cTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGa                                                                                                                                     <  1:7971296/150‑1 (MQ=255)
                 cTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATATTTGCCAGTTGCAGCGa                                                                                                                                     >  1:594119/1‑150 (MQ=255)
                    atgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAgtg                                                                                                                                  >  2:1446641/1‑150 (MQ=255)
                      gaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTTCAAATGGCAATTTTGACAGTTGCAGCGAGTGTa                                                                                                                                >  2:3329469/1‑150 (MQ=255)
                       aaCTTCCTGCCGGTGCTGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCAATACATGATTGAGTATTGAGTTTATGTGTGTTCATGAGCGAATCAAGTTTTCGTGGGGCTGAAACTGGATTTTTTGCCTGTTGCAGCGAGTGTAt                                                                                                                               <  1:6174635/150‑1 (MQ=255)
                        tcTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTGCAGTTGCAGCGAGTGTTTg                                                                                                                              >  1:2208087/2‑150 (MQ=255)
                         cTTCCTGCCGGTGCAGTAAGAGGCCTCGGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGTTTTGAGTTGTTGTGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGATTTTTTGCCTGTTGCAGCGAGTGTATGa                                                                                                                             <  1:4841234/150‑1 (MQ=255)
                                 cGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTTCCACTTGCAGCGAGTGTATGACCGCCGAg                                                                                                                     >  2:1867887/1‑150 (MQ=255)
                                 cGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAATGGCATTTTTGCCAGTTGAAGCGAGTGTATTACAGCCTAg                                                                                                                     >  2:5182330/1‑150 (MQ=255)
                                  ggTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGa                                                                                                                    <  1:5849060/150‑1 (MQ=255)
                                    tGCAGTAAGAGGCCTCGGCCGCACGCCCTTTGCGAAATGCCCCTACATGAATGCGGATTGAGTTGAGGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCGAGAAt                                                                                                                  <  2:6901355/150‑1 (MQ=255)
                                        ttAAGAGGCCTCTGCCGCACGCCCTTTGCTAAATGCCCATACATTATTGAGGTTTGATTTGATGTGTTTTCATGAGCTAATAAACGTTTCGTCGGGCTGAAACTGGCATTTTTGCCTGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              <  1:2493055/149‑1 (MQ=255)
                                            gagGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTTCAGCGAGTGTATGAAAGCCTAGAAATCAAGGGa                                                                                                          >  2:7464393/1‑150 (MQ=255)
                                            gagGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAATGGCATTTTTTCCAGTTGCAGCGAGTGTATGACAGCCAAGAAATAAAGGGa                                                                                                          >  2:4195798/1‑150 (MQ=255)
                                             agGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTTCCAGTTGCAGCGAGTGTATTACAGCCGAGAATTCAAGGGaa                                                                                                         >  1:633299/1‑150 (MQ=255)
                                                ccTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTGCAGTTGCAGCCAGTTTTTTACAGCCGAGAAATCAATGGaaaag                                                                                                      >  1:5145675/1‑148 (MQ=255)
                                                  tctGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTTCAGCGAGTGTATGACAGCCTAGAATTAAGGGGAAATGgt                                                                                                    >  2:5852454/1‑150 (MQ=255)
                                                   ctGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGAGTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCACCACGCAATTTTGCCAGTTTCAGCGCGTGTATTAAAGCCGAGAAATAAACGGAAAAGgtg                                                                                                   >  1:1922394/1‑150 (MQ=255)
                                                         cACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCAGCAAATGGCAATTTTTCCACTTGCAGCGAGTGTATGAAAGCCTAGAAATCAATGGAAAAGGTGTTatct                                                                                             >  2:5783803/1‑150 (MQ=255)
                                                          aCGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCCTTTTTTGCAGTTTCAGCGCGTGTATTAAAGCCGAGAAATAAATGGAAATGGTGTTgtctc                                                                                            >  1:168202/1‑149 (MQ=255)
                                                           cGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTaa                                                                                           <  2:4303187/150‑1 (MQ=255)
                                                            gCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCTAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAg                                                                                          <  2:1821704/150‑1 (MQ=255)
                                                            gCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTTCCAGTTGCAGCGAGTGTATGACAGCCGAGAATTCAAGGGCAATGGTGTTCTCTAAg                                                                                          >  2:6824904/1‑150 (MQ=255)
                                                            gCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAg                                                                                          >  2:5061843/1‑150 (MQ=255)
                                                            gCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAg                                                                                          <  2:787463/150‑1 (MQ=255)
                                                              ccTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTTCAACAGGCAATTTTGCCAGTTGCAGCGAGTGTATTACAGCCAAGAATTAAAGGGAAATTGTGTTATCTAAGGt                                                                                        >  1:1304442/1‑150 (MQ=255)
                                                              ccTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAAGGGCATGTTTTCCACTTGCAGCGCGTGTATTACAGCCTAGAAAAAAAGGGAAAAGGTGTTCTCTTAGGt                                                                                        >  1:4938970/1‑150 (MQ=255)
                                                                tttGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAATTTTTGCACGTTCCGCGAGTGTATGTAAGCCGAGAAATGAAGGGGAATGGTGGTGTCTGAGGTCa                                                                                      >  1:1097508/1‑150 (MQ=255)
                                                                    cGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAATTTTTCCACTTTCAGCGAGTGTATGACAGCCGAGAATTAAAGGGAAAATGTGTTGTCTAAGGTCATTaa                                                                                  >  2:382081/1‑150 (MQ=255)
                                                                        aggCCCATACTTGATTGAGTATTTTGTTGATGAGTGTTCATGAGCGAAACAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGt                                                                              <  1:4682338/148‑1 (MQ=255)
                                                                        aTGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCTAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGt                                                                              <  2:4091474/150‑1 (MQ=255)
                                                                        aTGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGt                                                                              <  1:1998411/150‑1 (MQ=255)
                                                                          ccccATACATCAATGAGGATTGAGTTGATCAGTGTTCATGACCGAAAAAACGTTTCGTCGGGCTGCAACTCGCATTTTTGCCAGTTCCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGa                                                                            <  1:4183194/149‑1 (MQ=255)
                                                                            ccATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAATTTTGCCAGTTTCAGCGAGTGTATTACAGCCAAGAATTAAAGGGAAATGGTGTTCTCTAAGGTCATTAATAGTGaaa                                                                          >  1:3520120/1‑150 (MQ=255)
                                                                             cATACATGTATGAGGATTGAGTTGAGGAGTGTTCATGAGCGAATCAAGGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAAc                                                                         <  2:8603167/150‑1 (MQ=255)
                                                                             cATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAAc                                                                         <  1:3323384/150‑1 (MQ=255)
                                                                              aTACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAAAGATTCGTCGGGCTGAAACAGGCAATTTTTCCAGATGCAACCAGATTCATACAGCCTACAATTAAAAGGAAATGGTGTTCTATAACATCAATAATCTTGaaaaa                                                                        >  2:5150948/1‑148 (MQ=255)
                                                                                aCATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGtt                                                                      <  1:8729520/150‑1 (MQ=255)
                                                                                 cATGACTGCGGTTTGTGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGGCGGGCTGCACCTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTc                                                                     <  2:1668252/150‑1 (MQ=255)
                                                                                 cATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTc                                                                     <  2:8030733/150‑1 (MQ=255)
                                                                                    gaatgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTATAAATCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCa                                                                  >  1:7300178/1‑150 (MQ=255)
                                                                                      atgtGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAg                                                                <  1:7776157/150‑1 (MQ=255)
                                                                                      atgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGAAACGGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAg                                                                <  1:1816833/150‑1 (MQ=255)
                                                                                        gaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCg                                                              <  2:1970657/150‑1 (MQ=255)
                                                                                         aGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGt                                                             <  2:1056707/150‑1 (MQ=255)
                                                                                          ggATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAATTTTGCCACTTTCAGCGAGTGTATGACAGCCTAGAATTAAAGGGAAATGGTGTTATCTAAGGTCATTAATAGTGAAACGTTAGCAAGCGTa                                                            >  2:3765571/1‑150 (MQ=255)
                                                                                           gATTTATTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTac                                                           <  1:7140940/150‑1 (MQ=255)
                                                                                                    gatgcGTGTTCATGCGCTAATAAACGTTGCGTCGGGCTGCAACTGGCATTTTTGCCCGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTAc                                                  <  1:3670219/150‑1 (MQ=255)
                                                                                                    gatgaGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAAAGGCATTTTTTCCAGTTGCAGCGCGTGTATGACAGCCAAGAAATCAAGGGAAAAGGTGTTCTATAAGGTCAATAAACGTGAAAAGTGCGCAAGCGGACACGCCCGAc                                                  >  2:8529200/1‑150 (MQ=255)
                                                                                                      tgaGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACgg                                                <  2:303787/150‑1 (MQ=255)
                                                                                                      tgaGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAATGGCATTTTTTCCAGTTGCAGCGAGTGTATGACAGCCGAGAATTCAAGGGAAATGGTGTTCTATAAGGTCATTAATAGTGAAAAGTTCGCAAGCGTACACGCCCTACgg                                                >  2:6483265/1‑150 (MQ=255)
                                                                                                              cATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTAttt                                        <  2:615946/150‑1 (MQ=255)
                                                                                                                tGAGCTAATCAACTTTTCGTCGGGCTGCAACTGGAATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTgg                                      <  1:1170894/150‑1 (MQ=255)
                                                                                                                 gAGCTAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTggg                                     <  2:3413097/150‑1 (MQ=255)
                                                                                                                    cGAATCAACGTTTCGGCGGGCTGCAACTGGAATTTTTGCCTGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTg                                  <  1:6660663/150‑1 (MQ=255)
                                                                                                                       aTCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATg                               <  2:4271465/150‑1 (MQ=255)
                                                                                                                        tCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATgg                              <  1:2081301/150‑1 (MQ=255)
                                                                                                                         cAACGTTTCGTCGGGCTGCAACATGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCAAGAATTAAAGGGAAATGGTGTTCTCTAATGTCATTAATCGTGAAAAGTTGGCAAGCGTACATGCCCTTCCGCGCTATTTGGGGtggtggg                             >  1:4320965/1‑149 (MQ=255)
                                                                                                                              tttCGTCGGGCTGCAACAGGCATTTTTGCCAGTTTCAGCGAGTGTATGACAGCCTAGAAATAAAGGGAAAAGGTGTTCACTAAGGTCCATAATAGTGAAAACTTAGCAAGCGTAAACCCCCTAACGCGCGATTTTGGGTGATGGCCAATg                        >  2:2083910/1‑150 (MQ=255)
                                                                                                                                  gtcgGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTcc                    <  1:4096541/150‑1 (MQ=255)
                                                                                                                                    cgGGCTTCAACTGGCATTTTTTCCACTTGCAGCGATTGTATTACACCCTAGAATTCAAGGGAAATGGTTTTATCTAAGGTCAATAATCGTTAAAAGTTAGCAAGCGTACACGCCCTACGGCGCTATATGGGGTGATGGCCATAGCACCgt                  >  2:2245032/1‑150 (MQ=255)
                                                                                                                                     gggCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCgtg                 <  2:8118033/150‑1 (MQ=255)
                                                                                                                                       gCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGtt               <  1:3442380/150‑1 (MQ=255)
                                                                                                                                         tGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGa             <  2:904708/150‑1 (MQ=255)
                                                                                                                                         tGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGa             >  2:2373428/1‑150 (MQ=255)
                                                                                                                                         tGCAACTGGCATTTTTGCCAGTTGAAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGa             <  1:7664611/150‑1 (MQ=255)
                                                                                                                                         tGCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTgt             >  2:5419262/1‑149 (MQ=255)
                                                                                                                                         tGCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGa             >  1:7581746/1‑150 (MQ=255)
                                                                                                                                         tGCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGa             >  1:7105298/1‑150 (MQ=255)
                                                                                                                                         tGCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGa             >  2:4917427/1‑150 (MQ=255)
                                                                                                                                           cAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAg           >  1:3773672/1‑150 (MQ=255)
                                                                                                                                           cAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAg           >  2:2762241/1‑150 (MQ=255)
                                                                                                                                           cAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAg           >  2:8476659/1‑150 (MQ=255)
                                                                                                                                            aaCTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTAttt                                        <  1:2735650/120‑1 (MQ=255)
                                                                                                                                             aCTGGAATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTa         <  1:4501106/150‑1 (MQ=255)
                                                                                                                                               tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACa       >  1:5588320/1‑150 (MQ=255)
                                                                                                                                               tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACa       >  2:7762764/1‑150 (MQ=255)
                                                                                                                                               tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACa       >  2:2570579/1‑150 (MQ=255)
                                                                                                                                                ggCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACaa      <  2:1253047/150‑1 (MQ=255)
                                                                                                                                                 gCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAAt     >  1:5277148/1‑150 (MQ=255)
                                                                                                                                                 gCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAAt     >  2:2925690/1‑150 (MQ=255)
                                                                                                                                                    aTTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATgg   >  2:4542660/1‑150 (MQ=255)
                                                                                                                                                     tttttGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTTCGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGt  <  2:2208087/150‑1 (MQ=255)
                                                                                                                                                     tttttGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGt  <  2:6780715/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGCGCAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGT  >  NC_002947/3583272‑3583569

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