Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,612,253 C→G 5.9% intergenic (+1215/‑892) hcpC‑II → / → PP_4084 hemolysin‑coregulated protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,612,2530CG 5.9% 215.0 / 3.4 68intergenic (+1215/‑892)hcpC‑II/PP_4084hemolysin‑coregulated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/30);  new base G (2/2);  total (36/32)
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.71e-01

GACCTCAATGCGCATTGGCGGGTCAACACGGTTACCCATAAAGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGACCGCGAAAGCCAGAACAACGAGTTCAGT  >  NC_002947/4612108‑4612398
                                                                                                                                                 |                                                                                                                                                 
gacCTCAATGCGCATTGGCGGGTCAACACGGTTACCCATAAAGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTgc                                                                                                                                                >  2:263270/1‑149 (MQ=255)
   cTCAATGCGCATTGGCGGGTCAACACGGTTACCCATAAAGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTgcgc                                                                                                                                              >  5:218640/1‑148 (MQ=255)
     cAATGCGCATTGGCGGGTCAACACGGTTACCCATAAAGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCAt                                                                                                                                           >  1:228307/1‑149 (MQ=255)
     cAATGCGCATTGGCGGGTCAACACGGTTACCCATAAAGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCAGGGCGCCAt                                                                                                                                           >  5:9831/1‑149 (MQ=255)
       aTGCGCATTGGCGGGTCAACACGGTTACCCATAAAGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCAGGGCGCCAttc                                                                                                                                         <  6:9831/149‑2 (MQ=255)
          cgcATTGGCGGTTCAACACGGTTACCCATAAATGCCACCAGTTCACCTGCCTGCAGTAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGAGGGCAGTGCTGGTCCCCGTTGGGAGCGATTGGCGCCGTGCGCCATTGaaa                                                                                                                                      <  4:315900/149‑1 (MQ=255)
          cgcATTGGCGGGTCAACACGGTTACCCATAAAGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGaa                                                                                                                                       <  5:253542/148‑1 (MQ=255)
                   gggTCAACACGGTTACCCATAAAGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCTCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGa                                                                                                                             <  1:69449/149‑1 (MQ=255)
                   gggTCAACACGGTTACCCATAAAGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACg                                                                                                                              >  5:161377/1‑148 (MQ=255)
                   gggTCAACACGGTTACCCATAAAGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGATGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCAc                                                                                                                               <  4:75563/147‑1 (MQ=255)
                   gggTCAACAAGGTTACCCATAAAGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCTGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCAc                                                                                                                               >  3:75563/1‑147 (MQ=255)
                        aaCACGGTTACCCATAAAGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGa                                                                                                                        >  5:298014/1‑149 (MQ=255)
                                   ccATAAAGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGAAGTGCTGGTCCCCGGTCGGACCGAATGCCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGcaca                                                                                                             >  8:155592/1‑149 (MQ=255)
                                   ccATAAAGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACAAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGcaca                                                                                                             >  5:14197/1‑149 (MQ=255)
                                         aGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCa                                                                                                       <  2:220528/149‑1 (MQ=255)
                                          ggCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCAc                                                                                                      <  6:298014/149‑1 (MQ=255)
                                              accaGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACcgtc                                                                                                  <  6:218640/149‑1 (MQ=255)
                                                   ttCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCgg                                                                                             <  6:175528/149‑1 (MQ=255)
                                                      aCCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGccc                                                                                          >  4:93169/1‑149 (MQ=255)
                                                      aCCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACAGTCCTCGGccc                                                                                          >  4:240493/1‑149 (MQ=255)
                                                        cAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGAc                                                                                                                       >  3:102636/1‑118 (MQ=255)
                                                        cAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGAc                                                                                                                       <  4:102636/118‑1 (MQ=255)
                                                           ccTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTg                                                                                     >  3:136666/1‑149 (MQ=255)
                                                            cTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGCACCGCCTgg                                                                                    >  4:25148/1‑149 (MQ=255)
                                                              gCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGGCGTCGGCCCGCCTGGCg                                                                                  >  7:295775/1‑149 (MQ=255)
                                                                aGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGTGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAg                                                                                >  4:106711/1‑149 (MQ=255)
                                                                  gaagaaGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACg                                                                                                                      <  7:307300/109‑1 (MQ=255)
                                                                  gaagaaGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACg                                                                                                                      >  8:307300/1‑109 (MQ=255)
                                                                   aagaagCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGa                                                                              >  4:73621/1‑148 (MQ=255)
                                                                    agaagCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAAGCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGaaa                                                                            <  1:263270/149‑1 (MQ=255)
                                                                     gaagCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAAt                                                                           <  6:305338/149‑1 (MQ=255)
                                                                      aagCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGTGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATc                                                                          <  3:175733/149‑1 (MQ=255)
                                                                        gccgccGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTTAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTa                                                                        >  6:298586/1‑149 (MQ=255)
                                                                        gccgccGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTa                                                                        <  3:93169/149‑1 (MQ=255)
                                                                          cgccgGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACt                                                                      <  7:155592/149‑1 (MQ=255)
                                                                                   cGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAAt                                                             <  2:228307/149‑1 (MQ=255)
                                                                                          aGCACGCGCTACAAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGcc                                                       <  6:14197/148‑1 (MQ=255)
                                                                                           gCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCgg                                                     <  6:161377/149‑1 (MQ=255)
                                                                                                 gcTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTaaaa                                               >  2:146155/1‑149 (MQ=255)
                                                                                                   tACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAAt                                                             >  4:20221/1‑133 (MQ=255)
                                                                                                   tACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAAt                                                             <  3:20221/133‑1 (MQ=255)
                                                                                                   tACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGt                                             >  6:196978/1‑149 (MQ=255)
                                                                                                     cGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCa                                           <  8:108532/149‑1 (MQ=255)
                                                                                                        gCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGtt                                        >  1:176833/1‑149 (MQ=255)
                                                                                                           gACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTcc                                     >  3:189805/1‑149 (MQ=255)
                                                                                                           gACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTcc                                     >  1:226860/1‑149 (MQ=255)
                                                                                                             ccGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTcccc                                   <  4:189805/149‑1 (MQ=255)
                                                                                                              cGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCt                                  >  5:138843/1‑149 (MQ=255)
                                                                                                              cGCAGTGCTGGTCCCCGGGCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGACGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATTGGGCCGGGTAAAAGTCAGTTTCCCCt                                  <  6:171450/149‑1 (MQ=255)
                                                                                                                cAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTgg                                <  3:73621/149‑1 (MQ=255)
                                                                                                                  gTGCTGGTCCCCGGTCGGACCGTATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGTCCGCACATGGCCAGCGTCGTCGGCCCGCCTGGCGATGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGa                              <  3:240493/149‑1 (MQ=255)
                                                                                                                   tGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGAc                             >  7:26001/1‑149 (MQ=255)
                                                                                                                       ggTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATAGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTAGTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGACcgc                          >  1:5531/1‑148 (MQ=255)
                                                                                                                           cccGGTCGGACCGAATGGCGCCCTGCGCCATTTAAAAAAACACGAATCGACGGGCCGCACATGGCCACCCTCGTCGGCCCGCCAGGCGAGGAAAACAACAGCGATGAATGGGGCCGGGTAAAAGTCAATTTACCAGGGGCCCGCGAAAg                     >  2:265521/1‑149 (MQ=18)
                                                                                                                           cccGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGACCGCGaaa                      <  6:138843/148‑1 (MQ=255)
                                                                                                                            ccGGTCGGACCGAATGGCGCCCTGTGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGTGTAAAAGTCAGTTTCCCCTGGGACCGCGAAAGc                    <  3:106711/149‑1 (MQ=255)
                                                                                                                              ggTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGACCGCGAAAGCCa                  <  8:26001/149‑1 (MQ=255)
                                                                                                                                   gACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGACCGCGAAAGCCAGaac              >  2:331010/1‑148 (MQ=255)
                                                                                                                                   gACCGAATGCCGCCCTGCGCAATTGAAAAAACCACGAATCGACGGGCCGCAAATGGCCACCGTCGTCGGCCCCCCTGGCGAGTAAATCTACACGGATGAAAGGGAACGGGTAAAAGTCAGTTACACATGGGACCCCCAAAGCCAGaac              >  6:215251/1‑148 (MQ=18)
                                                                                                                                    aCCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGACCGCGAAAGCCAGagcaa            >  3:264686/1‑149 (MQ=255)
                                                                                                                                    aCCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGACCGCGAAAGCCAGaacac            >  5:250254/1‑148 (MQ=255)
                                                                                                                                    aCCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGACCGCGAAAGCCAGaacaa            >  1:262576/1‑149 (MQ=255)
                                                                                                                                       gAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGACCGCGAAAGCCAGAACAACGa         <  5:196978/149‑1 (MQ=255)
                                                                                                                                        aaTGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGACCGCGAAAGCCAGAACAACGag        <  4:136666/149‑1 (MQ=255)
                                                                                                                                        aaTGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGACCGCGAAAGCCAGAACAACGag        >  6:231965/1‑149 (MQ=255)
                                                                                                                                        aaTGGCGCAGGGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGACCGCGAAAGCCAGAACAACGag        >  6:165528/1‑149 (MQ=255)
                                                                                                                                             cgcCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGACCGCGAAAGCCAGAACAACGagttcag   <  1:146155/149‑1 (MQ=255)
                                                                                                                                              gcCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGACCGCGAAAGCCAGAACAACGagttcagt  <  2:176833/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GACCTCAATGCGCATTGGCGGGTCAACACGGTTACCCATAAAGGCCACCAGTTCACCAGCCTGCAGGAAGAAGCCGCCGGCGCCGATGTAAGCACGCGCTACGAGCAGACCGCAGTGCTGGTCCCCGGTCGGACCGAATGGCGCCCTGCGCCATTGAAAAAACCACGAATCGACGGGCCGCACATGGCCACCGTCGTCGGCCCGCCTGGCGAGGAAATCTACTGCGATGAATGGGGCCGGGTAAAAGTCAGTTTCCCCTGGGACCGCGAAAGCCAGAACAACGAGTTCAGT  >  NC_002947/4612108‑4612398

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