Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,508,777 G→A 5.3% S113L (TCG→TTG) ‡ PP_2203 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,508,7770GA 5.3% 284.0 / 3.1 95S113L (TCG→TTG) ‡PP_2203hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (54/35);  new base A (2/3);  total (56/39)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.31e-01

CGGTAGCAGGGTGGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGA  >  NC_002947/2508630‑2508924
                                                                                                                                                   |                                                                                                                                                     
cGGTAGCAGGGTGGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGa                                                                                                                                                      <  6:69915/149‑1 (MQ=255)
cGGTAGCAGGGTGGAGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGa                                                                                                                                                      <  5:197908/149‑1 (MQ=255)
  gTAGCAGGGTGGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGc                                                                                                                                                    <  3:220474/149‑1 (MQ=255)
      cAGGGTGGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCAt                                                                                                                                                <  4:4252/149‑1 (MQ=255)
      cAGGGTGGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATg                                                                                                                                               <  7:214465/149‑1 (MQ=255)
          gTGGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGCTGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGc                                                                                                                                            <  1:256552/149‑1 (MQ=255)
          gTGGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCATGCCCATGGGc                                                                                                                                            <  7:271636/149‑1 (MQ=255)
              cggcACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGc                                                                                                                                       <  8:162118/149‑1 (MQ=255)
               ggcACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGcc                                                                                                                                      >  5:18341/1‑149 (MQ=255)
                gcACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGccc                                                                                                                                     <  7:184308/149‑1 (MQ=255)
                  aCCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGTCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCgg                                                                                                                                   <  6:209416/149‑1 (MQ=255)
                     aaaGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCATGCCCATGGGCTCGCCCGGCat                                                                                                                                <  1:238956/149‑2 (MQ=255)
                      aaGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCAcc                                                                                                                               >  4:122262/1‑149 (MQ=255)
                           gcTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAg                                                                                                                          >  7:6810/1‑149 (MQ=255)
                           gcTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAg                                                                                                                          <  5:190012/149‑1 (MQ=255)
                             tGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGc                                                                                                                        >  8:162235/1‑149 (MQ=255)
                                   aCAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGCCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGa                                                                                                                  <  2:226046/149‑1 (MQ=255)
                                   aCAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTCTCCTGGCCATCCTTGGTCACCCCGATGACCTGGAAGGCATCCTTCTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCATGcccgt                                                                                                                    >  8:304170/1‑145 (MQ=255)
                                    cAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAg                                                                                                                 <  3:122262/149‑1 (MQ=255)
                                    cAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGGCACCCATCTCCATGCCGGGCGAGCCCATGGTCATGCCCGGCACCGCCAGGCCCTTGAg                                                                                                                 <  3:54548/149‑1 (MQ=255)
                                      gcgcCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGt                                                                                                               >  8:204454/1‑149 (MQ=255)
                                       cgcCCCTTGGCTGAGCTACTGTAGTTGGCGACCAGGGTGTCCTGGCCTTCCTGGGTCACCCCGATGCCCTGGTAGGCATCCTTGTGGTCACCCATCTCAATGCCGGGCGAGCCCTTGGGCTTGCCCGGCACCGCCAGGCCCTTGAGGTc                                                                                                              <  4:53802/149‑1 (MQ=255)
                                        gcCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCCCCGCCATGCCCTTGAGGTCg                                                                                                             >  8:243604/1‑149 (MQ=255)
                                          cccTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAg                                                                                                                 >  7:138989/1‑143 (MQ=255)
                                                 cTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTg                                                                                                    <  7:207103/149‑1 (MQ=255)
                                                     gCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAg                                                                                                                 >  5:223609/1‑132 (MQ=255)
                                                      cATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGtt                                                                                                           >  7:296442/1‑137 (MQ=255)
                                                          agtagtTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAgg                                                                                           >  7:52073/1‑149 (MQ=255)
                                                          agtagtTGGCGACCACGGTGGCCTGCCCACCCGGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCAGCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAgg                                                                                           <  7:102242/149‑1 (MQ=255)
                                                            tagtTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTc                                                                                                              <  3:50806/128‑1 (MQ=255)
                                                            tagtTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTc                                                                                                              >  4:50806/1‑128 (MQ=255)
                                                                tGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCcagca                                                                                             >  4:283170/1‑141 (MQ=255)
                                                                tGGCGACCACGGGGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCcagca                                                                                             <  3:283170/141‑1 (MQ=255)
                                                                gggCGCCCCCGGGGTCCGGCCCATGCTGGGCCCCCCGTTTCACGTGGTAGGCACCATTGTGGCCACCAATCTCCATGCCGGGCGCGCCCTGGGGCTGGCCGGGCACCGCCAGGCCCTTGAGGTCTTTGCGTTTGGCCAGCAGGCGCAc                                                                                      <  6:18341/147‑1 (MQ=255)
                                                                  gCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGttt                                                                                                     >  5:276178/1‑131 (MQ=255)
                                                                  gCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTg                                                                                   >  1:44267/1‑149 (MQ=255)
                                                                  gCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTg                                                                                   >  3:135195/1‑149 (MQ=255)
                                                                   cGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGc                                                                                  >  3:151564/1‑149 (MQ=255)
                                                                   cGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCGTGAGGTCGTTGCGTTTGGCCAGCAGGCGCCCCTGc                                                                                  >  6:13407/1‑149 (MQ=255)
                                                                    gACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCt                                                                                 >  7:298305/1‑149 (MQ=255)
                                                                     aCCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTc                                                                                <  2:74120/149‑1 (MQ=255)
                                                                                cTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGt                                                                     >  7:240837/1‑149 (MQ=255)
                                                                                cTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGt                                                                     >  7:173553/1‑149 (MQ=255)
                                                                                 tGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTg                                                                   >  2:36767/1‑149 (MQ=255)
                                                                                 tGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTg                                                                   >  1:79683/1‑149 (MQ=255)
                                                                                 tGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTGGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTg                                                                   >  3:26796/1‑149 (MQ=255)
                                                                                   gCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGc                                                                 >  8:188593/1‑149 (MQ=255)
                                                                                         cTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGcc                                                                                                 <  5:255476/112‑1 (MQ=255)
                                                                                         cTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGcc                                                                                                 >  6:255476/1‑112 (MQ=255)
                                                                                            ggTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCAc                                                                                      >  6:114797/1‑120 (MQ=255)
                                                                                              tCACCCCGATGACGTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAAc                                                      <  2:199435/149‑1 (MQ=255)
                                                                                              cccccccCATCACCTGGTTGCCATCCTTTGGGCCACCCATCTCCATTCGGGGCGTGCCCTTGGGCTTCCCGGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAAc                                                      <  7:204454/146‑1 (MQ=255)
                                                                                                  cccGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTgc                                                  >  3:310172/1‑149 (MQ=255)
                                                                                                  cccGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTgc                                                  >  5:225969/1‑149 (MQ=255)
                                                                                                    cGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTgccg                                                >  4:202876/1‑149 (MQ=255)
                                                                                                           cTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCa                                         >  3:169434/1‑149 (MQ=255)
                                                                                                           cTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGCCCCCGCCAGGCCCTAGAGGTCGTTGCGTTTGGCCACCAGGAGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCTCCGATCa                                         >  4:236589/1‑149 (MQ=255)
                                                                                                            tGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTcgac                                                          <  5:77245/131‑1 (MQ=255)
                                                                                                            tGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTcgac                                                          >  6:77245/1‑131 (MQ=255)
                                                                                                            tGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCAc                                        >  5:184776/1‑149 (MQ=255)
                                                                                                              gTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCAccc                                      >  5:134300/1‑149 (MQ=255)
                                                                                                              gTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCAccc                                      <  7:188593/149‑1 (MQ=255)
                                                                                                              gTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCAccc                                      >  6:82170/1‑149 (MQ=255)
                                                                                                                  gCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGgt                                  >  3:293027/1‑149 (MQ=255)
                                                                                                                    aTCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGgtgt                                >  6:146856/1‑149 (MQ=255)
                                                                                                                    aTCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGgtgt                                >  5:265647/1‑149 (MQ=255)
                                                                                                                    aTCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGgcgt                                >  8:136954/1‑149 (MQ=255)
                                                                                                                       cTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGAACACCCCGGTGTgg                              <  2:44267/148‑1 (MQ=255)
                                                                                                                          gtgGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCac                           >  2:255762/1‑148 (MQ=255)
                                                                                                                             gTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCac                           <  5:150283/145‑1 (MQ=255)
                                                                                                                             gTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCac                           >  6:150283/1‑145 (MQ=255)
                                                                                                                              tCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCacgacg                       >  5:73496/1‑148 (MQ=255)
                                                                                                                                aCCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGc                      >  3:116849/1‑147 (MQ=255)
                                                                                                                                  ccATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAgg                  >  7:109391/1‑149 (MQ=255)
                                                                                                                                  ccATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAgg                  >  8:178666/1‑149 (MQ=255)
                                                                                                                                    aTCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAgg                  >  6:114368/1‑147 (MQ=255)
                                                                                                                                    aTCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAgg                  <  5:114368/147‑1 (MQ=255)
                                                                                                                                       tcCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTgg             >  2:36305/1‑149 (MQ=255)
                                                                                                                                        ccctGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTCGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTgg             <  1:255762/146‑1 (MQ=255)
                                                                                                                                         cATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCCCCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCACGGCGTGgcg           >  1:289168/1‑149 (MQ=255)
                                                                                                                                         cATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTgg                              <  1:241650/130‑1 (MQ=255)
                                                                                                                                         cATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTgg                              >  2:241650/1‑130 (MQ=255)
                                                                                                                                         cATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTgg             <  7:178666/147‑1 (MQ=255)
                                                                                                                                         cATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGgcg           >  1:290184/1‑149 (MQ=255)
                                                                                                                                          aTGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGgcgc          >  5:205248/1‑149 (MQ=255)
                                                                                                                                            gCCGGGCCGCCCCATGGGGCTGCCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGcc         <  4:151564/149‑1 (MQ=255)
                                                                                                                                            gCCGGGCATGCCCATGGGCTCGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTcgac                                                          <  7:237294/99‑1 (MQ=38)
                                                                                                                                            gCCGGGCATGCCCATGGGCTCGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTcgac                                                          >  8:237294/1‑99 (MQ=38)
                                                                                                                                             ccGGGCATGCCCATGGGCTCGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCAc       >  1:297764/1‑149 (MQ=255)
                                                                                                                                              cGGGCGCGCCCTGGGGCTGGCCGGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGCGTGGCGCCACGCCTTGCGTCGTGCCAca      <  4:310172/149‑2 (MQ=255)
                                                                                                                                               gggCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGc     >  6:211327/1‑149 (MQ=255)
                                                                                                                                                 gCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGttt                                                                                                     >  7:77256/1‑52 (MQ=255)
                                                                                                                                                 gCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGttt                                                                                                     <  8:77256/52‑1 (MQ=255)
                                                                                                                                                 gCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGCGTGGCGCCACGCCTTGCGTCGTGCCACACCg   <  7:136954/149‑1 (MQ=38)
                                                                                                                                                  cGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGCGTGGCGCCACGCCTTGCGTCGTGCCACACCgg  <  4:169434/149‑2 (MQ=38)
                                                                                                                                                   |                                                                                                                                                     
CGGTAGCAGGGTGGCGGCACCAAAGGCGCTGAGAAACAGCGCCCCTAGGCTGAGCATCAGTAGTTGGCGACCACGGTGTCCTGGCCATCCTGGGTCACCCCGATGACCTGGTAGGCATCCTTGTGGTCACCCATCTCCATGCCGGGCGAGCCCATGGGCATGCCCGGCACCGCCAGGCCCTTGAGGTCGTTGCGTTTGGCCAGCAGGCGCACCTGCTCGGCCGGCACGTGGCCTTCGACGAACTTGCCGCCGATCACCCCGGTGTGGCACGACGCAAGGCGTGGCGCCACGCCGA  >  NC_002947/2508630‑2508924

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: