Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 226,907 C→G 7.4% A61G (GCC→GGC)  PP_0175 → transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947226,9070CG 7.4% 213.9 / 5.8 67A61G (GCC→GGC) PP_0175transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (36/26);  new base G (3/2);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGCCGTGCCCTGGCCGAAGCCATTGAAGC  >  NC_002947/226761‑227054
                                                                                                                                                  |                                                                                                                                                   
tGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGccgccg                                                                                                                                                   >  8:196374/1‑149 (MQ=255)
          gTCGCTGCCCGCCACTGGCGGCGCATCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAg                                                                                                                                         >  7:301202/1‑149 (MQ=255)
             gCTGCCCGCCACTGGCGGCGCATCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGccc                                                                                                                                      >  1:99409/1‑149 (MQ=255)
             gCTGCCCGCCACTGGCGGCGCATCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGccc                                                                                                                                      >  5:343504/1‑149 (MQ=255)
                  ccGCCACTGGCGGCGCATCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGGCGGCCTggg                                                                                                                                            >  6:103758/1‑137 (MQ=255)
                  ccGCCACTGGCGGCGCATCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGGCGGCCTggg                                                                                                                                            <  5:103758/138‑2 (MQ=255)
                   cGCCACTGGCGGCGCATCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCt                                                                                                                                >  2:240289/1‑149 (MQ=255)
                    gCCACTGGCGGCGCATCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTg                                                                                                                               <  3:106426/149‑1 (MQ=255)
                    gCCACTGGCGGCGCATCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTg                                                                                                                               >  2:77294/1‑149 (MQ=255)
                      cACTGGCGGCGCATCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGt                                                                                                                             <  2:99409/149‑1 (MQ=255)
                             ggcgCATCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTg                                                                                                                      >  8:146954/1‑149 (MQ=255)
                              gcgcATCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTCTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGc                                                                                                                     <  1:240289/149‑1 (MQ=255)
                                  aTCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGa                                                                                                                 >  7:69538/1‑149 (MQ=255)
                                  aTCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCACCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCACCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGa                                                                                                                 >  2:70304/1‑149 (MQ=255)
                                      gTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTg                                                                                                                               <  6:154808/131‑1 (MQ=255)
                                      gTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTg                                                                                                                               >  5:154808/1‑131 (MQ=255)
                                       tCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGAc                                                                                                                <  2:156619/145‑1 (MQ=255)
                                         aTACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACtg                                                                                                          >  5:49382/1‑149 (MQ=255)
                                         aTACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACtg                                                                                                          >  5:299174/1‑149 (MQ=255)
                                           aCCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACtgtg                                                                                                        >  5:19413/1‑149 (MQ=255)
                                            ccGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACtgtg                                                                                                        <  4:261010/148‑1 (MQ=255)
                                             cGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGc                                                                                                       >  4:291111/1‑148 (MQ=255)
                                             cGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGc                                                                                                       >  2:120861/1‑148 (MQ=255)
                                               ccTTGACCGGCTACGGCATTTCCGAGGCCTGCGCAGCGCCTATGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTAGCTGGTACGCCTTCTTGACCAACTGTGCa                                                                                                      >  5:267125/1‑147 (MQ=255)
                                                     ccGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGc                                                                                              <  7:196374/149‑1 (MQ=255)
                                                     ccGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGc                                                                                              <  6:19413/149‑1 (MQ=255)
                                                        gCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTg                                                                                           >  6:111732/1‑149 (MQ=255)
                                                        gCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTg                                                                                           <  7:146954/149‑1 (MQ=255)
                                                             ggCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGtgtg                                                                                      >  8:66413/1‑149 (MQ=255)
                                                              gCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGc                                                                                     >  8:31955/1‑149 (MQ=255)
                                                               cATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCaa                                                                                                             <  6:48728/124‑1 (MQ=255)
                                                               cATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGcc                                                                                    <  8:301202/149‑1 (MQ=255)
                                                                      gAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCg                                                                              >  6:335760/1‑148 (MQ=255)
                                                                      gAGGCCTGCGCTGCGCCTATGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCg                                                                              >  2:266845/1‑148 (MQ=255)
                                                                       aGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGaa                                                                            <  6:343504/149‑1 (MQ=255)
                                                                       aGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGaa                                                                            <  8:69538/149‑1 (MQ=255)
                                                                         gCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGa                                                                          >  6:285683/1‑149 (MQ=255)
                                                                               gctgcgcCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACt                                                                    <  5:111732/149‑1 (MQ=255)
                                                                               gctgcgcCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACt                                                                    <  7:66413/149‑1 (MQ=255)
                                                                                ctgcgcCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTg                                                                   >  7:141609/1‑149 (MQ=255)
                                                                                    gcCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATc                                                               >  6:201210/1‑149 (MQ=255)
                                                                                       tttGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATcgcc                                                            >  4:33877/1‑149 (MQ=255)
                                                                                       tttGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATcgcc                                                            >  4:187141/1‑149 (MQ=255)
                                                                                         tGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATcgcc                                                            <  3:236376/147‑1 (MQ=255)
                                                                                          gCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCcgcg                                                         >  2:299386/1‑149 (MQ=255)
                                                                                          gCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCcgcg                                                         >  2:111138/1‑149 (MQ=255)
                                                                                           cTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCgcgc                                                        <  1:77294/149‑1 (MQ=255)
                                                                                                gatCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCa                                                      >  4:14845/1‑146 (MQ=255)
                                                                                                gatCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCa                                                      <  3:14845/146‑1 (MQ=255)
                                                                                                 atCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAgg                                                   <  6:267125/148‑1 (MQ=255)
                                                                                                       cgcCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTgc                                             <  5:285683/148‑1 (MQ=255)
                                                                                                          cTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCt                                         <  6:49382/149‑1 (MQ=255)
                                                                                                                         cACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACAAACTAGCCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGcc                           <  1:111138/148‑1 (MQ=255)
                                                                                                                            cAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGCCGTGc                       >  2:151098/1‑149 (MQ=255)
                                                                                                                                  gCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTgcgc                                           >  1:314074/1‑123 (MQ=255)
                                                                                                                                  gCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTgcgc                                           <  2:314074/123‑1 (MQ=255)
                                                                                                                                    tgtgGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTgcgc                                           <  6:67870/121‑1 (MQ=255)
                                                                                                                                    tgtgGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTgcgc                                           >  5:67870/1‑121 (MQ=255)
                                                                                                                                    tgtgGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGCCGTGCCCTGGcc                >  5:91842/1‑148 (MQ=255)
                                                                                                                                      cttcccAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGcc                           <  6:343450/132‑1 (MQ=255)
                                                                                                                                        ttccAGGCCGGCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGCCGTGCCCTGGCCGAAGc           >  3:212081/3‑149 (MQ=255)
                                                                                                                                         tccAGGCCGGCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAg                              <  5:182781/128‑1 (MQ=255)
                                                                                                                                         tccAGGCCGGCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAg                              >  6:182781/2‑129 (MQ=255)
                                                                                                                                          ccAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGCCGTGCCCTGGCCGAAGCCa         >  4:19306/1‑149 (MQ=255)
                                                                                                                                          ccAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGCCGTGCCCTGGCCGAAGCCa         >  6:149240/1‑149 (MQ=255)
                                                                                                                                           cAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGCCGTGCCCTGGcc                <  8:164984/141‑1 (MQ=255)
                                                                                                                                           cAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGCCGTGCCCTGGcc                >  7:164984/1‑141 (MQ=255)
                                                                                                                                             ggCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGCCGTGCCCTGGCCGAAGCCAtt       >  8:238288/1‑148 (MQ=255)
                                                                                                                                               ccgccgGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGCCGTGCCCTGGCCGAAGCCATTGaa    <  5:201210/149‑1 (MQ=255)
                                                                                                                                                cgccgGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGCCGTGCCCTGGCCGAAGCCATTGAAg   <  3:33877/149‑1 (MQ=255)
                                                                                                                                                cgccgGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGCCGTGCCCTGGCCGAAGCCATTGAAg   <  3:187141/149‑1 (MQ=255)
                                                                                                                                                 gccgGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGCCGTGCCCTGGCCGAAGCCATTGAAGc  <  1:120861/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGGCCTGGTAGTCGCTGCCCGCCACTGGCGGCGCATCTGTCATACCGCCTTGACCGGCTACGGCATTTCCGAGGCCTGCGCTGCGCCTTTGCTGATGATCGTGCGCCTGGGCGACGGCGTGCACCAGGTGGCTGTGGCCCAGGCCGCCGGCCTGGAAAGCCCGTCGCTGGTACGCCTGCTTGACCAACTGTGCAAAGCTGGCCTGGTGTGCCGCAGCGAAGACCCACTGGATCGCCGCGCCAAGGCACTGCGCCTGACTGCCGAAGGCCGTGCCCTGGCCGAAGCCATTGAAGC  >  NC_002947/226761‑227054

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: ‑