Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,557,923 C→A 7.3% intergenic (‑294/+37) PP_4888 ← / ← purA methyl‑accepting chemotaxis transducer/adenylosuccinate synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,557,9230CA 7.3% 206.9 / 6.2 77intergenic (‑294/+37)PP_4888/purAmethyl‑accepting chemotaxis transducer/adenylosuccinate synthetase
Reads supporting (aligned to +/- strand):  ref base C (32/38);  new base A (4/2);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

GCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGC  >  NC_002947/5557776‑5558068
                                                                                                                                                    |                                                                                                                                                 
gCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTtct                                                                                                                                                  <  1:260334/149‑1 (MQ=255)
 cTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTtctc                                                                                                                                                 <  3:103956/149‑1 (MQ=255)
    cTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAgg                                                                                                                                              >  5:343786/1‑149 (MQ=255)
    cTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAg                                                                                                                                               <  2:133365/148‑1 (MQ=255)
      aCAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGc                                                                                                                                            >  1:267380/1‑149 (MQ=255)
          cggccggcGCCCAACCATCGTTGGGTGCCCACTTCGCCGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCcgcg                                                                                                                                        <  4:211717/149‑1 (MQ=38)
             ccTGCGCCCAACCATCGTTGGGGCCCCACTTAGTAGGGGGTTTACCCGGGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCgcggcg                                                                                                                                     <  7:165260/149‑1 (MQ=38)
              cTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTGTCAGGGCCgcggcg                                                                                                                                     >  6:71492/1‑148 (MQ=255)
               tGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGt                                                                                                                                    >  5:167159/1‑148 (MQ=255)
               tGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGAGAAGGGCcgcg                                                                                                                                        >  1:270757/1‑144 (MQ=255)
               tGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGAGAAGGGCcgcg                                                                                                                                        <  2:270757/144‑1 (MQ=255)
                 cgcCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAg                                                                                                                                 <  2:267380/149‑1 (MQ=255)
                   cccAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCGGCGGCGTCAGAt                                                                                                                               <  6:167159/149‑1 (MQ=255)
                      aaCCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAg                                                                                                                            >  1:277907/1‑149 (MQ=255)
                       aCCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATc                                                                                                                              <  6:203016/146‑1 (MQ=255)
                       aCCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATc                                                                                                                              >  5:203016/1‑146 (MQ=255)
                         cATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCt                                                                                                                         >  4:182280/1‑149 (MQ=255)
                          aTCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGACAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTg                                                                                                                        >  1:98140/1‑149 (MQ=255)
                              ttGTGTGCCCAATTAGTTGGGGGTTTACCCGCGAATACGGTAGCGGCGGCAAGGGTGGATTACGGGTGGATATTGGCCAGCACATACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGaa                                                                                                                     <  6:343786/148‑1 (MQ=38)
                                 tgtgCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAgg                                                                                                                   >  5:184559/1‑147 (MQ=255)
                                 tgtcCCCCCATAGGAGCGTGTTTACCCGCCAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGc                                                                                                                  <  6:360785/148‑1 (MQ=255)
                                    gCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCGTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCCGAt                                                                                                              >  5:32559/1‑149 (MQ=255)
                                       cacaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAACATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCa                                                                                                            <  1:298509/149‑1 (MQ=255)
                                       cacaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAg                                                                                                           <  6:184559/149‑1 (MQ=255)
                                        acaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAgg                                                                                                                   >  4:92232/1‑140 (MQ=255)
                                        acaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAgg                                                                                                                   <  3:92232/140‑1 (MQ=255)
                                        acaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAgg                                                                                                          >  5:370155/1‑149 (MQ=255)
                                        acaTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAgg                                                                                                          <  1:7809/149‑1 (MQ=255)
                                          aTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCCGATCAGCCTGGAAGGCAGATCAGGCg                                                                                                        >  2:292608/1‑149 (MQ=255)
                                             ggAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAAc                                                                                                     <  6:362786/149‑1 (MQ=255)
                                              gAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACg                                                                                                    <  7:440/149‑1 (MQ=255)
                                              gAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGAGAAGGGCCGCGGCGtt                                                                                                                                   >  6:342294/1‑117 (MQ=37)
                                              gAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTGAGAAGGGCCGCGGCGtt                                                                                                                                   <  5:342294/118‑2 (MQ=37)
                                               cgCGGGTAAACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACgg                                                                                                   >  8:148123/2‑149 (MQ=255)
                                               aGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACgg                                                                                                   <  3:182280/149‑1 (MQ=255)
                                                 cGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGaa                                                                                                                     <  5:131739/129‑1 (MQ=255)
                                                 cGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGaa                                                                                                                     >  6:131739/1‑129 (MQ=255)
                                                    gTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGAc                                                                                              >  2:38674/1‑149 (MQ=255)
                                                       tACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACgc                                                                                            <  5:217733/148‑1 (MQ=255)
                                                         cccGCTAATACGTTAGCGGCGGAAACGGTGTATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACgcag                                                                                          <  8:359491/148‑1 (MQ=255)
                                                               aattgggTAGCGGCGGCAAGGGTGTTTTAGGGGTGTATATTGGCCAGCACACACATGCCAACAATTCAAACCGCCGGGGCCCTTCTCAGGGCGGCGGCGTCAGATCAGCGTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGAt                                                                                    <  6:35644/143‑1 (MQ=255)
                                                                   cGGTAGCGGCGGCAACGGTGGATTACCGGTGGATTTTGCCCAGCACAGACATCCCCACAAATGAAACCGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGtct                                                                               <  7:328598/149‑1 (MQ=255)
                                                                   cGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGtc                                                                                <  6:368577/148‑1 (MQ=255)
                                                                    ggTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACgcagc                                                                                         <  2:4295/138‑1 (MQ=255)
                                                                    ggTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACgcagc                                                                                         >  1:4295/1‑138 (MQ=255)
                                                                            cggcAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGt                                                                      >  4:123494/1‑149 (MQ=255)
                                                                               cAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCa                                                                                                                 <  7:30346/103‑1 (MQ=255)
                                                                               cAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTATCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCa                                                                                                                 >  8:30346/1‑103 (MQ=255)
                                                                                        aTTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCg                                                          <  1:38674/149‑1 (MQ=255)
                                                                                            ccGGTGGATATTGGCCAGCACAGCCTTGCCCAAAAATGAAACCGCCGCGGCCCTTCTCTGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGaaa                                                       <  6:32559/148‑1 (MQ=255)
                                                                                             cGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGcc                                                               >  8:147477/1‑139 (MQ=255)
                                                                                             cGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGcc                                                               <  7:147477/139‑1 (MQ=255)
                                                                                             cGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAat                                                      >  7:61208/1‑148 (MQ=255)
                                                                                             cGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAat                                                      >  6:226786/1‑148 (MQ=255)
                                                                                               gTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatga                                                    >  2:309355/1‑148 (MQ=255)
                                                                                               gTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAatga                                                    <  1:309355/148‑1 (MQ=255)
                                                                                                  gATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCGCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGt                                                 <  6:370155/148‑1 (MQ=255)
                                                                                                  catatTGGCCAGCACAGACATGCCCACAAAGGAAAACGCCGCGGCCGTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTGGCGGTCCGGGCCGGTCGAAATGATGt                                                 <  1:203533/147‑1 (MQ=255)
                                                                                                      ttGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGATCGACATCATTTCGac                                             <  2:277907/148‑2 (MQ=255)
                                                                                                          ccAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGa                                                                                     >  4:313101/1‑104 (MQ=255)
                                                                                                          ccAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGa                                                                                     <  3:313101/104‑1 (MQ=255)
                                                                                                           cAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGgcgc                                       >  5:177808/1‑149 (MQ=255)
                                                                                                                        cccACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGTCGGTCGAAATGATGt                                                 <  6:313310/126‑1 (MQ=255)
                                                                                                                        cccACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGTCGGTCGAAATGATGt                                                 >  5:313310/1‑126 (MQ=255)
                                                                                                                        cccACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGctct                          >  4:275871/1‑149 (MQ=255)
                                                                                                                        cccACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGctct                          >  4:119240/1‑149 (MQ=255)
                                                                                                                           acaAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAAAGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACAAGCTCTTc                        >  4:325243/1‑148 (MQ=255)
                                                                                                                             aaaTGAAAACGCCGCGGCCCGTCTCACGGCCGCGGCGTCAGATAACCCAGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGAt                     >  4:53513/1‑149 (MQ=255)
                                                                                                                             aaaTCAAAACGCCGCGGCCCTTCTCCGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGAt                     <  7:238209/149‑1 (MQ=255)
                                                                                                                              aaTGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATg                    <  3:119240/149‑1 (MQ=255)
                                                                                                                                   aaaaCGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCt                >  1:212326/1‑149 (MQ=255)
                                                                                                                                      ccgccgcGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGa             <  4:179918/148‑1 (MQ=255)
                                                                                                                                      aCGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGa             >  4:250924/1‑149 (MQ=255)
                                                                                                                                          cgcgGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTa         <  5:226786/149‑1 (MQ=255)
                                                                                                                                          cgcgGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTa         <  1:292608/149‑1 (MQ=255)
                                                                                                                                          cgcgGCCCTGAGAAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTa         >  3:90773/1‑149 (MQ=255)
                                                                                                                                                 cTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGaa                                                        <  6:321660/95‑1 (MQ=255)
                                                                                                                                                 cTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGaa                                                        >  5:321660/1‑95 (MQ=255)
                                                                                                                                                 cTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCAcgc  >  8:44516/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GCTCCTACAACGGCCTGCGCCCAACCATCGTTGTGTGCCCACATAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGGATTACCGGTGGATATTGGCCAGCACAGACATGCCCACAAATGAAAACGCCGCGGCCCTTCTCAGGGCCGCGGCGTCAGATCAGCCTGGAAGGCAGATCAGGCGAACGGGTGACGCAGCACGATGGTCTCGTTGCGGTCCGGGCCGGTCGAAATGATGTCGATCGGCGCGCCAACCAGCTCTTCGATGCGCTTGATGTAGTCACGC  >  NC_002947/5557776‑5558068

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