Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,233,594 A→T 8.3% D137V (GAC→GTC)  PP_4611 → ECF family RNA polymerase sigma‑70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,233,5940AT 8.3% 210.8 / 6.7 70D137V (GAC→GTC) PP_4611ECF family RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base A (33/31);  new base T (3/3);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTG  >  NC_002947/5233449‑5233740
                                                                                                                                                 |                                                                                                                                                  
cTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGc                                                                                                                                                >  1:698853/1‑150 (MQ=255)
 tGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGcc                                                                                                                                               <  2:1259712/150‑1 (MQ=255)
   cAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTg                                                                                                                                             <  2:1266007/150‑1 (MQ=255)
    aaGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTgg                                                                                                                                            <  2:698853/150‑1 (MQ=255)
    aaGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTgg                                                                                                                                            <  1:1286397/150‑1 (MQ=255)
         cTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCag                                                                                                                                       >  2:1266554/1‑149 (MQ=255)
          tCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGACAAAAATAATCGGGGGATTGCTAGCGCTAGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAg                                                                                                                                      >  2:349613/1‑150 (MQ=255)
          ccGCAAAGGCCCCCGAAGAGGTCGCCCCGTGTGAGGAACAGAAAGCGATCATCTTCGAGACATTGCTACTGCGAGACAACTCGCTGGACGGCCTTCCACCGTTGCTCAAACGAGCCTTCCGGATGGCCCAGGCCGTAAACCTGGGCCagc                                                                                                                                      <  2:494592/149‑3 (MQ=255)
           cGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAgg                                                                                                                                     <  1:1158372/150‑1 (MQ=255)
            gCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGCCATTGCTAGCGCGCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGc                                                                                                                                    >  1:884975/1‑150 (MQ=255)
              acaGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGa                                                                                                                                  >  1:532700/1‑150 (MQ=255)
               caGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCTAGGCCGTCTTCCTGGGCCagcatgaa                                                                                                                                 <  2:88201/150‑8 (MQ=255)
                aGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGaaa                                                                                                                                >  2:1352633/1‑150 (MQ=255)
                aGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGaaa                                                                                                                                <  2:532700/150‑1 (MQ=255)
                        gAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAAtcgctcgc                                                                                                                        >  2:1273219/1‑150 (MQ=255)
                               tCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACt                                                                                                                   >  2:1139049/1‑148 (MQ=255)
                               tCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACt                                                                                                                   <  1:1139049/148‑1 (MQ=255)
                                 gcgcCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGc                                                                                                               <  1:1352633/150‑1 (MQ=255)
                                          gAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTg                                                                                                      <  1:1273219/150‑1 (MQ=255)
                                             gAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGcc                                                                                                   >  1:274061/1‑150 (MQ=255)
                                                cAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACg                                                                                                <  1:755178/150‑1 (MQ=255)
                                                   aaGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTc                                                                                             <  1:212301/150‑1 (MQ=255)
                                                    aGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCa                                                                                            >  2:393570/1‑150 (MQ=255)
                                                             tcTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTAcc                                                                                   >  1:1083221/1‑150 (MQ=255)
                                                                 cGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGaa                                                                               <  1:780163/150‑1 (MQ=255)
                                                                 cGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCAACCCCCGTGGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTAACTGaa                                                                               >  2:410501/1‑150 (MQ=255)
                                                                  gagaCATTGCTAGCGCTCGACCACATGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGCGAGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTCGGCGTTTCGCTGGGCACGGTCAAGCGGTACCAGAAc                                                                              >  1:317944/1‑150 (MQ=255)
                                                                       aTTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGc                                                                         >  2:1143417/1‑150 (MQ=255)
                                                                          gCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGt                                                                                              >  2:576278/1‑126 (MQ=255)
                                                                          gCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGt                                                                                              <  1:576278/126‑1 (MQ=255)
                                                                              gcgcTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCGGCTGGGGCAGGTCGACGGCATGGGAAAAGGTAATATCGGCCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGCACAAGGCAGCCATg                                                                  >  1:1391154/1‑150 (MQ=255)
                                                                                 cTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATgcgc                                                               <  2:1083221/150‑1 (MQ=255)
                                                                                  tCGACCACTTGCTGGTCGGCCTTCCACCGGGGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATGGCTCGCGAACTGGGCATTTAGCTGGCCAGGGGCAAGCGCTACCTGAACGAGGCAGCCATGCgct                                                              >  2:943354/1‑150 (MQ=255)
                                                                                  tCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCAACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgct                                                              <  2:230803/150‑1 (MQ=255)
                                                                                    gcccACTTGCTGGACACGCTGTCACCGTTTGTGAAACGCGCCTTCCTGCTGGCCCAGGTCGAAGCGGCAAACCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGTCATGCgctgc                                                            <  1:667867/148‑1 (MQ=255)
                                                                                          ttGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAgg                                                                          <  2:485071/130‑1 (MQ=255)
                                                                                          ttGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAgg                                                                          >  1:485071/1‑130 (MQ=255)
                                                                                             cTGGACGGCCTTCCCCCGGTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAATGCAAAATCGCTCGCGAACTGGGCATTGCGCCGGCCACTGTCAAGCGCTACCTGAACAAGGCCGCCCTGCGCTGCTACTTcgcg                                                   >  2:84718/1‑150 (MQ=255)
                                                                                              tGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgctgc                                                            <  2:999740/140‑1 (MQ=255)
                                                                                              tGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgctgc                                                            >  1:999740/1‑140 (MQ=255)
                                                                                                 aCGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtgt                                               <  2:1391154/150‑1 (MQ=255)
                                                                                                 aCGGCCTTCCACCGTTGGTCAAACGCGCCTTCCCGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCtgt                                               <  1:672065/150‑1 (MQ=255)
                                                                                                     ccTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCAttt                                                                                                           <  1:682757/86‑1 (MQ=255)
                                                                                                     ccTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCAttt                                                                                                           >  2:682757/1‑86 (MQ=255)
                                                                                                            ccGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGt                                    <  2:884975/150‑1 (MQ=255)
                                                                                                                  gTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGATATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGtt                                   <  2:795243/145‑1 (MQ=255)
                                                                                                                  gTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgctgc                                                            <  2:951585/120‑1 (MQ=255)
                                                                                                                  gTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCgctgc                                                            >  1:951585/1‑120 (MQ=255)
                                                                                                                  gTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGtt                                   >  1:795243/1‑145 (MQ=255)
                                                                                                                  gTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTccc                              >  2:775217/1‑150 (MQ=255)
                                                                                                                   tCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTcccc                             >  2:496207/1‑150 (MQ=255)
                                                                                                                         gcgcCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCtggt                       <  1:393570/150‑1 (MQ=255)
                                                                                                                         gcgcCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTCCGCGCTGTGAACGCCTCGTTTTCCCCGCtggt                       >  2:206403/1‑150 (MQ=255)
                                                                                                                             cTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCg                   >  2:1175938/1‑150 (MQ=255)
                                                                                                                                  tgctgGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTa                                                                                     <  1:621237/79‑1 (MQ=255)
                                                                                                                                  tgctgGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTa                                                                                     >  2:621237/1‑79 (MQ=255)
                                                                                                                                   gctgGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCATGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAg             <  1:943354/150‑1 (MQ=255)
                                                                                                                                   cttgGCCCAGGCCGTCGACCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGc                    <  1:531033/141‑1 (MQ=255)
                                                                                                                                   cttgGCCCAGGCCGTCGACCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGc                    >  2:531033/3‑143 (MQ=255)
                                                                                                                                   cttgGCCCAGGCCGTCGACCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAg             >  2:621222/3‑150 (MQ=255)
                                                                                                                                   cttgGCCCAGGCCGTCGACCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAg             >  1:771121/3‑150 (MQ=255)
                                                                                                                                    ctgGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAgg            <  1:496207/150‑1 (MQ=255)
                                                                                                                                        cccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt        <  1:206403/150‑1 (MQ=255)
                                                                                                                                        cccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCAGGACGGGTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGt        <  1:1143417/150‑1 (MQ=255)
                                                                                                                                         ccATGTCGACGGCCTGGGCCAAGACTAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTg       >  2:22568/1‑150 (MQ=255)
                                                                                                                                         ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTg       >  1:639550/1‑150 (MQ=255)
                                                                                                                                         ccAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTg       >  1:209465/1‑150 (MQ=255)
                                                                                                                                          cAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGaa                <  2:654182/140‑1 (MQ=255)
                                                                                                                                          cAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGaa                >  1:654182/1‑140 (MQ=255)
                                                                                                                                              tCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTg  >  2:392899/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CTGCAAGCGCTCGCACAGGTCCCGGAAGAGGTCGCGCCGTGCGAGGAACAGAAGGCGATCATCTTCGAGACATTGCTAGCGCTCGACCACTTGCTGGACGGCCTTCCACCGTTGGTCAAACGCGCCTTCCTGCTGGCCCAGGTCGACGGCCTGGGCCAAGGCGAAATCGCTCGCGAACTGGGCATTTCGCTGGCCACGGTCAAGCGCTACCTGAACAAGGCAGCCATGCGCTGCTACTTCGCGCTGTGAACGCCTCGTTTTCCCCGCTGGTGGCGGAACAGGCCGTGCATTG  >  NC_002947/5233449‑5233740

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