Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 133,399 C→G 5.9% H59Q (CAC→CAG PP_0128 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947133,3990CG 5.9% 263.9 / 4.6 82H59Q (CAC→CAGPP_0128hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (39/38);  new base G (3/2);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTGGCCCACTTGGGG  >  NC_002947/133251‑133546
                                                                                                                                                     |                                                                                                                                                   
ttGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGAATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCAc                                                                                                                                                     >  4:196510/1‑149 (MQ=255)
    ccTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACctgc                                                                                                                                                 <  6:169403/149‑1 (MQ=255)
       gCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCtgctgc                                                                                                                                              >  2:73193/1‑149 (MQ=255)
        cACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGcc                                                                                                                                             >  7:91964/1‑149 (MQ=255)
         aCGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTAACGCCATTACCTGACCCGTAGCTGGCAGCACCTGCTGCCg                                                                                                                                            >  5:205096/1‑149 (MQ=255)
                   ggTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGc                                                                                                                                  >  7:130863/1‑149 (MQ=255)
                    gTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCg                                                                                                                                            <  1:335373/138‑1 (MQ=255)
                    gTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCGGCTGCCg                                                                                                                                            >  2:335373/1‑138 (MQ=255)
                     tCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGt                                                                                                                                <  3:362752/149‑1 (MQ=255)
                     tCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGt                                                                                                                                <  2:191850/149‑1 (MQ=255)
                           aGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACCGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCgg                                                                                                                           <  1:73193/149‑1 (MQ=255)
                             cATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCAGGTGCTGcc                                                                                                                                             >  6:176252/1‑128 (MQ=255)
                             cATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCAGGTGCTGcc                                                                                                                                             <  5:176252/128‑1 (MQ=255)
                               tCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAAcc                                                                                                                      >  8:26883/1‑149 (MQ=255)
                               tCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAAcc                                                                                                                      >  8:21340/1‑149 (MQ=255)
                                   gggCAGGCCCGGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGCGCGTGCCGGTAACCTGc                                                                                                                   >  8:340235/3‑148 (MQ=255)
                                   cTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGc                                                                                                                                  >  2:350630/1‑133 (MQ=255)
                                   cTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGc                                                                                                                                  >  2:356396/1‑133 (MQ=255)
                                   cTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGc                                                                                                                                  <  1:350630/133‑1 (MQ=255)
                                   cTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGc                                                                                                                                  <  1:356396/133‑1 (MQ=255)
                                      cAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCag                                                                                                                 >  6:7352/1‑147 (MQ=255)
                                                      ccGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGAGGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGACGGTAACCTGCAGAAGATCGGCAGGctgctg                                                                                               >  7:246461/1‑149 (MQ=255)
                                                           gATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGa                                                                                          >  6:94994/1‑149 (MQ=255)
                                                                       cGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCt                                                                                                                     <  5:119862/110‑1 (MQ=255)
                                                                       cGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCtg                                                                                <  8:91964/147‑1 (MQ=255)
                                                                       cGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCCCTACCGCCATTACCTGACCCCTAGCTGGCAGCACCTTCTGCCGCATAACGGGCGTGCCGGTAAcca                                                                                                                     >  6:119862/1‑109 (MQ=255)
                                                                             ctCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGccc                                                                         <  6:205096/148‑1 (MQ=255)
                                                                             ctCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCt                                                                        >  3:141644/1‑149 (MQ=255)
                                                                                 gTTTCAACATCCAGGTGGGCATCATCACCGAGCGCTACCGCCATTACCTTACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAg                                                                    <  7:340235/149‑1 (MQ=255)
                                                                                 gTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTc                                                                                      <  4:108626/131‑1 (MQ=255)
                                                                                 gTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTc                                                                                      >  3:108626/1‑131 (MQ=255)
                                                                                  tttCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAg                                                                    >  6:244386/1‑148 (MQ=255)
                                                                                   ttCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAgg                                                                   >  7:227882/1‑148 (MQ=255)
                                                                                    tCAACCTCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAgg                                                                   >  7:346077/1‑147 (MQ=255)
                                                                                     cAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAg                                                                >  6:135102/1‑149 (MQ=255)
                                                                                      aaCATCCCGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACTCGTCGCTGTCAGCACCTGCTGCCGCATAACGTGCGTGCCGGTAACCTGGAGATGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGc                                                               <  6:200529/149‑1 (MQ=255)
                                                                                      aaCATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTg                                                                       <  8:157160/141‑1 (MQ=255)
                                                                                      aaCATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTg                                                                       >  7:157160/1‑141 (MQ=255)
                                                                                      aaCATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGc                                                               <  8:130863/149‑1 (MQ=255)
                                                                                       aCATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAgg                                                                   <  6:266648/144‑1 (MQ=255)
                                                                                       aCATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAgg                                                                   >  5:266648/1‑144 (MQ=255)
                                                                                         aTCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGa                                                            <  7:206473/149‑1 (MQ=255)
                                                                                         aTCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGa                                                            <  5:94994/149‑1 (MQ=255)
                                                                                         aTCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGa                                                            >  5:69132/1‑149 (MQ=255)
                                                                                                     aTCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCt                                                >  8:336433/1‑149 (MQ=255)
                                                                                                            ccGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCg                                         <  5:244386/149‑1 (MQ=255)
                                                                                                            ccGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCg                                         >  4:356368/1‑149 (MQ=255)
                                                                                                            ccGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCg                                         <  4:141644/149‑1 (MQ=255)
                                                                                                             cGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGATGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCgg                                        >  7:46157/1‑149 (MQ=255)
                                                                                                              gAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCggg                                       <  3:356368/149‑1 (MQ=255)
                                                                                                                  gcTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTAc                                   <  7:26883/149‑1 (MQ=255)
                                                                                                                  gcTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTAc                                   <  7:21340/149‑1 (MQ=255)
                                                                                                                     aCCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCtggtgg                                                                            >  6:131027/1‑105 (MQ=255)
                                                                                                                     aCCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCtggtgg                                                                            <  5:131027/105‑1 (MQ=255)
                                                                                                                       cGCCTTAACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCAGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCaa                              >  2:112531/1‑149 (MQ=255)
                                                                                                                       cGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCa                                                    >  8:159496/1‑127 (MQ=255)
                                                                                                                       cGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCa                                                    <  7:159496/127‑1 (MQ=255)
                                                                                                                        gCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAAc                             <  5:7352/149‑1 (MQ=255)
                                                                                                                         ccATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGAc                                                                                   >  1:319644/1‑94 (MQ=255)
                                                                                                                         ccATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGAc                                                                                   <  2:319644/94‑1 (MQ=255)
                                                                                                                         ccATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAAc                             >  5:8794/1‑148 (MQ=255)
                                                                                                                          cATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCa                           <  8:346077/149‑1 (MQ=255)
                                                                                                                           aTTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAg                          <  5:135102/149‑1 (MQ=255)
                                                                                                                           aTTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAg                          <  8:246461/149‑1 (MQ=255)
                                                                                                                           aTTACCTGACGCGTAGCTGGCAGCACCTGCTGCCACATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAg                          >  5:180990/1‑149 (MQ=255)
                                                                                                                            ttACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAgg                         >  3:289061/1‑149 (MQ=255)
                                                                                                                              aCCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGGCCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGt                        <  8:39988/148‑1 (MQ=255)
                                                                                                                                cTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACgg                                                                                                                                    >  4:183335/1‑38 (MQ=255)
                                                                                                                                cTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACgg                                                                                                                                    <  3:183335/38‑1 (MQ=255)
                                                                                                                                cTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCg                      >  3:333586/1‑148 (MQ=255)
                                                                                                                                 tGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGCAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGaa                    >  6:79508/1‑149 (MQ=255)
                                                                                                                                      gcgTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCtggtgg                                                                            <  1:254476/89‑1 (MQ=255)
                                                                                                                                      gcgTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCtggtgg                                                                            >  2:254476/1‑89 (MQ=255)
                                                                                                                                         tAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTgg             <  1:362684/149‑1 (MQ=255)
                                                                                                                                         tAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTgg             <  8:46157/149‑1 (MQ=255)
                                                                                                                                           gCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTggtt           >  8:107406/1‑147 (MQ=255)
                                                                                                                                            gcggCAGCAGGTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGc                                                                                                        <  2:147938/53‑1 (MQ=255)
                                                                                                                                            gcggCAGCAGGTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGc                                                                                                        >  1:147938/3‑55 (MQ=38)
                                                                                                                                            cTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTggttg          <  7:107406/149‑4 (MQ=255)
                                                                                                                                             cggCAGCAGGTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTGGCCCa         >  2:38400/2‑149 (MQ=255)
                                                                                                                                                 agcaCCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTGGcc           <  8:108129/143‑1 (MQ=255)
                                                                                                                                                 agcaCCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTGGcc           >  7:108129/1‑143 (MQ=255)
                                                                                                                                                    aCCTGCTGCCACATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTGGCCCACTTgggg  <  6:180990/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTGGCCTGCACGAGCCGCAGGTCAACGAGCATCACCTGCAGGCCCCGGGCCTGCCCGAGGATGGCCGTCTGCGGCTGCTCAGTTTCAACATCCAGGTGGGCATCAGCACCGAGCGCTACCGCCATTACCTGACGCGTAGCTGGCAGCACCTGCTGCCGCATAACGGGCGTGCCGGTAACCTGCAGAAGATCGGCAGGCTGCTGGGCGACTTCGACCTGGTGGCCCTGCAGGAAGCCGATGGTGGCAGCCTGCGCTCGGGCTACGTCAACCAGGTCGAACACCTGGCCCACTTGGGG  >  NC_002947/133251‑133546

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: