Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,061,026 T→C 5.7% intergenic (‑107/+42) 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,0260TC 5.7% 258.1 / 3.9 89intergenic (‑107/+42)PP_5312/hupAhypothetical protein/DNA‑binding protein HU‑alpha
Reads supporting (aligned to +/- strand):  ref base T (44/39);  new base C (2/3);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.54e-01

AGCCCCATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTCC  >  NC_002947/6060885‑6061173
                                                                                                                                             |                                                                                                                                                    
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)
   cccATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCtt                                                                                                                                                      >  8:282747/1‑139 (MQ=255)
   cccATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCtt                                                                                                                                                      <  7:282747/139‑1 (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                                                           >  1:133637/1‑135 (MQ=255)
                                                                                                 gCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGc                                                           <  2:133637/135‑1 (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   <  2:378713/149‑1 (MQ=255)
                                                                                                                                           cTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTc   >  7:357805/1‑149 (MQ=255)
                                                                                                                                            ttGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTcc  >  2:201830/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
AGCCCCATGGCCCAGAGAAGAAAGCGGAACTTCATCAACGCGCCTCGAATGAAAAGTGACTAAGGCGCGCAGTTTATCCTTTTCCAGAACTATTTATGCAAGTTCGAAACATTTTGCCTGCAGAAAGCACAACGGGCGCCTTGGGGCGCCCGCTTGGCCGACACTGGGGTTGTCAGAGAGCGATCACTTGCCTTTCTTGGCTGGCTTGGCGGGGGCGTTAACGCTGTCACGCAGGTTCTTGCCAGGCTTGAAGGCCACGGTATTGCTGGCCTTGATCTTGACCGGTTCC  >  NC_002947/6060885‑6061173

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: