Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,103,207 G→A 6.6% G14S (GGC→AGC) ‡ PP_5353 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,103,2070GA 6.6% 271.0 / 5.7 90G14S (GGC→AGC) ‡PP_5353hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (45/39);  new base A (2/4);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.97e-01

AACTTTGTAACATTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGC  >  NC_002947/6103059‑6103352
                                                                                                                                                    |                                                                                                                                                 
aaCTTTGTAACATTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCg                                                                                                                                                   <  8:239718/149‑1 (MQ=255)
 aCTTTGTAACATTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCgg                                                                                                                                                  <  8:249185/149‑1 (MQ=255)
  cTTTGTAACATTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGc                                                                                                                                                 <  3:390393/149‑1 (MQ=255)
  cTTTGTAACATCCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGc                                                                                                                                                 <  7:305927/149‑1 (MQ=255)
   tttGTAACATTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGcc                                                                                                                                                >  1:298648/1‑149 (MQ=255)
        aaCATTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACg                                                                                                                                           >  2:302869/1‑149 (MQ=255)
           aTTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCt                                                                                                                                         >  5:210568/1‑148 (MQ=255)
           aTTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTg                                                                                                                                        <  1:137007/149‑1 (MQ=255)
            ttCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCt                                                                                                                                         >  5:25504/1‑147 (MQ=255)
             tCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAg                                                                                                                                      >  4:382907/1‑149 (MQ=255)
               aaGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGc                                                                                                                                                 >  3:218751/1‑136 (MQ=255)
               aaGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGc                                                                                                                                                 <  4:218751/136‑1 (MQ=255)
               aaGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGc                                                                                                                                     <  1:302869/148‑1 (MQ=255)
               aaGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCAGCGTCAGGCCGACGCAGAGcg                                                                                                                                    <  1:358822/149‑2 (MQ=255)
                aGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGc                                                                                                                                     <  7:380887/147‑1 (MQ=255)
                aGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGc                                                                                                                                     >  8:380887/1‑147 (MQ=255)
                aGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCAGGCCGACGCAGAGcgc                                                                                                                                   <  5:202318/149‑3 (MQ=255)
                aGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCAGCGTCAGGCCGACGCAGAGcga                                                                                                                                   >  5:167655/1‑147 (MQ=255)
                aGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCAGCGTCAGGCCGACGCAGAGc                                                                                                                                     <  6:167655/147‑1 (MQ=255)
                 gCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTg                                                                                                                                  >  8:197429/1‑149 (MQ=255)
                 gCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCAGCGTCAGGCCGACGCAGAGc                                                                                                                                     >  1:131517/1‑146 (MQ=255)
                 gCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCAGCGTCAGGCCGACGCAGAGc                                                                                                                                     <  2:131517/146‑1 (MQ=255)
                   aaaaGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGa                                                                                                                                >  8:220813/1‑149 (MQ=255)
                    aaaGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGAc                                                                                                                               <  3:382907/149‑1 (MQ=255)
                     aaGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACg                                                                                                                              >  8:194164/1‑149 (MQ=255)
                      aGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGc                                                                                                                             <  4:359108/149‑1 (MQ=255)
                           gAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAgcg                                                                                                                        >  2:331750/1‑149 (MQ=255)
                           gAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAgcg                                                                                                                        >  8:332051/1‑149 (MQ=255)
                           gAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAgcg                                                                                                                        >  4:388294/1‑149 (MQ=255)
                           gAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCAGGACGCCAgcg                                                                                                                        >  3:248582/1‑149 (MQ=255)
                            aaaCACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAgcgc                                                                                                                       >  3:376471/1‑149 (MQ=255)
                            aaaCACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAgcgc                                                                                                                       >  3:241139/1‑149 (MQ=255)
                              acacTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGcc                                                                                                                      >  6:285474/1‑148 (MQ=255)
                                  tATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGAGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCCGACGCTGa                                                                                                                                       >  7:284969/1‑127 (MQ=255)
                                  tATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGAGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCCGACGCTGa                                                                                                                                       <  8:284969/127‑1 (MQ=255)
                                    tGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGtt                                                                                                               <  1:331750/149‑1 (MQ=255)
                                      tCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGttt                                                                                                              <  4:241139/148‑1 (MQ=255)
                                       cGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCaa                                                                                                                    <  8:119672/141‑1 (MQ=255)
                                       cGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCaa                                                                                                                    >  7:119672/1‑141 (MQ=255)
                                        gCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTcgg                                                                                                           >  3:63105/1‑149 (MQ=255)
                                            gTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTcggcgg                                                                                                        <  5:285474/148‑1 (MQ=255)
                                               aTGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCggcggcggc                                                                                                    <  3:142928/149‑1 (MQ=255)
                                                 gtgCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCggcggcggc                                                                                                    >  8:342259/1‑147 (MQ=255)
                                                 gtgCCCGGGAACCTGTGCACCATGACGCGAGTGCAATCCACAACAGGGGATATCCAACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGCCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCggcggcggc                                                                                                    >  6:28280/1‑147 (MQ=255)
                                                  tgCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCggcggcggc                                                                                                    >  5:331404/1‑146 (MQ=255)
                                                  tgCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCaa                                                                                                  <  2:298648/148‑1 (MQ=255)
                                                    cccGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAaga                                                                                                <  6:25504/148‑1 (MQ=255)
                                                      cGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCt                                                                                             >  2:124486/1‑149 (MQ=255)
                                                        ggAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCg                                                                                           <  4:376471/149‑1 (MQ=255)
                                                          aaCCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCggg                                                                                         <  7:220813/149‑1 (MQ=255)
                                                             cTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTc                                                                                      >  7:205069/1‑149 (MQ=255)
                                                                tgCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTcc                                                                                     <  7:285452/147‑1 (MQ=255)
                                                                  cACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCgcgc                                                                                 <  8:227134/149‑1 (MQ=255)
                                                                   aCCATGACGCGCGTGCAATCCACAACAGGGGATATCACACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGcc                                                                                >  6:145132/1‑149 (MQ=255)
                                                                     cATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCTGCAAGAGCTCGGGCTCCGCGCCTa                                                                              <  2:224378/149‑1 (MQ=255)
                                                                     cATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTa                                                                              >  4:160227/1‑149 (MQ=255)
                                                                      aTGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCaa                                                                                                                    <  8:270665/110‑1 (MQ=255)
                                                                      aTGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTAAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCaa                                                                                                                    >  7:270665/1‑110 (MQ=255)
                                                                          cgcgcgTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATccac                                                                         <  7:194164/149‑1 (MQ=255)
                                                                          cgcgcgTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATccac                                                                         <  8:205069/149‑1 (MQ=255)
                                                                                 gCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGAccc                                                                  >  4:66201/1‑149 (MQ=255)
                                                                                   aaTCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGc                                                                >  6:334853/1‑149 (MQ=255)
                                                                                      ccACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTTGGCGGCGGCAAGAGCTCGGGCTCCgcgc                                                                                 >  6:105607/1‑129 (MQ=255)
                                                                                      ccACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTTGGCGGCGGCAAGAGCTCGGGCTCCgcgc                                                                                 <  5:105607/129‑1 (MQ=255)
                                                                                       cacaACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCa                                                              <  1:369244/147‑1 (MQ=255)
                                                                                       cacaACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCa                                                              >  2:369244/1‑147 (MQ=255)
                                                                                        acaacaGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATccacca                                                                       <  4:314477/137‑1 (MQ=255)
                                                                                        acaacaGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATccacca                                                                       >  3:314477/1‑137 (MQ=255)
                                                                                        acaacaGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGc                                                           >  6:170639/1‑149 (MQ=255)
                                                                                          aacaGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGATTCGGCTCCGGCAAGAGCTCGGGCTCCGCGCATATCCACCAGACCCGCCAGGCTa                                                         >  5:401643/1‑149 (MQ=255)
                                                                                               gggATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAgg                                                            <  2:362031/141‑1 (MQ=255)
                                                                                               gggATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAgg                                                            >  1:362031/1‑141 (MQ=255)
                                                                                                 gATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAAc                                                  <  7:197429/149‑1 (MQ=255)
                                                                                                      cccACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGc                                             >  8:15452/1‑149 (MQ=255)
                                                                                                              tGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTTACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCtgctgc                                      <  3:160227/148‑1 (MQ=255)
                                                                                                                cAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCcgc                                   >  5:93721/1‑149 (MQ=255)
                                                                                                                 aaCGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGTGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCcgcg                                  >  4:255825/1‑149 (MQ=255)
                                                                                                                 aaCGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCcgcg                                  >  8:236717/1‑149 (MQ=255)
                                                                                                                 aaCGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCcgcg                                  >  5:50864/1‑149 (MQ=255)
                                                                                                                    gTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACg                                                       >  6:306922/1‑125 (MQ=255)
                                                                                                                    gTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACg                                                       <  5:306922/125‑1 (MQ=255)
                                                                                                                     ttttCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGa                              <  6:331404/149‑1 (MQ=255)
                                                                                                                                tCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGc                   >  6:412979/1‑149 (MQ=255)
                                                                                                                                  gCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGcc                  <  7:332051/148‑1 (MQ=255)
                                                                                                                                   cACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGt                <  1:271026/149‑1 (MQ=255)
                                                                                                                                    aCTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTg               <  3:66201/149‑1 (MQ=255)
                                                                                                                                    aCTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTg               <  1:124486/149‑1 (MQ=255)
                                                                                                                                        gcgcTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCAcc           >  6:189910/1‑149 (MQ=255)
                                                                                                                                          gcTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCgg         <  5:334853/149‑1 (MQ=255)
                                                                                                                                          gcTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCgg         >  6:141142/1‑149 (MQ=255)
                                                                                                                                                cGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTgc     <  7:342259/147‑1 (MQ=255)
                                                                                                                                                 gTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTgccgc  <  4:63105/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
AACTTTGTAACATTCAAGCAAAAGCTCGAAACACTATGTCGCTGGTCATGTGCCCGGGAACCTGTGCACCATGACGCGCGTGCAATCCACAACAGGGGATATCCCACTTATGCAACGTTTTCTTAGCATCGCACTGGCGCTCTGCGTCGGCCTGACGCTGAGCCTGGACGCCAGCGCCAAGCGTTTCGGCGGCGGCAAGAGCTCGGGCTCCGCGCCTATCCACCAGACCCGCCAGGCTACGCCAACCACGCCTGCTGCCGCGCCGACCGTACCTGGCCGTGCACCGGCTGCCGC  >  NC_002947/6103059‑6103352

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: