Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,562,171 C→A 8.6% intergenic (‑81/‑88) pctA ← / → PP_1372 methyl‑accepting chemotaxis protein PctA/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,562,1710CA 8.6% 227.6 / 9.5 78intergenic (‑81/‑88)pctA/PP_1372methyl‑accepting chemotaxis protein PctA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/33);  new base A (4/3);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.47e-01

AGGCGAGTATCACGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGGCCAGCCAGTCGGGCAAG  >  NC_002947/1562024‑1562319
                                                                                                                                                   |                                                                                                                                                    
aGGCGAGTATCACGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGAACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTccc                                                                                                                                                     >  3:252899/1‑149 (MQ=255)
  gCGAGTATCACGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCtt                                                                                                                                                   >  5:237833/1‑149 (MQ=255)
     aGTATCACGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTc                                                                                                                                                <  6:237833/149‑1 (MQ=255)
      gTATCACGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCg                                                                                                                                               <  4:218998/149‑1 (MQ=255)
      gTATCACGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCg                                                                                                                                               <  1:69632/149‑1 (MQ=255)
      gTATCACGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGGGACTTTTCg                                                                                                                                               <  7:12565/149‑1 (MQ=255)
       tATCACGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGGGACTTTTCgg                                                                                                                                              <  7:178066/149‑1 (MQ=255)
        aTCACGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCgg                                                                                                                                              <  4:63343/148‑1 (MQ=255)
         tCACGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCt                                                                                                                                            >  4:31574/1‑149 (MQ=255)
           aCGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGt                                                                                                                                          >  4:106595/1‑149 (MQ=255)
             gATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGtt                                                                                                                                         >  4:162539/1‑148 (MQ=255)
             gATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGtt                                                                                                                                         >  4:167941/1‑148 (MQ=255)
                      tGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTAtt                                                                                                                               <  6:277404/149‑1 (MQ=255)
                       gccgccAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTAtt                                                                                                                               >  2:227384/1‑148 (MQ=255)
                        ccgccAAAAGAATATTGTTGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTaa                                                                                                                             <  3:207145/149‑1 (MQ=255)
                          gccAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCgg                                                                                                                                              >  7:117462/1‑130 (MQ=255)
                          gccAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCgg                                                                                                                                              <  8:117462/130‑1 (MQ=255)
                             aaaaGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCg                                                                                                                         >  2:259757/1‑148 (MQ=255)
                             aaaaGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGa                                                                                                                        <  4:307258/149‑1 (MQ=255)
                             aaaaGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATACCGAAAAGTACCTTTTAGGCTGTTACTCTCTATTAATGCg                                                                                                                         >  3:292856/1‑148 (MQ=255)
                              aaaGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGaa                                                                                                                       <  3:122218/149‑1 (MQ=255)
                              aaaGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGaa                                                                                                                       <  2:301001/149‑1 (MQ=255)
                              aaaGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGaa                                                                                                                       <  1:259757/149‑1 (MQ=255)
                               aaGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAAt                                                                                                                      <  3:162539/149‑1 (MQ=255)
                               aaGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAAt                                                                                                                      <  3:167941/149‑1 (MQ=255)
                                        gtgGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACAcc                                                                                                                <  5:273666/146‑1 (MQ=255)
                                        gtgGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCca                                                                                                              <  3:31574/148‑1 (MQ=255)
                                           gCTGAAACGAAGGCTTTTGTTCATTGCGTAATCGGTTCCGCTAAGGTTTTATATCAGTCGTTCGCACGAGTGGATCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAg                                                                                                           <  5:304403/148‑1 (MQ=255)
                                                   gAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCACGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGcc                                                                                                  >  5:337237/1‑149 (MQ=255)
                                                         ttttGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGGGACTTTTc                                                                                                                                                >  6:315931/1‑97 (MQ=255)
                                                         ttttGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGGGACTTTTc                                                                                                                                                <  5:315931/97‑1 (MQ=255)
                                                                caaggCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGAGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTAATAAGGCGAATCACACCCACAGCTATCGGCAGAATTAAGCTaaa                                                                                     >  6:86149/4‑149 (MQ=255)
                                                                cATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTa                                                                                       >  7:256735/1‑147 (MQ=255)
                                                                  tGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGc                                                                                   <  2:204289/149‑1 (MQ=255)
                                                                           cGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCctctct                                                                                                                                  >  5:249386/1‑93 (MQ=255)
                                                                           cGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCctctct                                                                                                                                  <  6:249386/93‑1 (MQ=255)
                                                                            ggTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCaga                                                                          <  3:22192/148‑1 (MQ=255)
                                                                                   cTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGc                                                                                   <  2:348849/132‑1 (MQ=255)
                                                                                   cTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGc                                                                                   >  1:348849/1‑132 (MQ=255)
                                                                                   cTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTaaa                                                                  <  1:227384/149‑1 (MQ=255)
                                                                                   cTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTaaa                                                                  >  8:147927/1‑149 (MQ=255)
                                                                                   cTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTa                                                                    >  1:58642/1‑147 (MQ=255)
                                                                                    tAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTa                                                                    >  6:204027/1‑146 (MQ=255)
                                                                                      aGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCTAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTaaa                                                                  <  5:204027/146‑1 (MQ=255)
                                                                                      aGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAagag                                                               <  2:58642/149‑1 (MQ=255)
                                                                                         ttttATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGa                                                                                                <  6:5049/113‑1 (MQ=255)
                                                                                         ttttATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGa                                                                                                >  5:5049/1‑113 (MQ=255)
                                                                                             atatCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCg                                                         <  8:256735/148‑1 (MQ=255)
                                                                                             atatCAGGCGTTCGAACGAGTGGAGCCGTAGGTTGTGCGACATCCCGAAAAGTCACATTTCGTCTGTTCCTCTCTATTAAg                                                                                                                            >  4:287915/1‑81 (MQ=255)
                                                                                                tCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAAt                                                     >  2:180290/1‑149 (MQ=255)
                                                                                                  aGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGACGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAAttt                                                   >  2:291815/1‑149 (MQ=255)
                                                                                                    gCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTc                                                 >  4:115430/1‑149 (MQ=255)
                                                                                                    gCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTc                                                 <  6:337237/149‑1 (MQ=255)
                                                                                                      gTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCgg                                                                                                                                              >  4:222340/1‑54 (MQ=255)
                                                                                                      gTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCgg                                                                                                                                              <  3:222340/54‑1 (MQ=255)
                                                                                                      gTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGa                                               >  3:75602/1‑149 (MQ=255)
                                                                                                         cGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTaa                                                                                           >  6:151039/1‑102 (MQ=255)
                                                                                                          gCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCAt                                           >  3:177695/1‑149 (MQ=255)
                                                                                                              gAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGtt                                       <  7:147927/149‑1 (MQ=255)
                                                                                                                   gAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAgg                                  >  8:268795/1‑149 (MQ=255)
                                                                                                                   gAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTATTCGGCGGTTCCTCGCTATTAAGGGGAATAACACCCACAGGGATCGGCCGCATTAAGCGAAAGCTTGAGCAGAGGATTAAAGAGCTGTCGGAATTTTCGACCATAGTTGTAgg                                  >  1:235644/1‑149 (MQ=255)
                                                                                                                         tAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCT‑AAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGc                           >  3:199480/1‑149 (MQ=255)
                                                                                                                                tgCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGAt                     <  7:268795/149‑1 (MQ=255)
                                                                                                                                 gCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACGATCGTTGTAGGCGCGGGCCCGGATg                    >  4:147701/1‑149 (MQ=255)
                                                                                                                                   aaCATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATgg                   <  4:292856/148‑1 (MQ=255)
                                                                                                                                     cATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGgcca                >  8:246307/1‑149 (MQ=255)
                                                                                                                                      aTCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGgccag               >  6:61648/1‑149 (MQ=255)
                                                                                                                                      aTCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGgc                  <  3:65857/146‑1 (MQ=255)
                                                                                                                                      aTCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGgc                  >  4:65857/1‑146 (MQ=255)
                                                                                                                                       tCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGgccagc              >  7:142987/1‑149 (MQ=255)
                                                                                                                                         ccGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCTAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGgccagc              >  1:300517/1‑147 (MQ=255)
                                                                                                                                         ccGAAAAGGGACTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGgccagc              >  6:242274/1‑147 (MQ=255)
                                                                                                                                         ccGAAAAGGGACTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGgccagc              >  6:237058/1‑147 (MQ=255)
                                                                                                                                          cGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCCGCCGAATTAAGCTAAAGCTTGAGCAGAGGa                                                                       <  1:29742/89‑1 (MQ=255)
                                                                                                                                          cGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACACACAGGTATCCGCCGAATTAAGCTAAAGCTTGAGCAGAGGa                                                                       >  2:29742/1‑89 (MQ=255)
                                                                                                                                           gAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGGCCAGCCAGt          >  1:44314/1‑149 (MQ=255)
                                                                                                                                                gTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGGCCAGCCAGTCGGgc     <  1:180290/149‑1 (MQ=255)
                                                                                                                                                 tCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGTGCGGGCCCGGATGGCCAGCCAGTCGGgc     >  7:270834/1‑148 (MQ=255)
                                                                                                                                                 tCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGGCCAGCCAGTCGGgc     >  2:274367/1‑148 (MQ=255)
                                                                                                                                                   ccTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGGCCAGCCAGTCGGgcaag  >  1:358077/1‑149 (MQ=255)
                                                                                                                                                   ccTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGGCCAGCCAGTCGGgcaag  >  6:48770/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AGGCGAGTATCACGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGGCCAGCCAGTCGGGCAAG  >  NC_002947/1562024‑1562319

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑