Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 81,984 G→C 8.0% E98D (GAG→GAC)  PP_0071 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00294781,9840GC 8.0% 262.6 / 10.9 87E98D (GAG→GAC) PP_0071hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/45);  new base C (3/4);  total (38/49)
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.27e-01

ACGTAGCCCTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTC  >  NC_002947/81836‑82130
                                                                                                                                                    |                                                                                                                                                  
aCGTAGCCCTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAg                                                                                                                                                    >  3:150775/1‑149 (MQ=255)
  gTAGCCCTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAGCCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTc                                                                                                                                                  >  7:331451/1‑149 (MQ=255)
  gTAGCCCTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTc                                                                                                                                                  <  5:183367/149‑1 (MQ=255)
     gCCCTATCCCTGCGCAGCGCTTCCGCAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGAGCGAGGCGATACTCGATGGCATCGAGTCGAc                                                                                                                                               <  5:347908/149‑1 (MQ=255)
     gCCCTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGAc                                                                                                                                               <  6:381581/149‑1 (MQ=255)
      cccTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGAc                                                                                                                                               <  7:268251/148‑1 (MQ=255)
      cccTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACg                                                                                                                                              <  7:295564/149‑1 (MQ=255)
       ccTAAACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCTATGGCATCGAGTCGACGc                                                                                                                                             >  4:20569/1‑149 (MQ=255)
            cccTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACAACATGCGGTAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCGTCGACTCGATGCCATCg                                                                                                                                        <  8:209069/148‑1 (MQ=255)
             ccTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGTGTGTGATCGAGGCGATACTCGATGGCGTCGACTCGATGCCATCga                                                                                                                                       <  8:185108/149‑2 (MQ=255)
             ccTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGAGTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCGTCGACTCGATGCCATCga                                                                                                                                       <  7:74471/149‑2 (MQ=255)
             ccTGCGAATCGCTTACGAAACCGGCTACAAATCTTTTGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCTAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCTATACTCGATTGCGTCGACTCGATGCCATCga                                                                                                                                       <  1:374607/149‑2 (MQ=255)
              cTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGcc                                                                                                                                      >  1:27581/1‑149 (MQ=255)
              cTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGcc                                                                                                                                      >  3:218561/1‑149 (MQ=255)
                gcgcAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGAGTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGt                                                                                                                                    >  6:22726/1‑149 (MQ=255)
                  gcAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCa                                                                                                                                  <  4:178442/149‑1 (MQ=255)
                   cAGCGCTTACGAAACCGGATACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCAc                                                                                                                                 >  1:341722/1‑149 (MQ=255)
                       gcTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGcc                                                                                                                              >  8:202344/1‑148 (MQ=255)
                        cTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCa                                                                                                                            <  8:253311/149‑1 (MQ=255)
                         ttACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAg                                                                                                                           >  5:115145/1‑149 (MQ=255)
                           aCGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGcc                                                                                                                         >  1:45556/1‑149 (MQ=255)
                           aCGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGc                                                                                                                          >  2:375528/1‑148 (MQ=255)
                             gAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTc                                                                                                                       >  4:2931/1‑149 (MQ=255)
                                aCCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCaa                                                                                                                    <  5:42197/149‑1 (MQ=255)
                                aCCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCaa                                                                                                                    <  7:202414/149‑1 (MQ=255)
                                      tACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCa                                                                                                                >  1:190277/1‑147 (MQ=255)
                                        cAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAgc                                                                                                            <  2:94158/149‑1 (MQ=255)
                                         aaaGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAgcg                                                                                                           <  6:128510/149‑1 (MQ=255)
                                         aaaGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCACCAAGGCATAgcg                                                                                                           >  4:57195/1‑149 (MQ=255)
                                          aaGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAgcgc                                                                                                          <  8:263624/149‑1 (MQ=255)
                                           aGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCgct                                                                                                         >  6:220079/1‑149 (MQ=255)
                                             cTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCgctg                                                                                                        <  1:375528/148‑1 (MQ=255)
                                                  gacgtaaTGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGTTACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATc                                                                                                   <  6:259093/142‑1 (MQ=255)
                                                  gACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCTTGCTATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCGAATGCATAGCGCTGCTATc                                                                                                   <  3:214885/148‑1 (MQ=255)
                                                  gACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATc                                                                                                   <  3:182946/148‑1 (MQ=255)
                                                   aCTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGa                                                                                                 <  3:20569/149‑1 (MQ=255)
                                                   aCTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGa                                                                                                 <  3:6094/149‑1 (MQ=255)
                                                    cTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAg                                                                                                <  6:285906/149‑1 (MQ=255)
                                                       cATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAgg                                                                                               <  3:284105/147‑1 (MQ=255)
                                                        aTGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAgcg                                                                                                           <  2:298581/134‑1 (MQ=255)
                                                        aTGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAgcg                                                                                                           >  1:298581/1‑134 (MQ=255)
                                                        aTGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGcc                                                                                             >  1:101286/1‑148 (MQ=255)
                                                         tGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGa                                                                                           <  8:324869/149‑1 (MQ=255)
                                                         tGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGa                                                                                           <  4:330144/149‑1 (MQ=255)
                                                           cGTGAAGAAGGCGGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACt                                                                                         <  8:83077/149‑1 (MQ=255)
                                                           cGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACt                                                                                         <  7:202344/149‑1 (MQ=255)
                                                             tGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCa                                                                                       <  3:99194/149‑1 (MQ=255)
                                                                     gCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAgg                                                                                                                  <  5:161993/114‑1 (MQ=255)
                                                                     gCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGACCAGCCTCCAAgg                                                                                                                  >  6:161993/1‑114 (MQ=255)
                                                                                  cTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTg                                                                  >  7:37927/1‑149 (MQ=255)
                                                                                        gAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATTGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCcgcg                                                            <  4:218561/149‑1 (MQ=255)
                                                                                        gAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCcgc                                                             >  5:197630/1‑148 (MQ=255)
                                                                                        gAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCcgc                                                             >  1:39790/1‑148 (MQ=255)
                                                                                        gAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCcgc                                                             <  3:44432/148‑1 (MQ=255)
                                                                                        gAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAACGCATAGCGCTCGTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCcgc                                                             >  7:311567/1‑148 (MQ=255)
                                                                                         aaCCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCcgcg                                                            >  1:180519/1‑148 (MQ=255)
                                                                                            cTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAaca                                                        >  3:292804/1‑149 (MQ=255)
                                                                                               cATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAacac                                                       <  2:190277/147‑1 (MQ=255)
                                                                                               cATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACAcc                                                      <  2:101286/148‑1 (MQ=255)
                                                                                                 tgggATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTTTCGAGGCCGACTCAGTACCTCTCCGAGCGCGCGAGCCCGCGAACACCgg                                                    >  1:110488/4‑148 (MQ=255)
                                                                                                   cGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCa                                                  <  3:2931/148‑1 (MQ=255)
                                                                                                     aTTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAg                                               <  4:292804/149‑1 (MQ=255)
                                                                                                          gaTGGCAGCGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGt                                          <  7:72071/149‑1 (MQ=255)
                                                                                                          gaTGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGt                                          <  5:22726/149‑1 (MQ=255)
                                                                                                           aTGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTg                                         >  7:196281/1‑149 (MQ=255)
                                                                                                              gCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGcc                                       <  5:232229/148‑1 (MQ=255)
                                                                                                                agagGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGccatc                                    >  7:28413/1‑149 (MQ=255)
                                                                                                                 gagGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGc                                                   >  1:220315/1‑133 (MQ=255)
                                                                                                                 gagGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGc                                                   <  2:220315/133‑1 (MQ=255)
                                                                                                                 gagGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCTTTGCCGGTgttcgcg                                   <  2:27581/149‑7 (MQ=39)
                                                                                                                          aTCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGt                           <  8:256869/148‑1 (MQ=255)
                                                                                                                             gAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGct                       >  5:48816/1‑149 (MQ=255)
                                                                                                                               ggCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGcttc                     >  3:349551/1‑149 (MQ=255)
                                                                                                                                gCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGcttct                    >  6:32977/1‑149 (MQ=255)
                                                                                                                                 cGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGcttctt                   <  8:37927/149‑1 (MQ=255)
                                                                                                                                     aCTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTcgcg               >  7:114416/1‑149 (MQ=255)
                                                                                                                                     aCTCGATGCCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGAGGTGGCTTCTTcgcg               >  2:221258/1‑149 (MQ=255)
                                                                                                                                      cTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCgg              <  8:114416/149‑1 (MQ=255)
                                                                                                                                        cGATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCggg             >  1:290417/1‑148 (MQ=255)
                                                                                                                                         gATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCACGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGt            >  6:3197/1‑148 (MQ=255)
                                                                                                                                          aTGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTaa          <  2:180519/149‑1 (MQ=255)
                                                                                                                                           tGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTaaa         >  1:232672/1‑149 (MQ=255)
                                                                                                                                             gCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAAcc       <  5:220079/149‑1 (MQ=255)
                                                                                                                                             gCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAAc        <  2:341722/148‑1 (MQ=255)
                                                                                                                                              cATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAAccc      >  5:32582/1‑149 (MQ=255)
                                                                                                                                                  gAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTc  >  8:39333/1‑149 (MQ=255)
                                                                                                                                                   aGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGt                           >  1:376160/1‑123 (MQ=255)
                                                                                                                                                   aGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGt                           <  2:376160/123‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACGTAGCCCTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTC  >  NC_002947/81836‑82130

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