Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 421,165 T→C 5.5% H928R (CAC→CGC)  PP_0348 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947421,1650TC 5.5% 268.9 / 4.3 91H928R (CAC→CGC) PP_0348hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (42/43);  new base C (3/2);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.08e-01

GCACGCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAATCTCCAAGTGG  >  NC_002947/421017‑421311
                                                                                                                                                     |                                                                                                                                                   
gCACGCCACCCATCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGt                                                                                                                                                     <  7:104956/149‑1 (MQ=255)
  aCGCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTgg                                                                                                                                                   >  1:382836/1‑149 (MQ=255)
       aCCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCcgg                                                                                                                                             <  4:323608/149‑1 (MQ=255)
       aCCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCcgg                                                                                                                                             <  4:287004/149‑1 (MQ=255)
       aCCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCcgg                                                                                                                                             >  4:195993/1‑149 (MQ=255)
         ccAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCcggcg                                                                                                                                           <  5:386832/149‑1 (MQ=255)
         ccAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCcggcg                                                                                                                                           >  3:9183/1‑149 (MQ=255)
          cagctcagctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCcggcgg                                                                                                                                          >  6:303554/1‑149 (MQ=255)
           agctcagctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCggt                                                                                                                                         >  6:236610/1‑149 (MQ=255)
              tcagctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCTGGCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTg                                                                                                                                                    >  1:282557/1‑137 (MQ=255)
              tcagctTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTg                                                                                                                                                    <  2:282557/136‑1 (MQ=255)
                    tCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCt                                                                                                                                >  1:77074/1‑149 (MQ=255)
                     cGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCt                                                                                                                                <  1:78767/148‑1 (MQ=255)
                        cccGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTgg                                                                                                                              >  3:254762/1‑147 (MQ=255)
                        cccGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCa                                                                                                                            <  1:136468/149‑1 (MQ=255)
                          cGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAgt                                                                                                                          <  4:308280/149‑1 (MQ=255)
                          cGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAgt                                                                                                                          <  7:327416/149‑1 (MQ=255)
                             ccTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTc                                                                                                                       <  4:156667/149‑1 (MQ=255)
                             ccTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTc                                                                                                                       >  1:153577/1‑149 (MQ=255)
                                 tGCCCGTGACCGGGCACGGCTGCT‑GCGCATCTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGc                                                                                                                    <  1:391529/148‑1 (MQ=255)
                                   cccGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAg                                                                                                                 >  8:315922/1‑149 (MQ=255)
                                   cccGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGC‑CCACTGACAGCGGAGCcggcg                                                                                                                                           <  5:205702/123‑1 (MQ=255)
                                   cccGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGC‑CCACTGACAGCGGAGCcggcg                                                                                                                                           >  6:205702/1‑123 (MQ=255)
                                       tGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCAc                                                                                                             <  7:83425/149‑1 (MQ=255)
                                         aCCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCAct                                                                                                            <  7:350464/148‑1 (MQ=255)
                                          ccGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCActgc                                                                                                          >  5:253826/1‑149 (MQ=255)
                                           cGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCActgct                                                                                                         <  8:157746/149‑1 (MQ=255)
                                               cACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTgcgc                                                                                                     <  2:153577/149‑1 (MQ=255)
                                                aCGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGAGCTGCCAGGGACAGCCAGGCACTGATGcccc                                                                                                    >  6:84785/1‑149 (MQ=255)
                                                aCGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGACCGCCGGC‑CCACTGACAGCGGAGCcggcggc                                                                                                                                         <  3:138122/112‑2 (MQ=37)
                                                aCGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGACCGCCGGC‑CCACTGACAGCGGAGCcggcggc                                                                                                                                         >  4:138122/1‑111 (MQ=37)
                                                       ct‑gcGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCGGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGGCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCcctgcct                                                                                              <  2:339426/148‑1 (MQ=255)
                                                        t‑gcGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTgg                                                                                            >  1:246130/1‑149 (MQ=255)
                                                          gcgcAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGc                                                                                           <  7:315922/149‑1 (MQ=255)
                                                              aaCTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGccagt                                                                                       <  2:146830/149‑5 (MQ=255)
                                                                cTGGCCCACTGGCAGTGCCTTGCAGCAGCGGCTGGACCGGCTGTTGCCAAAACTGGGCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAAcgc                                                                                     <  4:254762/149‑1 (MQ=255)
                                                                 tGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAAcgcg                                                                                    >  4:15158/1‑149 (MQ=255)
                                                                     ccACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTAt                                                                                 >  4:124728/1‑148 (MQ=255)
                                                                      cACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCa                                                                               <  4:9183/149‑1 (MQ=255)
                                                                        cTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGc                                                                             >  2:183738/1‑149 (MQ=255)
                                                                            cAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACt                                                                         >  1:258337/1‑149 (MQ=255)
                                                                             aGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGACTGGCTGACActgcca                                                                                                        >  6:378783/1‑115 (MQ=255)
                                                                             aGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGACTGGCTGACActgcca                                                                                                        <  5:378783/117‑3 (MQ=255)
                                                                                 ccTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCg                                                                    >  1:361754/1‑149 (MQ=255)
                                                                                     gCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCGGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGttt                                                                 <  7:136073/148‑1 (MQ=255)
                                                                                     gCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGtttt                                                                >  1:79882/1‑149 (MQ=255)
                                                                                       aTCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCg                                                              >  1:370411/1‑149 (MQ=255)
                                                                                        tCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGt                                                             >  2:151802/1‑149 (MQ=255)
                                                                                        tCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGt                                                             >  3:136506/1‑149 (MQ=255)
                                                                                          aGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTTGCGAACGCGTATCAGCCACTGGCCGTTTTCGTgg                                                           <  2:387909/149‑1 (MQ=255)
                                                                                           gCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGgc                                                          >  8:141195/1‑149 (MQ=255)
                                                                                           gCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGgc                                                          <  1:183738/149‑1 (MQ=255)
                                                                                                tGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCgcagc                                                     >  2:387058/1‑149 (MQ=255)
                                                                                                  gACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCgcagca                                                    <  1:20803/148‑1 (MQ=255)
                                                                                                    ccGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGcc                                                                     <  4:84711/129‑1 (MQ=255)
                                                                                                    ccGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGcc                                                                     >  3:84711/1‑129 (MQ=255)
                                                                                                    ccGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCAcca                                                 <  3:360487/149‑2 (MQ=255)
                                                                                                      ggCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGcc                                               <  3:124728/149‑1 (MQ=255)
                                                                                                             gCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGca                                        >  4:362826/1‑149 (MQ=255)
                                                                                                             gCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGca                                        >  3:194557/1‑149 (MQ=255)
                                                                                                              ccAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGcac                                       <  2:382836/149‑1 (MQ=255)
                                                                                                              ccAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGcac                                       >  8:111531/1‑149 (MQ=255)
                                                                                                               cAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGcacc                                      >  8:53296/1‑149 (MQ=255)
                                                                                                               cAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGcacc                                      >  2:41561/1‑149 (MQ=255)
                                                                                                               cAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGcac                                       >  4:147304/1‑148 (MQ=255)
                                                                                                                aaaaCTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGcacca                                     >  8:128097/1‑149 (MQ=255)
                                                                                                                aaaaCTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGACGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGcacca                                     >  6:415376/1‑149 (MQ=255)
                                                                                                                       gTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACcgcg                                <  6:287811/147‑1 (MQ=255)
                                                                                                                       gTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGc                              <  3:362826/149‑1 (MQ=255)
                                                                                                                        tCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCg                             >  5:360153/1‑149 (MQ=255)
                                                                                                                          gACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGt                            <  2:258337/148‑1 (MQ=255)
                                                                                                                          gACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCACTGTCAGCCAGTCACTGCTGCGCCCTGCCTGACGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGt                            <  5:415376/148‑1 (MQ=255)
                                                                                                                           aCTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTg                          <  5:303554/149‑1 (MQ=255)
                                                                                                                            cTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGc                         <  3:269327/149‑1 (MQ=255)
                                                                                                                             tGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCt                        >  8:41074/1‑149 (MQ=255)
                                                                                                                               aGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCt                        <  2:77074/147‑1 (MQ=255)
                                                                                                                               aGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTg                       <  2:246130/148‑1 (MQ=255)
                                                                                                                               aGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTg                       <  5:236610/148‑1 (MQ=255)
                                                                                                                                ggCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGc                     <  7:128097/149‑1 (MQ=255)
                                                                                                                                caccgccgGC‑CCCCTGACAGCGGAGCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTgg                      >  1:257650/3‑148 (MQ=255)
                                                                                                                                   cgccgGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAAc                  <  6:253826/149‑1 (MQ=255)
                                                                                                                                   cgccgGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAAc                  >  5:335999/1‑149 (MQ=255)
                                                                                                                                          tCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACgg                <  7:182834/144‑1 (MQ=255)
                                                                                                                                          tCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACgg                >  8:182834/1‑144 (MQ=255)
                                                                                                                                          tCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAAtc           >  2:315229/1‑149 (MQ=255)
                                                                                                                                            cGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAAtctc         >  7:152724/1‑149 (MQ=255)
                                                                                                                                            cGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAAtctc         >  6:15682/1‑149 (MQ=255)
                                                                                                                                               tGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAATCTCCaa      <  2:370411/149‑1 (MQ=255)
                                                                                                                                                 tCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCTACGGCACCACCCCGGCGTTGCTGGCAACGGCAATCTCCAAg     <  2:79882/148‑1 (MQ=255)
                                                                                                                                                 tCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAATCTCCaa      <  5:241028/147‑1 (MQ=255)
                                                                                                                                                  cAGGGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAATCTCCAAGTg   <  5:84785/149‑1 (MQ=255)
                                                                                                                                                   aGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAATCTCCAAGTgg  >  2:190128/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCACGCCACCCAGCTCAGCTTCGCCCCGGCCTATGCCCGTGACCGGGCACGGCTGCT‑GCGCAACTGGCCCACTGGCAGTGCCTTGCATCAGCGGCTGGACCGGCTGTTGCCAAAACTGGTCGACTGAGGCCGCCGGCTCCGCTGTCAGTGG‑GCCGGCGGTGGTTGATCTGGCAGTGTCAGCCAGTCACTGCTGCGCCCTGCCTGGCGAACGCGTATCAGCCACTGGCCGTTTTCGTGGCGCAGCACCGCCCACGGCACCACCCCGGCGTTGCTGGCAACGGCAATCTCCAAGTGG  >  NC_002947/421017‑421311

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