Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,628,993 T→G 8.6% M168R (ATG→AGG) ‡ mucB → sigma factor AlgU regulator MucB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,628,9930TG 8.6% 205.6 / 9.5 70M168R (ATG→AGG) ‡mucBsigma factor AlgU regulator MucB
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base G (4/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

ATGGGTTGGTACGATTTGAGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAATGATGCC  >  NC_002947/1628845‑1629132
                                                                                                                                                    |                                                                                                                                           
aTGGGTTGGTACGATTTGAGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGAt                                                                                                                                             <  6:331360/149‑1 (MQ=255)
  gggTTGGTACGATTTGAGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGAACAGCATAGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGc                                                                                                                                           >  5:337022/1‑149 (MQ=255)
         tACGATTTGAGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGa                                                                                                                                    >  6:366635/1‑149 (MQ=255)
             aTTTGAGTGTGGCAGGCAAGTCACGGGTCGCCGGGTGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAg                                                                                                                                >  4:63329/1‑149 (MQ=255)
               ttGAGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCATTGATGCTCAACGACAAggg                                                                                                                              >  1:230265/1‑149 (MQ=255)
                tGAGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGAGCATCAACGaacgcag                                                                                                                              <  6:336533/148‑7 (MQ=255)
                 gAGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAgggg                                                                                                                             <  2:12370/148‑1 (MQ=255)
                  aGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGGTCGTTGATGCTCAACGACAAGGGGCa                                                                                                                           <  7:131971/149‑1 (MQ=255)
                        gcaggcaAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGc                                                                                                                      <  6:282017/148‑1 (MQ=255)
                         caggcaAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTg                                                                                                                    <  5:207775/149‑1 (MQ=255)
                         caggcaAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTg                                                                                                                    <  4:125175/149‑1 (MQ=255)
                          agggaAGTCACTGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTgg                                                                                                                   <  8:155596/149‑1 (MQ=255)
                           ggcaAGTCACGGGTCTCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGa                                                                                                                  <  5:50529/149‑1 (MQ=255)
                            gcaAGGCACGGGTCGCCGGTCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGaa                                                                                                                 <  3:214585/149‑1 (MQ=255)
                                 tCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCtt                                                                                                            <  5:137038/149‑1 (MQ=255)
                                   aCGGGTCGCCGGGCGCGATGCCGTGATTATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTc                                                                                                           >  5:262081/1‑148 (MQ=255)
                                     gggTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGAt                                                                                                                                             >  2:296264/1‑112 (MQ=255)
                                     gggTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGAt                                                                                                                                             <  1:296264/112‑1 (MQ=255)
                                      ggTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGa                                                                                                       <  2:298629/149‑1 (MQ=255)
                                       gTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGAt                                                                                                      >  6:331172/1‑149 (MQ=255)
                                          gCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGa                                                                                                    <  5:366635/148‑1 (MQ=255)
                                          gCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGAc                                                                                                   <  5:322877/149‑1 (MQ=255)
                                                cgcgATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGtct                                                                                             <  1:98074/149‑1 (MQ=255)
                                                       ccGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATAc                                                                                       <  3:63329/148‑1 (MQ=255)
                                                        cGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATAcc                                                                                      >  4:306102/1‑148 (MQ=255)
                                                         gTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGa                                                                                    >  6:334853/1‑149 (MQ=255)
                                                             tCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGa                                                                                 <  7:158886/148‑1 (MQ=255)
                                                                ttACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATa                                                                             >  2:352752/1‑149 (MQ=255)
                                                                    gCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATAcc                                                                           >  5:2537/1‑147 (MQ=255)
                                                                    gCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCa                                                                         <  5:331172/149‑1 (MQ=255)
                                                                        acgccacgcGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGa                                                                     >  6:110206/1‑149 (MQ=255)
                                                                         cgccacgcGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCATTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGAc                                                                    <  2:230265/149‑1 (MQ=255)
                                                                          gccacgcGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACg                                                                   >  5:97034/1‑149 (MQ=255)
                                                                            cacgcGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACg                                                                   >  7:180/1‑147 (MQ=255)
                                                                                 gACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACtt                                                            <  8:180/149‑1 (MQ=255)
                                                                                 gACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACtc                                                            >  4:338111/1‑148 (MQ=255)
                                                                                         tCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAg                                                    <  1:352752/149‑1 (MQ=255)
                                                                                             tctGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAg                                                    <  1:320652/143‑1 (MQ=255)
                                                                                               tGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGtgct                                               <  3:306102/148‑1 (MQ=255)
                                                                                                 cGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGtgctgc                                             <  3:67251/148‑1 (MQ=255)
                                                                                                 cGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGtgctgc                                             <  2:338086/148‑1 (MQ=255)
                                                                                                     cGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTGCAGATGACCCGTCTCGATACCGATCAAATACCCAGCGACGCAGACTTGCGCCCAGGTGctgcctgc                                         >  3:214644/1‑148 (MQ=255)
                                                                                                      gAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCaa                                       <  6:337022/149‑1 (MQ=255)
                                                                                                            cACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGt                                 >  5:364129/1‑149 (MQ=255)
                                                                                                               cTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTg                                >  5:353005/1‑147 (MQ=255)
                                                                                                                  gACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTgcagca                           >  4:231640/1‑149 (MQ=255)
                                                                                                                  gACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTgcagca                           >  2:313054/1‑149 (MQ=255)
                                                                                                                  gACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTCGCGCCCAGGTGCTGCCTGCAAGCCAGTgcagc                            <  3:338111/148‑1 (MQ=255)
                                                                                                                   aCCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCAt                          >  3:185458/1‑149 (MQ=255)
                                                                                                                    ccGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGCCCCGTCTCGATACCGATGAAATACCAAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCaacc                                     >  6:10364/1‑135 (MQ=255)
                                                                                                                    ccGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCAtg                         <  1:313054/149‑1 (MQ=255)
                                                                                                                    ccGTCACTCGGGTGAGCCGCTGCGTGCGTTCATCCTCACCGACAAGGGGCAGTTGCGGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCCACGCAGACTTGCGCCCAGGAGCTGCCTGCAAGc                                     <  5:10364/137‑1 (MQ=255)
                                                                                                                      gTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCAtgtg                       >  1:48054/1‑149 (MQ=255)
                                                                                                                          cTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAg                   >  8:51247/1‑149 (MQ=255)
                                                                                                                           tCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGc                  >  5:326144/1‑149 (MQ=255)
                                                                                                                           tCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGc                  >  6:347859/1‑149 (MQ=255)
                                                                                                                             gggTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCg                >  5:19984/1‑149 (MQ=255)
                                                                                                                             gggTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCg                >  2:188842/1‑149 (MQ=255)
                                                                                                                             gggTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCg                >  5:14714/1‑149 (MQ=255)
                                                                                                                             gggTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCg                >  2:370544/1‑149 (MQ=255)
                                                                                                                                    ccGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCaa         <  6:2537/149‑1 (MQ=255)
                                                                                                                                      gCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCACGATACCGATGAAATAACCAGCGACGCAGCCTTGCGCCCAGGTGCTGCCTGCAAGCGAGTGCAGCATGTGGCAGCCGGCAGCAat        >  7:249076/1‑148 (MQ=255)
                                                                                                                                        ccttgtcgttGAGCATCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAatgat     >  3:332718/6‑149 (MQ=255)
                                                                                                                                         cttgtcgttGAGCATCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAatgatg    >  1:157212/5‑149 (MQ=255)
                                                                                                                                          ttgtcgttGAGCATCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTgg                      <  4:225610/127‑1 (MQ=255)
                                                                                                                                          ttgtcgttGAGCATCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTgg                      >  3:225610/4‑130 (MQ=255)
                                                                                                                                          ttgtcgttGAGCATCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAATGATGc   >  5:155424/4‑149 (MQ=255)
                                                                                                                                          cgttcgttGATGCTCAACGACAAGGGGCAGTTGCTGTAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAatgat     <  6:262081/147‑1 (MQ=255)
                                                                                                                                           gttcgttGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAATGATGcc  <  6:353005/149‑1 (MQ=255)
                                                                                                                                               gttGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTgca                              <  7:321650/117‑1 (MQ=255)
                                                                                                                                               gttGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTgca                              >  8:321650/1‑117 (MQ=255)
                                                                                                                                                  gATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGgcag            <  3:272606/132‑1 (MQ=255)
                                                                                                                                                  gATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGgcag            >  4:272606/1‑132 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
ATGGGTTGGTACGATTTGAGTGTGGCAGGCAAGTCACGGGTCGCCGGGCGCGATGCCGTGATCATTACGCTCACGCCACGCGACCAGCATCGCTATGCGTTCGAATTGCACCTGGACCGTCACTCGGGTTTGCCGCTGCGTTCGTTGATGCTCAACGACAAGGGGCAGTTGCTGGAACGCTTCCAGATGACCCGTCTCGATACCGATGAAATACCCAGCGACGCAGACTTGCGCCCAGGTGCTGCCTGCAAGCCAGTGCAGCATGTGGCAGCCGGCAGCAATGATGCC  >  NC_002947/1628845‑1629132

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