Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,580,393 T→A 11.7% 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,580,3930TA11.7% 488.4 / 17.1 191intergenic (+105/‑219)PP_3158/PP_5589hypothetical protein/rhs repeat‑associated core domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (73/87);  new base A (24/3);  total (101/90)
Fisher's exact test for biased strand distribution p-value = 1.87e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 204 total aligned reads displayed.

TACGCGCAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGT  >  NC_002947/3580244‑3580542
                                                                                                                                                     |                                                                                                                                                     
tACGCGCAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAt                                                                                                                                                       >  1:3973434/1‑150 (MQ=255)
    cgcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAttttt                                                                                                                                                   <  1:2038705/150‑1 (MQ=255)
    cgcAGAGCCTGGTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAttttt                                                                                                                                                   <  2:59746/150‑1 (MQ=255)
     gcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCCACAGGCAtttttt                                                                                                                                                  >  2:3488330/1‑149 (MQ=255)
     gcAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTg                                                                                                                                                  >  2:4047294/1‑150 (MQ=255)
      cAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGc                                                                                                                                                 <  2:3518307/150‑1 (MQ=255)
      cAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGc                                                                                                                                                 <  2:1811742/150‑1 (MQ=255)
      cAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGc                                                                                                                                                 <  1:2928894/150‑1 (MQ=255)
       agagCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGGTGAAACAGGCAATTTTGcc                                                                                                                                                >  1:3591484/1‑150 (MQ=255)
        gagCCTGTTGCTGATGAACTTCCTGCCGGTTCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCa                                                                                                                                               <  1:255107/150‑1 (MQ=255)
        gagCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCa                                                                                                                                               <  2:2016850/150‑1 (MQ=255)
           ccTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGtt                                                                                                                                            <  1:3785482/150‑1 (MQ=255)
            cTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAAAAATGCCAGTTg                                                                                                                                           <  1:344113/150‑1 (MQ=255)
              gTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCa                                                                                                                                         <  2:21739/150‑1 (MQ=255)
              gTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCa                                                                                                                                         <  1:1374117/150‑1 (MQ=255)
                   tgatgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCGGCAACAGGCCATTTTGCCACTTGCAGCGAg                                                                                                                                    >  1:3184013/1‑150 (MQ=255)
                   tgatgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCGGCAACAGGCAATTTTGCCACTTGCAGCGAg                                                                                                                                    >  1:212568/1‑150 (MQ=255)
                     atgaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAGTTTTTCCAGGTGCAGCGCgtg                                                                                                                                  >  1:2767788/1‑150 (MQ=255)
                       gaACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTGCAGTTGCAGCGAGTGTa                                                                                                                                >  1:2008190/1‑150 (MQ=255)
                        accTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGCGGATTGAGTTGATGAGTGTTCATGAGCGAGTCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTAt                                                                                                                               <  1:4462694/148‑1 (MQ=255)
                         aCTTCCTGCCGGTGCATTAAGAGGCCTCTGCCGCACGCCCTTTGCGAATTGCCCATAAATGAATGAGTTTTGTGTTGATGCGTGTTCATGAGCTAATCAACGTTTCGTCGGGCTGCAACTGGAATTTTTGCCTGTTGCAGCGAGTGTATg                                                                                                                              <  1:945663/150‑1 (MQ=255)
                          cTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATAGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACGTGCAATTTTTCCACTTGCACCGAGAGTAAGa                                                                                                                             >  2:34453/1‑150 (MQ=255)
                           ttCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAATGGCAATTTTTACACTTGCAGCGAGGGTATggc                                                                                                                            >  2:1174807/1‑148 (MQ=255)
                           ttCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCGGCAAAAGGCAATTTTTCCAGTTTCAGCGAGTGTATGAc                                                                                                                            >  1:109764/1‑150 (MQ=255)
                           ttCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGAGTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCCTTTTTGCCAGGTGCAGCGAGTGTATGAc                                                                                                                            >  2:28959/1‑150 (MQ=255)
                              cTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACCGc                                                                                                                         >  2:2724557/1‑150 (MQ=255)
                               tGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGGGGGCACAGGCCTTTTTTCCAGGTTCAGCGAGTTTAttaaaacc                                                                                                                        >  1:409123/1‑143 (MQ=255)
                               tGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGcc                                                                                                                        >  1:3628580/1‑150 (MQ=255)
                                 ccGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTTCAACTGGCCTTTTTGCCAGTTGCAGCGAGTGTATGACCGccga                                                                                                                      >  2:2419126/1‑148 (MQ=255)
                                 ccGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTa                                                                                                                      <  2:4325846/150‑1 (MQ=255)
                                 ccGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTa                                                                                                                      >  2:4344090/1‑150 (MQ=255)
                                    gTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCTAAAGGCCCATACATGAATGAGGATTGATTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGaa                                                                                                                   <  1:2534644/150‑1 (MQ=255)
                                    gTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGaa                                                                                                                   <  2:362820/150‑1 (MQ=255)
                                         gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGGCTTTTTTCCAGTTTCAGCGAGTGTATGAAAGCCAAGAATTaaa                                                                                                              >  1:615417/1‑150 (MQ=255)
                                         gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              <  2:1684847/150‑1 (MQ=255)
                                         gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              <  2:3229296/150‑1 (MQ=255)
                                         gTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACGGGAATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCaa                                                                                                              <  1:4006482/150‑1 (MQ=255)
                                          tAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCCTTTTTGGCAGTTTCAGCGCGTGGATGTCAGCCGAGAATTAAAg                                                                                                             >  1:3532022/1‑150 (MQ=255)
                                          tAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATTACAGCCAAGAATTCAAg                                                                                                             >  1:3506229/1‑150 (MQ=255)
                                           aaGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGCGGATTGAGTTGAGGAGTGTTCTTGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAgg                                                                                                            <  2:4173416/150‑1 (MQ=255)
                                                 ccTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGTTTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATg                                                                                                      <  1:2607375/150‑1 (MQ=255)
                                                  ctctGCCGCACGCCCTTTGCGAAATGCCCATACATGATTGAGTATTGTTTTTTTGCGTTTTAATGAGCTAATCAACGTTTCGTCGGGCTGAAACTGGAATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATgg                                                                                                     <  1:1943948/150‑1 (MQ=255)
                                                   tctGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGGGAGGGTTTTTCAGCCGAGGATTAAAGGGCAATGgt                                                                                                    >  1:4346729/1‑150 (MQ=255)
                                                    ctGCCGCACGCCCTTGGCGAAATGCCCATACATGAATGAGGATTGAGTTTATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGCAATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGgtg                                                                                                   <  2:3185925/150‑1 (MQ=255)
                                                     tGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGGTGCAACAGGGATTTTTTGCAGTTTCCGCGAGTTTATGACAGCCAAGAATTAAAGGGAAAATgtgt                                                                                                  >  1:1643315/1‑150 (MQ=255)
                                                          cACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTctct                                                                                             <  2:4000035/150‑1 (MQ=255)
                                                           aCGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGTGTTGATGTGTGTTCATGAGCGAATCAACTTTTCGTCGGGCTGCAACTGGCATTTTTGCCTGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTa                                                                                            <  1:632589/150‑1 (MQ=255)
                                                            cGCCCTTTGCGAAATGCCCATACATGATTGTGTTTTGTGTTGATGAGTGTTCATGCGCGAATCAACTTTTCGTCGGGCTGCAACTGGAATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTaa                                                                                           <  1:2694886/150‑1 (MQ=255)
                                                            cGCCCTTTGCGAAATGCCCATACATGAATGAGTATTGAGTTGTTGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGATTTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTaa                                                                                           <  1:3148705/150‑1 (MQ=255)
                                                             gCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAg                                                                                          <  2:2671906/150‑1 (MQ=255)
                                                               ccTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAATTTTGCCAGTTGCAGCGAGTGTATGACAGCCGAGAATTCAAGGGAAAAGGTGTTATCTAAggg                                                                                        >  1:1025023/1‑149 (MQ=255)
                                                               acTTTGCTAAATGCCAATACATCAATGATCATTCATTTGATGAGTTTTCATCACCGAATAAACTTTTCGTCGGGCTGAAACTGGCATTTTTGCCATTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGt                                                                                        <  1:1399638/149‑1 (MQ=255)
                                                                cTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTCCAGTTTCAGCGAGTGTATTACAGCCAAGAATTAAAAGGAAAATGTGTTGTCTAAGGta                                                                                       >  1:2326056/1‑149 (MQ=255)
                                                                 tttGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAGTTTTGCCAGTTGCAGCGAGTGTGTGACAGCCGAGAATTAAAGGGGAATGGTGTTCTCTAAGGtgg                                                                                      >  2:774651/1‑148 (MQ=255)
                                                                 tttGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCACTTTCAGCGAGTTTATTACAGCCAAGAAATAAATGGAAAATGTGTTATATAATGTCa                                                                                      >  1:4227967/1‑150 (MQ=255)
                                                                 tttGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCGGCAACATGCAATTTTTGCAGTTGCAGCGAGGGTTTTTCAGCCGAGAAATCAAGGGAAAATGTGTTGTGTAAgggta                                                                                      >  1:1230669/1‑147 (MQ=255)
                                                                       aaaTGCCCATACATGACTGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATc                                                                                <  2:468603/150‑1 (MQ=255)
                                                                           gCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCGAGAATTAAAGGGAAATGGTGTTCTTTAAGGTCATTAATCGTGa                                                                            >  2:3354386/1‑150 (MQ=255)
                                                                                 aCATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTTCCAGTTTCAACGCGTGTATTAAAGCCGAGAAATAAAGGGAAAAGTTGTTGTATAAGGTAATTAATAGTGaaaactt                                                                      >  1:3911034/1‑146 (MQ=255)
                                                                                 aCATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAAAAGGCCTTTTTTCCAGTTGCAGCGAGTGTATTAAAGCCAAGAAATAAAGGGAAATGGTTTTGTATAAGGTCATTAATAGTGaaaactt                                                                      >  1:909317/1‑146 (MQ=255)
                                                                                  cATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTc                                                                     <  2:1273868/150‑1 (MQ=255)
                                                                                   atgaatgaGGATTGGGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCACTTTCAGCGAGTGTATGACAGCCTAGAATTAAAGGGAAATGGTTTTATGTAATGTCATTAATAGTGaaaacttag                                                                    >  2:2456010/1‑144 (MQ=255)
                                                                                   atgaatgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACATGCATTTTTTACACTTTCAGCCAGTGTATGAAAGCCAAGAATTAAAGGGAAAAGGTGTTATATAATGTCATTAAACGTGAAAAGttag                                                                    >  1:2617820/1‑148 (MQ=255)
                                                                                     gaatgaGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCAGCAACATGCATTTTTGCCAGGTGCAGCGAGTGTATGACAGCCGAGAATACAAGGGAAAAGGTGGTGTCTAAGGTCATTAAAAGTGAAAAGTTCGCa                                                                  >  2:2091212/1‑150 (MQ=255)
                                                                                      aatgaGTATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCaa                                                                 <  1:1304475/150‑1 (MQ=255)
                                                                                            gATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCCTTTTTGCCAGTTGCAGCGAGTGTATGACACCCTAGAATTCAAGGGAAATGGTGTTATCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTac                                                           >  2:2596765/1‑150 (MQ=255)
                                                                                            gATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAAAGGTGTTCTATAAGGTCATTAATCGTGAAACGTTCGCAAGCGTac                                                           >  2:3681746/1‑150 (MQ=255)
                                                                                              ttGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCAATTTTTCCAGTTTCAGCCAGTGTATGAAAGCCTAGAAATCAAGGGAAATGGTGTTCTATAAGGTCATTAATAGTGAAAACTTCCAAAGCGTacac                                                         >  2:3896081/1‑150 (MQ=255)
                                                                                              ttGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAATTTTTCCAGTTGCAGCGAGTGTATGACACCCTATAATTCAAGGGAAATTGTGTTCTCTAAGGTCAATAATAGTGAAAAGTTAGCAAGCGTacac                                                         >  2:3585432/1‑150 (MQ=255)
                                                                                                    ttatgaGTGTTCTTGAGCTAATCAACGTTTGGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTa                                                   <  1:758409/148‑1 (MQ=255)
                                                                                                      atgaGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCCACGGGCAGTTTTGCCAGTTTCAGCGAGTGTATGACAGCCGAGAATTCAAGGGCAATTGTGTTGTCTAAGGTGATTAATAGTGTAACGTTCGCAAGCGTAACCGCCCTTCg                                                 >  2:2397376/1‑150 (MQ=255)
                                                                                                        gaGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCGAGAATTCAAGGGAACTGGTGTTCTCTAAGGTCATTAATCGTGAAAAGTTCGCAAGCGTACACGCCCAACGgc                                               >  2:712798/1‑150 (MQ=255)
                                                                                                                aTGAGCTAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCTGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTg                                       <  1:4047294/150‑1 (MQ=255)
                                                                                                                aTGAGCGAATCAACGTTTCGTCGGGCTGCAACAGGCAATTTTTCCCCTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTAAGTAAGGTCAATAAAAGGGTAAAGTTAGCAAGCCTACAAGCCCTACCGCGCGAGTTg                                       >  2:4229509/1‑150 (MQ=255)
                                                                                                                 tGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGGTTCAGCGAGTGTATGACACCCTAGAATTCAAGGGAAATGGTGTTGTCTAAGGTCATTAAACGTGAAACGTTAGCAATCGTACAAGCCCAAACGCGCGATTTgg                                      >  2:3262175/1‑150 (MQ=255)
                                                                                                                 tGAGCGAATCAACGTTTCGTCGGGCTGCAAAAGGCATTTTTTCCACTTTCAGCGAGTGTAATAAAGCCTAGAATTAAAGGGAAAAGGTGTTGTATAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCAAAGgcgcggttttg                                      >  2:3020401/1‑143 (MQ=255)
                                                                                                                  ttgCGAATCAACTTTTCGTCGGGCTGAAACTGGCATTTTTGCCTTTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTggg                                     <  2:409123/148‑1 (MQ=255)
                                                                                                                  gAGCGATTCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTggg                                     <  2:2278290/150‑1 (MQ=255)
                                                                                                                       attCAACGTTTCGTCGGGATGAAACTGTCTTTTTTGCCTGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGAt                                <  1:148189/148‑1 (MQ=255)
                                                                                                                       aaTCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGAt                                <  1:3624211/150‑1 (MQ=255)
                                                                                                                       aaTCAACGTTTCGTCGGGCTGCAACAGGCAATTTTTCCACTTGCAGCGAGTGTATGGAAGCCAAGAAATAAAAGGAAAATGTGGTATCTAAGGTCATTAAAAGGGAAACCTTCGCAAGCGTACACGCCCAACGGCGCGATTTGGGGTGAt                                >  1:1443907/1‑150 (MQ=255)
                                                                                                                        aTCAACGTTTCGTCGGGCTGCAACTGGCATTTTTTCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTGAATGGAAAATGTGTTATCTAAGGTCATTAATCGTGAAAAGTTAGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATg                               >  2:4386518/1‑150 (MQ=255)
                                                                                                                            aCGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCa                           <  2:1035612/150‑1 (MQ=255)
                                                                                                                                     cgGGCTGCAACTGGCTTTTTGGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCgt                  <  1:72024/150‑1 (MQ=255)
                                                                                                                                      gggCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCgtg                 >  2:2798013/1‑150 (MQ=255)
                                                                                                                                         cTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTg              >  2:1062952/1‑150 (MQ=255)
                                                                                                                                          tGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGa             >  2:1348179/1‑150 (MQ=255)
                                                                                                                                          tGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGa             >  2:2835955/1‑150 (MQ=255)
                                                                                                                                          tGCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGa             >  1:2359547/1‑150 (MQ=255)
                                                                                                                                          tGCAACTGGCAAAAATGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGa             >  2:1635720/1‑150 (MQ=255)
                                                                                                                                           gattcTGGCTTTTTTGCTTTTTGTTGTGTGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGaa            <  1:2752493/146‑1 (MQ=255)
                                                                                                                                           gCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGaa            <  1:2217638/150‑1 (MQ=255)
                                                                                                                                           gCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGaa            <  1:4058537/150‑1 (MQ=255)
                                                                                                                                           gCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGaa            <  1:1350373/150‑1 (MQ=255)
                                                                                                                                           gCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGaa            <  2:1323915/150‑1 (MQ=255)
                                                                                                                                           gCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGaa            <  1:1010045/150‑1 (MQ=255)
                                                                                                                                                tGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACa       >  2:3799923/1‑150 (MQ=255)
                                                                                                                                                tGGCATGTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACa       >  2:2636025/1‑150 (MQ=255)
                                                                                                                                                 ggCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACaa      <  2:407181/150‑1 (MQ=255)
                                                                                                                                                    aTTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATg                                                                                                      >  2:2892931/1‑51 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACGCGCAGAGCCTGTTGCTGATGAACTTCCTGCCGGTGCAGTAAGAGGCCTCTGCCGCACGCCCTTTGCGAAATGCCCATACATGAATGAGGATTGAGTTGATGAGTGTTCATGAGCGAATCAACGTTTCGTCGGGCTGCAACTGGCATTTTTGCCAGTTGCAGCGAGTGTATGACAGCCTAGAATTCAAGGGAAATGGTGTTCTCTAAGGTCATTAATCGTGAAACGTTCGCAAGCGTACACGCCCTACGGCGCTATTTGGGGTGATGGCCATTGCTCCGTGTTGAAGTACAATGGT  >  NC_002947/3580244‑3580542

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: