Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,061,027 G→A 5.6% intergenic (‑108/+41) PP_5312 ← / ← hupA hypothetical protein/DNA‑binding protein HU‑alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,061,0270GA 5.6% 264.0 / 2.8 88intergenic (‑108/+41)PP_5312/hupAhypothetical protein/DNA‑binding protein HU‑alpha
Reads supporting (aligned to +/- strand):  ref base G (44/39);  new base A (2/3);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.97e-01

AGCCCCATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTCCC  >  NC_002947/6060885‑6061174
                                                                                                                                              |                                                                                                                                                    
aGCCCCATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGgcgc                                                                                                                                                >  3:347673/1‑149 (MQ=255)
aGCCCCATGGCAAAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGAGCAGTTTATCCTTTTCCAGAACTATTTATGAAAGTTCGAAAGATTATGCCTGCAGAAAGCACAACGGCCGCCTTGGGgcgc                                                                                                                                                >  4:19271/1‑149 (MQ=255)
  ccccATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGccc                                                                                                                                              <  1:161982/149‑1 (MQ=255)
  ccccATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGccc                                                                                                                                              >  3:257192/1‑149 (MQ=255)
      aTGGCCCAGAGAAGAACGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCCCAAGGCGCCCgttg                                                                                                                                          <  1:271253/149‑4 (MQ=255)
       tGGCCCAGAGAAGAAACCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTg                                                                                                                                         <  1:169277/149‑1 (MQ=255)
       ccggtgagagAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTg                                                                                                                                         <  8:289312/143‑1 (MQ=255)
            cAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACGAAGGCGCGCAGTTTAACCGTTTCCAGAACTAGTTATGCAAGTTCGAAACATTTTGCCTGCAGCACGCACAACGGGCGCCATGGGGCGCCCGCTTGGCCGa                                                                                                                                    >  7:327880/1‑149 (MQ=255)
             agagAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGa                                                                                                                                    >  1:183540/1‑148 (MQ=255)
              gagaAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGaca                                                                                                                                  >  5:148650/1‑149 (MQ=255)
              gagaAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCCCAAGGCGCCCgttgtgc                                                                                                                                       >  5:301602/1‑138 (MQ=255)
              gagaAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCCCAAGGCGCCCgttgtgc                                                                                                                                       <  6:301602/144‑7 (MQ=255)
                 aagaaAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTg                                                                                                                               >  8:292527/1‑149 (MQ=255)
                 aacaacgCGGAACTTCATCAACGCGCCTCGAATGAAACGTCACTAAGGCGCGCAGTTTATCCTTTTCCACAACTATTTATGCAAGTTCTAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTg                                                                                                                               <  4:165516/143‑1 (MQ=255)
                    aaaGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTgggg                                                                                                                            >  4:252456/1‑149 (MQ=255)
                     cctatctaCTTCATCAACTCGCGTCTAATGAAATGTGAATAAGGCGCTCTGTTTATCCTTTTCCATAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACACCGGGCGCCTGGGGGCGCCCGCTTGGCCGACACTgggg                                                                                                                            <  4:225417/141‑1 (MQ=255)
                         ggAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTc                                                                                                                       <  2:253161/149‑1 (MQ=255)
                          taaCTACCTCAACGCGCCTCTAATGAAATTTGACTTAGGCGCGCAGTTTATCCTTTTCCATAACTATTTATGCATGTTCGAAACATTTTGCCTGCACAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCa                                                                                                                      <  4:347673/148‑1 (MQ=255)
                            aCTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATTCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCaga                                                                                                                    <  7:140566/149‑1 (MQ=255)
                            aCTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCaga                                                                                                                    <  3:119524/149‑1 (MQ=255)
                             cTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTgggg                                                                                                                            >  3:225417/1‑140 (MQ=255)
                              ttCATCAACGCGCCTCTATTTAAAATTTAATAAGCCGGGCATTTTATCCTTTACCAGAAATTTTTTTGCAAGTTCTAAAAATTTTGCGTGCAGAAAGCAAAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCagaga                                                                                                                  <  4:11004/149‑1 (MQ=255)
                              ttCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCagaga                                                                                                                  <  3:252456/149‑1 (MQ=255)
                               tattcaACGCGCCGCGAAGGAAATGGTACTACGGCGCGCAGTTTATCCTTTTCCAGAACAATTTTTGCAAGTTCGAAACTTTTGGCCTGCCGAAAGCACACCGGGCGCGTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCagaga                                                                                                                  <  8:327880/145‑1 (MQ=255)
                                  tcaACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGa                                                                                                              >  7:94591/1‑149 (MQ=255)
                                   caACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGa                                                                                                              >  6:49194/1‑148 (MQ=255)
                                    aaCGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTggg                                                                                                                                                    >  4:33720/1‑109 (MQ=255)
                                    aaCGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTggg                                                                                                                                                    <  3:33720/109‑1 (MQ=255)
                                    aaCGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGAt                                                                                                             >  6:56608/1‑148 (MQ=255)
                                     aCGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCa                                                                                                           >  6:335678/1‑149 (MQ=255)
                                     aCGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGCGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATc                                                                                                            >  7:280716/1‑148 (MQ=255)
                                         gcCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTg                                                                                                      <  7:292527/149‑1 (MQ=255)
                                         gcCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCCCAAGGCGCCCgtt                                                                                                                                           <  7:165797/113‑3 (MQ=255)
                                         gcCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCCCAAGGCGCCCgtt                                                                                                                                           >  8:165797/1‑111 (MQ=255)
                                           cTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGcc                                                                                                    >  6:152594/1‑149 (MQ=255)
                                              gAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCttt                                                                                                 >  4:54739/1‑149 (MQ=255)
                                               aaTGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTc                                                                                                >  3:267014/1‑149 (MQ=255)
                                                aTGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCGTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCt                                                                                               <  5:355334/149‑1 (MQ=255)
                                                    aaaGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTTTTTATCCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTtggc                                                                                           <  2:199724/149‑1 (MQ=255)
                                                    aaaGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTtggc                                                                                           <  5:49194/149‑1 (MQ=255)
                                                           cTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTg                                                                                    >  5:120159/1‑149 (MQ=255)
                                                           cTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTg                                                                                    >  4:291351/1‑149 (MQ=255)
                                                              aGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCg                                                                                 <  7:92156/149‑1 (MQ=255)
                                                              aGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAAAATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACAATGGGGTAGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCg                                                                                 >  7:296613/1‑149 (MQ=255)
                                                                gcgcgcAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGCCTTGGCggt                                                                               >  2:262011/1‑148 (MQ=255)
                                                                gcgcgcAGTTTATCCCTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCggg                                                                               >  7:44429/1‑149 (MQ=255)
                                                                   cgcAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTc                                                                                                >  1:188576/1‑129 (MQ=255)
                                                                   cgcAGTTTATCCTTATCCCGAACTCTTTATGCAAGTGCGAAACATTTTGCCTGCCGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGc                                                                            <  4:257192/149‑1 (MQ=255)
                                                                       gTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGtt                                                                         >  3:325925/1‑148 (MQ=255)
                                                                         ttATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAAc                                                                      >  6:352024/1‑149 (MQ=255)
                                                                            tCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGc                                                                    <  5:352024/148‑1 (MQ=255)
                                                                            tCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCt                                                                   <  8:94591/149‑1 (MQ=255)
                                                                            tCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCt                                                                   <  5:335678/149‑1 (MQ=255)
                                                                              cTTTTCCAGAACTATTTAGGCAAGTTCGAATCATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGt                                                                 <  5:178197/149‑1 (MQ=255)
                                                                               ttttCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCggg                                                                               <  7:236424/134‑1 (MQ=255)
                                                                               ttttCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCggg                                                                               >  8:236424/1‑134 (MQ=255)
                                                                               ttttCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGt                                                                 >  2:166124/1‑148 (MQ=255)
                                                                               ttttCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGt                                                                 >  1:378713/1‑148 (MQ=255)
                                                                                     agagCTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACtt                                                                                                        <  1:63049/104‑1 (MQ=255)
                                                                                     aGATCTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACtt                                                                                                        >  2:63049/1‑104 (MQ=255)
                                                                                      gAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTCGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAg                                                          <  1:79199/149‑1 (MQ=255)
                                                                                          tATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCt                                                     >  1:217035/1‑149 (MQ=255)
                                                                                               aTGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGc                                                  >  3:369316/1‑147 (MQ=255)
                                                                                                tGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCa                                                <  6:120159/148‑1 (MQ=255)
                                                                                                 gCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGc                                                            <  2:133637/135‑1 (MQ=255)
                                                                                                 gCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGc                                                            >  1:133637/1‑135 (MQ=255)
                                                                                                     gTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTg                                          <  3:291351/149‑1 (MQ=255)
                                                                                                         gAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAg                                       <  5:152594/148‑1 (MQ=255)
                                                                                                         gAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGCGCGCCCGCTTGGCCGACACTCGGGTTCTCAGAGACCGATCACTTGCCTTACTTGGCTGGCTTGGCGCGCGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAgg                                      >  8:358164/1‑149 (MQ=255)
                                                                                                               ttttGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACgg                                <  1:166124/149‑1 (MQ=255)
                                                                                                                  tGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTa                              <  8:44429/148‑1 (MQ=255)
                                                                                                                  tGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTa                              <  4:369316/148‑1 (MQ=255)
                                                                                                                   gcgtGCAGAAAGCACAACGGGCGCTTTGGGGCGCCCGCTTGGCCGACACTTGGGTTGTCAGAGAGCGATGACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACTCTTTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTAt                             <  8:296613/145‑1 (MQ=255)
                                                                                                                      tGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTtggctggc                                                                                       >  6:360071/1‑87 (MQ=255)
                                                                                                                        cAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTgg                       >  2:11859/1‑149 (MQ=255)
                                                                                                                          gAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGATTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGcc                     >  1:275546/1‑149 (MQ=255)
                                                                                                                           aaaGCACAACGGGCGCCTTGGGGGGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCt                    <  3:19271/149‑1 (MQ=255)
                                                                                                                           aaaGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCt                    <  3:173896/149‑1 (MQ=255)
                                                                                                                           aaaGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCt                    >  5:260297/1‑149 (MQ=255)
                                                                                                                               cacaACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGAt                >  6:192678/1‑149 (MQ=255)
                                                                                                                                 caACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCt              >  5:24125/1‑149 (MQ=255)
                                                                                                                                  aaCGGGGGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGc                                                                                                                <  8:79521/51‑1 (MQ=255)
                                                                                                                                  aaCGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGc                                                                                                                >  7:79521/1‑51 (MQ=255)
                                                                                                                                  aaCGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCtt             >  5:226518/1‑149 (MQ=255)
                                                                                                                                  aaCGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCtt             <  6:148650/149‑1 (MQ=255)
                                                                                                                                     gggCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGAc          <  3:54739/149‑1 (MQ=255)
                                                                                                                                        cgcCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCgg       >  3:114975/1‑149 (MQ=255)
                                                                                                                                           cTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTc    >  7:357805/1‑149 (MQ=255)
                                                                                                                                           cTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTc    <  2:378713/149‑1 (MQ=255)
                                                                                                                                            ttGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTcc   >  2:201830/1‑149 (MQ=255)
                                                                                                                                              ggggCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTccc  <  5:192678/148‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
AGCCCCATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTCCC  >  NC_002947/6060885‑6061174

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: