Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,038,399 C→G 5.4% A66A (GCC→GCG PP_3559 → glycine betaine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,038,3990CG 5.4% 241.1 / 4.1 74A66A (GCC→GCGPP_3559glycine betaine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (35/35);  new base G (2/2);  total (37/37)
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.62e-01

TGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTAC  >  NC_002947/4038260‑4038544
                                                                                                                                           |                                                                                                                                                 
tGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGc                                                                                                                                           >  3:336650/1‑148 (MQ=255)
 ggTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGc                                                                                                                                           >  4:353384/1‑147 (MQ=255)
  gTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGGCTGCTGGCCTGGCACGCCAg                                                                                                                                        <  5:14137/149‑1 (MQ=255)
   tCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCGTgccaggccag                                                                                                                                        >  3:56017/1‑148 (MQ=255)
   tCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCGTGCCAGGCCAGc                                                                                                                                       <  4:56017/149‑1 (MQ=255)
        tggttggtATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAg                                                                                                                                  >  5:224949/1‑149 (MQ=255)
          gttggtATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGcc                                                                                                                                >  8:191581/1‑149 (MQ=255)
          gttggtATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGcc                                                                                                                                >  7:314035/1‑149 (MQ=255)
           ttggtATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCct                                                                                                                               >  5:179191/1‑149 (MQ=255)
                 ttGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCtgctgc                                                                                                                          <  5:267293/148‑1 (MQ=255)
                   gAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTgcg                                                                                                                         >  4:250613/1‑147 (MQ=255)
                   gAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTgcg                                                                                                                         >  6:212352/1‑147 (MQ=255)
                           gCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTg                                                                                                               >  4:82618/1‑149 (MQ=255)
                              accaccTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTg                                                                                                            >  3:227462/1‑149 (MQ=255)
                                  ccTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCGTGCCAGGCCAGca                                                                                                                                      <  8:337982/119‑2 (MQ=255)
                                  ccTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCGTGCCAGGCCAGca                                                                                                                                      >  7:337982/1‑118 (MQ=255)
                                   cTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTcgc                                                                                                       >  3:66121/1‑149 (MQ=255)
                                    tGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTcgcg                                                                                                      >  4:290591/1‑149 (MQ=255)
                                    tGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTcgcg                                                                                                      >  8:165611/1‑149 (MQ=255)
                                      cgcAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTcgcg                                                                                                      <  6:224949/147‑1 (MQ=255)
                                        cAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCgctgc                                                                                                  <  1:252315/149‑1 (MQ=255)
                                           tttGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTgg                                                                                               >  2:164877/1‑149 (MQ=255)
                                            ttGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTgg                                                                                               <  3:250613/148‑1 (MQ=255)
                                              gTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGctggcctgg                                                                                                                                                <  7:288375/97‑1 (MQ=255)
                                              gTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGctggcctgg                                                                                                                                                >  8:288375/1‑97 (MQ=255)
                                              gTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCGGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCgctg                                                                                                   >  7:311281/1‑142 (MQ=255)
                                              gTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCGGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCgctg                                                                                                   <  8:311281/142‑1 (MQ=255)
                                              gTCTCACCACTTGATGAAACTTCTGGTCGGCCTTGAGAACCTGTTGCGCCTTTTTCCTTGGTGGCTTTTGTTGTTGCTGGTTGTTCTGCTGGCCTTTCACGCCAGCCGCAGCCTTCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGcc                                                                                             <  2:322474/148‑1 (MQ=255)
                                               tCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGccc                                                                                            <  3:82618/148‑1 (MQ=255)
                                                 tGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGa                                                                                         >  5:9302/1‑149 (MQ=255)
                                                   aCCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATc                                                                                       >  6:255577/1‑149 (MQ=255)
                                                     cAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCgg                                                                                     <  3:290591/149‑1 (MQ=255)
                                                      aGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGt                                                                                    <  3:353384/149‑1 (MQ=255)
                                                      aGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGt                                                                                    <  4:336650/149‑1 (MQ=255)
                                                       gatgctgcAACTGCTGGTCGGCCTGCAGCACCTGTGGCGCCTGTCGCCTTGGTGGCTGTCGTTGCTGCGGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTa                                                                                   <  8:314035/147‑1 (MQ=255)
                                                            tgcAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTc                                                                              >  3:340279/1‑149 (MQ=255)
                                                             gcAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTgg                                                                                               <  4:160054/131‑1 (MQ=255)
                                                             gcAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTgg                                                                                               >  3:160054/1‑131 (MQ=255)
                                                             gcAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCg                                                                             >  5:69811/1‑149 (MQ=255)
                                                              cAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTc                                                                              <  7:165611/147‑1 (MQ=255)
                                                                       gTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTgg                                                                     <  1:172907/147‑1 (MQ=255)
                                                                       gTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTgg                                                                     <  6:69811/147‑1 (MQ=255)
                                                                        tCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGa                                                                   >  1:174433/1‑148 (MQ=255)
                                                                          ggCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGcc                                                                                             <  6:67900/120‑1 (MQ=255)
                                                                          ggCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGcc                                                                                             >  5:67900/1‑120 (MQ=255)
                                                                          ggCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTGGCTGGTGTGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGTTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACaa                                                                <  3:325239/149‑1 (MQ=255)
                                                                            ccTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGt                                                              <  4:340279/149‑1 (MQ=255)
                                                                                gagaACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGt                                                                                    <  8:199796/123‑1 (MQ=255)
                                                                                gagaACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGt                                                                                    >  7:199796/1‑123 (MQ=255)
                                                                                  gaACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTtgctgc                                                        <  5:255577/149‑1 (MQ=255)
                                                                                   aaCCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTtgctg                                                         >  5:299337/1‑147 (MQ=255)
                                                                                     ccTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCAGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCaca                                                     >  3:199436/1‑147 (MQ=255)
                                                                                       tGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGAcc                                                   >  7:253456/1‑149 (MQ=255)
                                                                                             ggcTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGGCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGc                                                                                <  3:245247/112‑1 (MQ=255)
                                                                                             gcCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGAccctggcc                                             >  5:247543/1‑149 (MQ=255)
                                                                                             gcCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGAccctggcc                                             <  2:174433/149‑1 (MQ=255)
                                                                                             gcCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCCAGCTGCTGGCCCTGATCGGTAAGc                                                                                >  4:245247/1‑114 (MQ=255)
                                                                                                 gTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGAccctggccctgg                                         >  3:279059/1‑149 (MQ=255)
                                                                                                  ttGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGAccctggccctg                                          <  1:164877/147‑1 (MQ=255)
                                                                                                      cTTGGTGGCTGGGGTGGTTGCTGGTGGGGCGGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGGCGCGCGGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTg                                    <  4:66121/149‑1 (MQ=255)
                                                                                                          gtggCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCt                                                                                                             <  3:337917/72‑1 (MQ=255)
                                                                                                          gtggCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCt                                                                                                             >  4:337917/1‑72 (MQ=255)
                                                                                                          gtggCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTg                                    >  8:310775/1‑145 (MQ=255)
                                                                                                          gtggCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTg                                    <  7:310775/145‑1 (MQ=255)
                                                                                                          gtggCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCa                                <  4:279059/149‑1 (MQ=255)
                                                                                                          gtggCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGACa                                <  7:191581/149‑1 (MQ=255)
                                                                                                             gCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCTCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTa                             <  6:299337/149‑1 (MQ=255)
                                                                                                                   gttgttgCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTACTGACCAGCAccagg                        <  5:212352/148‑4 (MQ=255)
                                                                                                                    ttgttgCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGtt                       <  8:253456/148‑1 (MQ=255)
                                                                                                                      gttgCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTtgt                    >  7:308061/1‑149 (MQ=255)
                                                                                                                      gttgCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTtgt                    >  5:276747/1‑149 (MQ=255)
                                                                                                                       ttgCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGGATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTGAGTACCGGTTtgtg                   >  4:244361/1‑149 (MQ=255)
                                                                                                                             gtggGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTa               <  6:179191/147‑1 (MQ=255)
                                                                                                                                  tCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCGCTGTGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACtggtgg        <  4:227462/149‑1 (MQ=255)
                                                                                                                                   ctgctgGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCAGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTtgct                                                          <  2:356557/98‑1 (MQ=255)
                                                                                                                                   ctgctgGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCAGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTtgct                                                          >  1:356557/1‑98 (MQ=255)
                                                                                                                                         gcctggcACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTAc  >  4:275567/1‑148 (MQ=255)
                                                                                                                                           ctggcACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCgg                         <  2:43101/123‑1 (MQ=255)
                                                                                                                                           ctggcACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCgg                         >  1:43101/1‑123 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
TGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTAC  >  NC_002947/4038260‑4038544

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: