Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,139,671 A→T 6.2% R165R (CGT→CGA PP_5481 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,139,6710AT 6.2% 309.6 / 7.5 108R165R (CGT→CGAPP_5481hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (49/50);  new base T (3/4);  total (54/54)
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

Only 100 of 118 total aligned reads displayed.

GCTACGGCATTCTAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAG  >  NC_002947/2139522‑2139820
                                                                                                                                                     |                                                                                                                                                     
gCTACGGCATTCTAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGcga                                                                                                                                                       <  1:520867/150‑1 (MQ=255)
  tACGGCATTCTAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGcgacg                                                                                                                                                     >  1:1856102/1‑150 (MQ=255)
     ggCATTCTAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGacgacg                                                                                                                                                  >  2:1320191/1‑150 (MQ=255)
        tttCTGGGGGTGTGGGGTCTCAGCATTGGCTCGGCAATCCCCTTCGTGGGGGTGGCCAGTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCa                                                                                                                                               >  2:1089578/2‑150 (MQ=255)
             aGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGcc                                                                                                                                          <  1:2219106/150‑1 (MQ=255)
               gggTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCAGTGGCGTCGTCGCCAATGCCgg                                                                                                                                        <  1:1650743/150‑1 (MQ=255)
                    tCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGcgac                                                                                                                                                      >  2:1704789/1‑131 (MQ=255)
                    tCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGcgac                                                                                                                                                      <  1:1704789/131‑1 (MQ=255)
                    tCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGcccac                                                                                                                                   >  1:569420/1‑150 (MQ=255)
                    tCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCAGTGGCGTCGTCGCCAATGCCggtgcaa                                                                                                                                   <  2:1120677/150‑6 (MQ=255)
                      ggATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGCGTGGAGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGcccaccc                                                                                                                                 >  1:2155070/1‑150 (MQ=255)
                            cAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGcgac                                                                                                                                                      >  2:369979/1‑123 (MQ=255)
                            cAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGcgac                                                                                                                                                      <  1:369979/123‑1 (MQ=255)
                             aGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAgcg                                                                                                                          >  1:438314/1‑150 (MQ=255)
                              gCATTGGCACGGCGATCACCTTAATGCTGGTGGCGACTGCGGTGCTCTACGTCGTTTGCGTTCGCTTGTTGGTCAGCGGCCGCGGCATACGGGCCTGAAGCCTTGCACCGGCCTTGTCGACGACGCCCCTGCCGGCCCCCCCACCAgcgc                                                                                                                         >  2:346501/1‑150 (MQ=255)
                                aTTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCa                                                                                                                       >  1:437828/1‑150 (MQ=255)
                                aTTGGCACGGCGATCACCTTCATGCTGGTGGCAAATGCGGTTGTTTACATTGTTGGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCa                                                                                                                       >  1:1573503/1‑150 (MQ=255)
                                   ggCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTa                                                                                                                    <  1:2158189/150‑1 (MQ=255)
                                    gCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGc                                                                                                                                                 <  2:1184673/120‑1 (MQ=255)
                                      aCGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTCCGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGaa                                                                                                                 >  1:2227995/1‑150 (MQ=255)
                                           gATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACAtgg                                                                                                            <  2:1894297/150‑1 (MQ=255)
                                            aTCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACAtggt                                                                                                           <  2:1518626/150‑1 (MQ=255)
                                             tCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACAtggtg                                                                                                          <  2:569420/150‑1 (MQ=255)
                                               aCCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGga                                                                                                        <  2:2012985/150‑1 (MQ=255)
                                               aCCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGga                                                                                                        >  1:1012634/1‑150 (MQ=255)
                                                ccTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCgg                                                                                                                                        <  2:2092903/117‑1 (MQ=255)
                                                ccTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCgg                                                                                                                                        >  1:2092903/1‑117 (MQ=255)
                                                ccTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGgag                                                                                                       <  2:1014144/150‑1 (MQ=255)
                                                 cTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGgaga                                                                                                      >  1:1338479/1‑150 (MQ=255)
                                                 cTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAAAACATGGTGgaga                                                                                                      >  1:132224/1‑150 (MQ=255)
                                                 cTTCATGCTGGTGGCCACTGCGGTGCTCCACGTTGCATGAGCTAGCTCGCTGGACAGCGGCCGCGGCAGGCGGCCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGgaga                                                                                                      >  1:1799190/1‑150 (MQ=255)
                                                  ttCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGaa                                                                                                     <  2:1472051/150‑1 (MQ=255)
                                                  ttCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGaa                                                                                                     <  2:1732738/150‑1 (MQ=255)
                                                  ttCATGCTGGTGGCCACGGCGGTGCTCTACGTTGTTGGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGAGTCGACGCCACTGCCGGCCCTCCCACCAGCGCCAGTAGAACATGGTGGAGaa                                                                                                     <  1:1008193/150‑1 (MQ=255)
                                                      tGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAg                                                                                                 >  2:1698049/1‑150 (MQ=255)
                                                       gCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAgg                                                                                                <  1:1080303/150‑1 (MQ=255)
                                                        cTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGc                                                                                               <  2:1338479/150‑1 (MQ=255)
                                                        cTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGc                                                                                               <  1:1618706/150‑1 (MQ=255)
                                                         tggtggCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTa                                                                                                  <  2:520872/146‑1 (MQ=255)
                                                          ggtggCCACTGTGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCaa                                                                                             <  2:1799190/150‑1 (MQ=255)
                                                          ggtggCCACTGCGGTGCTCTACGTTGTTTGCCTCCGACTTCGTGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCaa                                                                                             <  2:2301199/150‑1 (MQ=255)
                                                              gCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGAt                                                                                         >  1:1511723/1‑150 (MQ=255)
                                                                cACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGCCGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGc                                                                                       >  2:1791522/1‑150 (MQ=255)
                                                                cACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGc                                                                                       <  2:1511723/150‑1 (MQ=255)
                                                                cACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGc                                                                                       >  1:137903/1‑150 (MQ=255)
                                                                cACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGc                                                                                       <  1:1089578/150‑1 (MQ=255)
                                                                 aCTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGcc                                                                                      >  1:960814/1‑150 (MQ=255)
                                                                 aCTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGcc                                                                                      >  2:1223126/1‑150 (MQ=255)
                                                                 aCTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGcc                                                                                      <  1:1223126/150‑1 (MQ=255)
                                                                  cTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGccc                                                                                     <  1:1698049/150‑1 (MQ=255)
                                                                         gCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGcacca                                                                              <  2:2155070/150‑1 (MQ=255)
                                                                          ctctACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGcaccac                                                                             >  1:481019/1‑150 (MQ=255)
                                                                           tctACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCAcg                                                                            >  1:437183/1‑150 (MQ=255)
                                                                            ctACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGTCGGGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCCCCCCCCACCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACTAcgc                                                                           <  2:163630/150‑1 (MQ=255)
                                                                                 ttgttTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGcac                                                                      >  1:1222150/1‑150 (MQ=255)
                                                                                    tttGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCAcacc                                                                   >  1:1818490/1‑150 (MQ=255)
                                                                                         gTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCAccga                                                              >  1:1143711/1‑150 (MQ=255)
                                                                                         gTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCAccga                                                              >  1:1736101/1‑150 (MQ=255)
                                                                                          ttCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGACACTGTCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCAccgag                                                             >  1:1603197/1‑150 (MQ=255)
                                                                                           tCGCTGGGTGGGCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCGGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCAccgagc                                                            <  2:2241106/150‑1 (MQ=255)
                                                                                             gCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGGCGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCAccgagccg                                                          >  2:196842/1‑150 (MQ=255)
                                                                                              cTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCGCCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCAcaccac                                                                 <  1:1448547/142‑1 (MQ=255)
                                                                                              cTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCAcaccac                                                                 >  2:1448547/1‑142 (MQ=255)
                                                                                               ttGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGccc                                                                                     <  1:1984060/121‑1 (MQ=255)
                                                                                               ttGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGccc                                                                                     >  2:1984060/1‑121 (MQ=255)
                                                                                                 gCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAACCTGGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGNANTNGGCAAGGATGCCCANNACCACGCCGCACACCACCGAGCCGAGCa                                                      <  1:1180914/150‑1 (MQ=255)
                                                                                                     gTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGaa                                                  <  2:1856102/150‑1 (MQ=255)
                                                                                                      tCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAg                                                 >  2:370749/1‑150 (MQ=255)
                                                                                                       cAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAAAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAgg                                                <  2:132224/150‑1 (MQ=255)
                                                                                                            gCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGc                                           <  2:137903/150‑1 (MQ=255)
                                                                                                                  gCATGGGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGTAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGc                                     <  1:346501/150‑1 (MQ=255)
                                                                                                                    aTGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGt                                   <  2:1143711/150‑1 (MQ=255)
                                                                                                                    aTGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGt                                   <  2:107946/150‑1 (MQ=255)
                                                                                                                     tGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCATTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTc                                  <  2:437828/150‑1 (MQ=255)
                                                                                                                       cGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGa                                <  2:858491/150‑1 (MQ=255)
                                                                                                                       cGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGa                                <  2:437183/150‑1 (MQ=255)
                                                                                                                       cGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGa                                <  2:2042485/150‑1 (MQ=255)
                                                                                                                           ccTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGc                            >  1:393860/1‑150 (MQ=255)
                                                                                                                           ccTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGc                            >  2:243746/1‑150 (MQ=255)
                                                                                                                            cTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCt                           <  1:370749/150‑1 (MQ=255)
                                                                                                                              gAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGAGGAGTATTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTgg                         >  2:2200557/1‑150 (MQ=255)
                                                                                                                                    ttGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGcc                                                                                      <  2:2174871/83‑1 (MQ=255)
                                                                                                                                    ttGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGcc                                                                                      >  1:2174871/1‑83 (MQ=255)
                                                                                                                                    ttGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTg                   >  1:1932184/1‑150 (MQ=255)
                                                                                                                                    ttGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTg                   >  1:763064/1‑150 (MQ=255)
                                                                                                                                    ttGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTg                   >  1:2129853/1‑150 (MQ=255)
                                                                                                                                        aCCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGAt                                                                                         <  1:341440/76‑1 (MQ=255)
                                                                                                                                        aCCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGAt                                                                                         >  2:341440/1‑76 (MQ=255)
                                                                                                                                          cGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGaa                                                                                                     <  2:668161/62‑1 (MQ=255)
                                                                                                                                          cGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGaa                                                                                                     >  1:668161/1‑62 (MQ=255)
                                                                                                                                          cGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAg                                                 >  1:376586/1‑114 (MQ=255)
                                                                                                                                          cGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAg                                                 <  2:376586/114‑1 (MQ=255)
                                                                                                                                          cGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAg                                                 <  2:376745/114‑1 (MQ=255)
                                                                                                                                          cGGCAGTGGCGTCGTCGCCAATGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCAc             >  1:626083/1‑150 (MQ=255)
                                                                                                                                            gCAGTGGCGTCGTCGCCAATGCCGGCCCACCCACCAGCGCCAGTAGAACGTGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACcgagccgagc                                                       <  1:2077429/106‑1 (MQ=255)
                                                                                                                                            gCAGTGGCGTCGTCGCCAATGCCGGCCCACCCACCAGCGCCAGTAGAACGTGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACcgagccgagc                                                       >  2:2077429/1‑106 (MQ=255)
                                                                                                                                             cagtGGCGTCGTCGCCAATGCCGGTGCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACAATGCCAGGCATTAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTAATCCACTgg          >  2:1815134/4‑150 (MQ=255)
                                                                                                                                               ttGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGaa        >  1:607514/1‑150 (MQ=255)
                                                                                                                                                    gacgacGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCa   <  2:607514/150‑1 (MQ=255)
                                                                                                                                                     acgacgCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAg  >  1:226636/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTACGGCATTCTAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAGCGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAG  >  NC_002947/2139522‑2139820

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: