Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,291,245 G→C 5.4% P281P (CCG→CCC PP_4664 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,291,2450GC 5.4% 232.0 / 2.1 73P281P (CCG→CCCPP_4664LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (33/36);  new base C (2/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

TGGTCGCGGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCC  >  NC_002947/5291097‑5291384
                                                                                                                                                    |                                                                                                                                            
tGGTCGCGGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCg                                                                                                                                              <  1:255120/149‑1 (MQ=255)
gggTCGCGGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCg                                                                                                                                              <  7:221128/148‑1 (MQ=255)
  gTCGCGGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCggg                                                                                                                                            <  5:421901/149‑1 (MQ=255)
  gTCGCGGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCgg                                                                                                                                             <  4:177815/148‑1 (MQ=255)
   tCGCGGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGc                                                                                                                                           >  2:134536/1‑149 (MQ=255)
    cgcgGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGc                                                                                                                                           >  7:254473/1‑148 (MQ=255)
      cgggcgggGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGCGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTAt                                                                                                                                        <  6:425556/149‑1 (MQ=255)
         gcgggGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGcc                                                                                                                                     <  8:254473/149‑1 (MQ=255)
             ggCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAg                                                                                                                                 <  5:108902/149‑1 (MQ=255)
              gCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAg                                                                                                                                 >  3:190442/1‑148 (MQ=255)
               ccGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGaa                                                                                                                               >  1:27166/1‑149 (MQ=255)
                cGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAAc                                                                                                                              <  2:203270/149‑1 (MQ=255)
                cGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAAc                                                                                                                              >  1:186016/1‑149 (MQ=255)
                cGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAAc                                                                                                                              <  8:114583/149‑1 (MQ=255)
                 ggAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACt                                                                                                                             <  7:325160/149‑1 (MQ=255)
                  gAATCGGGGTGGTGCCACTGTCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTg                                                                                                                            >  5:342278/1‑149 (MQ=255)
                       ggggTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGCCGCaaaatcc                                                                                                                            <  3:89456/144‑4 (MQ=255)
                       ggggTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGCCGCaaaatcc                                                                                                                            >  4:89456/1‑141 (MQ=255)
                         ggtggtGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCa                                                                                                                     <  4:190442/149‑1 (MQ=255)
                          gtggtgCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCAc                                                                                                                    <  4:199429/149‑1 (MQ=255)
                              tgCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCAc                                                                                                                    <  2:102849/145‑1 (MQ=255)
                              tgCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCAc                                                                                                                    >  1:102849/1‑145 (MQ=255)
                               gCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCa                                                                                                                     >  4:287355/1‑143 (MQ=255)
                               gCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCa                                                                                                                     <  3:287355/143‑1 (MQ=255)
                                 cACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTagc                                                                                                             <  2:186016/149‑3 (MQ=255)
                                 cACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTAcc                                                                                                             <  1:134536/149‑1 (MQ=255)
                                  aCTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCg                                                                                                            >  3:222646/1‑149 (MQ=255)
                                      gcagcaGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGc                                                                                                          <  1:244410/147‑1 (MQ=255)
                                      gcagcaGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCt                                                                                                         <  8:374735/148‑1 (MQ=255)
                                         gcagTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCGACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCAGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGAcg                                                                                                     >  7:57485/1‑149 (MQ=255)
                                          cagTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGAcgc                                                                                                    <  2:27166/149‑1 (MQ=255)
                                             tACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACgcgc                                                                                                  >  8:100786/1‑148 (MQ=255)
                                                aGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAgg                                                                                              >  5:299783/1‑149 (MQ=255)
                                                  cgcTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGgc                                                                                            >  5:411986/1‑149 (MQ=255)
                                                    cTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGgcgc                                                                                          <  1:383083/149‑1 (MQ=255)
                                                       gCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGcaca                                                                                       >  2:273912/1‑149 (MQ=255)
                                                       gCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGCCg                                                                                                                                    >  4:109743/1‑104 (MQ=255)
                                                       gCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATAGCCCGGCAATGCCg                                                                                                                                    <  3:109743/104‑1 (MQ=255)
                                                          agtcATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCGATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCAcacca                                                                                    <  6:342278/146‑1 (MQ=255)
                                                              aTGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCa                                                                                                                                  <  3:233005/99‑1 (MQ=255)
                                                              aTGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCa                                                                                                                                  >  4:233005/1‑99 (MQ=255)
                                                                gctgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTc                                                                              >  4:209810/1‑148 (MQ=255)
                                                                  tgcCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAg                                                                           >  8:328741/1‑149 (MQ=255)
                                                                    ccGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAg                                                                           <  7:100786/147‑1 (MQ=255)
                                                                      gTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGTGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGa                                                                            >  1:248206/1‑144 (MQ=255)
                                                                      gTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGTGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGa                                                                            <  2:248206/144‑1 (MQ=255)
                                                                      gTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTa                                                                                                                         <  1:562/100‑1 (MQ=255)
                                                                      gTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTa                                                                                                                         >  2:562/1‑100 (MQ=255)
                                                                        tGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAAc                                                                     >  5:258630/1‑149 (MQ=255)
                                                                         gCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAAcc                                                                    >  3:3061/1‑149 (MQ=255)
                                                                               ggTGGCACTGCAGGAGCCATGGGCCGACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCt                                                                <  8:57485/147‑1 (MQ=255)
                                                                                 tGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCa                                                            >  8:264288/1‑149 (MQ=255)
                                                                                    cACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAg                                                         >  8:390411/1‑149 (MQ=255)
                                                                                         cAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGcgcg                                                     <  4:325928/148‑1 (MQ=255)
                                                                                         cAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGcgcg                                                     >  3:325928/1‑148 (MQ=255)
                                                                                           ggAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGc                                                  <  8:304312/149‑1 (MQ=255)
                                                                                           ggAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGc                                                  <  7:390411/149‑1 (MQ=255)
                                                                                                 atttgCCAACCGACGTTTGCTGGTGAGCGCACGGTATTTTGCGGCATTTCCGGGCTATGCCGAAGAACTGGTACCCACCTTTACAGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGTGCAGacac                                            >  2:297638/5‑149 (MQ=255)
                                                                                                  tGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGacac                                            <  1:273912/149‑1 (MQ=255)
                                                                                                  tGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACAcc                                           >  2:47254/1‑149 (MQ=255)
                                                                                                    ggCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCt                                          <  5:294189/148‑1 (MQ=255)
                                                                                                     gCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTg                                         >  8:282813/1‑148 (MQ=255)
                                                                                                      ccAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGa                                       >  2:168944/1‑149 (MQ=255)
                                                                                                       cAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGAc                                      >  5:233145/1‑149 (MQ=255)
                                                                                                           cgccgGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCg                                  <  7:282813/149‑1 (MQ=255)
                                                                                                                gTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGccc                             <  7:264288/149‑1 (MQ=255)
                                                                                                                   gCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAg                           <  8:55016/148‑1 (MQ=255)
                                                                                                                        tgtgCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACt                     >  1:38590/1‑149 (MQ=255)
                                                                                                                         gtgCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTg                    <  3:209810/149‑1 (MQ=255)
                                                                                                                           gcgcACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGc                                                              >  6:422343/1‑105 (MQ=255)
                                                                                                                           gcgcACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGc                                                              <  5:422343/105‑1 (MQ=255)
                                                                                                                            cgcACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgcgg                 <  4:222646/149‑1 (MQ=255)
                                                                                                                            cgcACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgcg                  <  7:328741/148‑1 (MQ=255)
                                                                                                                               aCGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgcggcgc              >  2:45138/1‑149 (MQ=255)
                                                                                                                                 gggATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTgcggcgcgg            <  6:258630/149‑1 (MQ=255)
                                                                                                                                         gCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTgg     <  1:45138/148‑1 (MQ=255)
                                                                                                                                           ggCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGcc  >  5:60842/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
TGGTCGCGGGCGGGGCCGGAATCGGGGTGGTGCCACTGGCAGCAGTACAGCGCTGGCAAGGCATGCTGCCGTTGCACTGGGTGGCACTGCAGGAGCCATGGGCCAACCGCCGGTTGCTGGTGTGCGCACGGGATTTTGCGGCATTGCCGGGCTATGCCGCAGAACTGGTAGCCACCTTTACCGGCTGACGCGCAAGGGCGCACACCACTGTCGAGGCAAACCGCCTGCCATAGCGCGGGCAGACACCTGGACGTCGAGCCCAGGCAACTGCGGCGCGGCCGTTGGGCC  >  NC_002947/5291097‑5291384

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: