Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 4,713,585 A→C T591P (ACC→CCC)  PP_4171 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,713,5850AC98.8% 564.9 / ‑3.1 166T591P (ACC→CCC) PP_4171hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base C (91/72);  total (93/73)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 168 total aligned reads displayed.

CTCGACCTGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCACCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTAT  >  NC_002947/4713441‑4713731
                                                                                                                                                |                                                                                                                                                  
 tcgacctGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTgc                                                                                                                                              >  1:2765219/1‑150 (MQ=255)
 tcgacctGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTgc                                                                                                                                              >  1:1453860/1‑150 (MQ=255)
 cctaccgGAGCGTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTgc                                                                                                                                              <  2:905518/147‑1 (MQ=255)
  cgacctGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTgcg                                                                                                                                             >  1:1739542/1‑150 (MQ=255)
       tGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGAcc                                                                                                                                        <  2:402912/150‑1 (MQ=255)
           ccTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACcgcct                                                                                                                                    >  2:2138149/1‑150 (MQ=255)
               gATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACcgcctgcgc                                                                                                                                <  1:1627428/150‑1 (MQ=255)
                   tgAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACcgcctgcgc                                                                                                                                >  2:379616/1‑146 (MQ=255)
                   tgAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACcgcctgcgc                                                                                                                                <  1:379616/146‑1 (MQ=255)
                      aGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCt                                                                                                                         >  2:2993055/1‑150 (MQ=255)
                       gCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTa                                                                                                                        >  1:1557940/1‑150 (MQ=255)
                       gCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTa                                                                                                                        >  2:3045060/1‑150 (MQ=255)
                        ctcttCTACAACGAGATCACCGGCATGAACAGCGAGCGTAATAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTAt                                                                                                                       <  2:1453860/147‑1 (MQ=255)
                          gTTCTACGCAGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGc                                                                                                                     <  1:165741/150‑1 (MQ=255)
                          gTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACCGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGc                                                                                                                     <  2:3397704/150‑1 (MQ=255)
                           ttCTATGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGcc                                                                                                                    >  2:56030/1‑150 (MQ=255)
                           ttCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGcc                                                                                                                    <  2:324397/150‑1 (MQ=255)
                            tCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCC‑CCCTGCGCGACCGCCTGCGCCTGACCTATGCCTa                                                                                                                  <  2:3306097/150‑1 (MQ=255)
                             cTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTa                                                                                                                  <  2:3612041/150‑1 (MQ=255)
                             cTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTa                                                                                                                  >  2:21037/1‑150 (MQ=255)
                                 gacgaGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTc                                                                                                              <  2:1557940/150‑1 (MQ=255)
                                 gacgaGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTc                                                                                                              >  1:829203/1‑150 (MQ=255)
                                  acgaGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCg                                                                                                             >  2:1446086/1‑150 (MQ=255)
                                   cgaGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTAt                                                                                                                       <  2:28247/139‑1 (MQ=255)
                                    gagaTCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGAt                                                                                                           >  1:821904/1‑150 (MQ=255)
                                    gagaTCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGAt                                                                                                           >  2:2619488/1‑150 (MQ=255)
                                      gaTCACCGGCATGATCAGCGAGCCTCTGAGGAACGATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGc                                                                                                         >  1:2924996/1‑150 (MQ=255)
                                        tCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCt                                                                                                       <  2:2579350/150‑1 (MQ=255)
                                          aCCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTgg                                                                                                     <  1:2877006/150‑1 (MQ=255)
                                           ccGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGTGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTggg                                                                                                    >  2:2197833/1‑150 (MQ=255)
                                             ggCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGcc                                                                                                  >  1:3445881/1‑150 (MQ=255)
                                               cATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAg                                                                                                <  1:2138149/150‑1 (MQ=255)
                                                aTGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGc                                                                                               <  1:3045060/150‑1 (MQ=255)
                                                aTGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGc                                                                                               <  1:1383195/150‑1 (MQ=255)
                                                  gATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCaa                                                                                             <  1:1775893/150‑1 (MQ=255)
                                                    tCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAAcc                                                                                           >  2:1131047/1‑150 (MQ=255)
                                                    tCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAAcc                                                                                           >  1:3203434/1‑150 (MQ=255)
                                                      aGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAAcccc                                                                                         >  1:1098153/1‑150 (MQ=255)
                                                         gAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCcgac                                                                                      <  1:546386/150‑1 (MQ=255)
                                                             ctctGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGAcc                                                                                  <  2:2309774/150‑1 (MQ=255)
                                                              tctGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCcga                                                                                       <  2:653945/144‑1 (MQ=255)
                                                                tGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTc                                                                               >  1:355058/1‑150 (MQ=255)
                                                                tGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTc                                                                               <  2:3391564/150‑1 (MQ=255)
                                                                tGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTc                                                                               <  1:2836258/150‑1 (MQ=255)
                                                                  aGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCa                                                                                                 >  2:2857628/1‑130 (MQ=255)
                                                                    gAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGCCCTAGGCCTACGTCGATGCCTGGGCCAGCAACCCCGCCGACGGTCTccg                                                                           <  2:1739542/150‑2 (MQ=255)
                                                                     aaCAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTc                                                                          >  1:757626/1‑150 (MQ=255)
                                                                          tCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTc                                                                          <  2:1333082/145‑1 (MQ=255)
                                                                          tCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTc                                                                          >  1:1333082/1‑145 (MQ=255)
                                                                          tCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAgcgc                                                                     >  1:732794/1‑150 (MQ=255)
                                                                            aGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGCCACACGCGCCAGCAACCCCGACGACCGTCGCCTCAGCGccc                                                                   >  2:2045429/1‑150 (MQ=255)
                                                                             gggcATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCACCTCAGCGCCCg                                                                  >  2:1991406/4‑150 (MQ=255)
                                                                              tACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACcgcctgcg                                                                                                                                 >  2:2482579/1‑86 (MQ=255)
                                                                              tACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACcgcctgcg                                                                                                                                 <  1:2482579/86‑1 (MQ=255)
                                                                                 aTCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGTTCAACCGCCTGCGACTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGGGCCCGCAAc                                                              >  1:2007945/1‑150 (MQ=255)
                                                                                 aTCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAAc                                                              >  1:2394430/1‑150 (MQ=255)
                                                                                 aTCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAAc                                                              <  1:2507693/150‑1 (MQ=255)
                                                                                  tCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACa                                                             >  1:3390857/1‑150 (MQ=255)
                                                                                  tCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACa                                                             >  1:1503959/1‑150 (MQ=255)
                                                                                      cAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCgg                                                         >  1:1152489/1‑150 (MQ=255)
                                                                                      cAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCgg                                                         >  1:1876069/1‑150 (MQ=255)
                                                                                       aGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGc                                                        >  1:1989718/1‑150 (MQ=255)
                                                                                             gCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGcc                                                  >  2:420781/1‑150 (MQ=255)
                                                                                             gCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGcc                                                  >  1:2026032/1‑150 (MQ=255)
                                                                                              ccAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCg                                                 >  1:3299164/1‑150 (MQ=255)
                                                                                               caacaaGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGAcc                                                                                  <  2:684596/116‑1 (MQ=255)
                                                                                               caacaaGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGAcc                                                                                  >  1:684596/1‑116 (MQ=255)
                                                                                               caacaaGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCgg                                                >  1:483839/1‑150 (MQ=255)
                                                                                                 acaaGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGt                                              <  2:3445881/150‑1 (MQ=255)
                                                                                                  caaGGCCCGTTTCAGCGGCAGCGAGACGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGTGCCAGCAACACCGACGACCGTCGCATCAGCGCCCGCGACAGCGGCTCGGCCGCGTg                                             >  2:1256043/1‑150 (MQ=255)
                                                                                                   aaGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTgg                                            >  2:1678230/1‑150 (MQ=255)
                                                                                                     ggCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGAt                                          >  2:1001841/1‑150 (MQ=255)
                                                                                                      gCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGt                                                                                                               <  1:3294017/80‑1 (MQ=255)
                                                                                                      gCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCaa                                                               <  1:361241/128‑1 (MQ=255)
                                                                                                      gCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATg                                         >  1:3267054/1‑150 (MQ=255)
                                                                                                       cccGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATgc                                        >  1:3176639/1‑150 (MQ=255)
                                                                                                        ccGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTggg                                                                                                    >  1:3437309/1‑89 (MQ=255)
                                                                                                         cGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATgcgc                                      >  2:357224/1‑150 (MQ=255)
                                                                                                               aGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGTCTGCGCCTGCCCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGc                                                                             <  2:1818422/105‑1 (MQ=255)
                                                                                                               aGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATgg                                <  1:21037/150‑1 (MQ=255)
                                                                                                               aGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATgg                                >  1:627026/1‑150 (MQ=255)
                                                                                                               aGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGACTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGc                                                                             >  1:1818422/1‑105 (MQ=255)
                                                                                                                gCGGCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATggg                               <  2:3293909/150‑1 (MQ=255)
                                                                                                                  ggCAGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGt                             <  2:2924996/150‑1 (MQ=255)
                                                                                                                     aGCGAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGcc                          <  2:1989718/150‑1 (MQ=255)
                                                                                                                        gAGGCGCAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGc                       <  2:2015803/150‑1 (MQ=255)
                                                                                                                           gcgcAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTgg                    >  1:2672040/1‑150 (MQ=255)
                                                                                                                           gcgcAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTgg                    <  1:1300125/150‑1 (MQ=255)
                                                                                                                             gcAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATg                                                                                                          <  1:1908312/62‑1 (MQ=255)
                                                                                                                             gcAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCATCCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTc                  >  1:234810/1‑150 (MQ=255)
                                                                                                                             gcAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTc                  <  1:1001841/150‑1 (MQ=255)
                                                                                                                             gcAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTc                  <  2:3176639/150‑1 (MQ=255)
                                                                                                                             gcAAATGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTc                  <  1:56030/150‑1 (MQ=255)
                                                                                                                               aaaTGGACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGa                <  2:2389661/150‑1 (MQ=255)
                                                                                                                                   ggACTGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAgcgc            <  2:82522/150‑1 (MQ=255)
                                                                                                                                       tGGCGCCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTg        >  1:1631484/1‑150 (MQ=255)
                                                                                                                                          cgcCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGCCGGCCGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTc     <  1:216592/150‑1 (MQ=255)
                                                                                                                                           gcCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCcga                                                                                       >  2:2737850/1‑67 (MQ=255)
                                                                                                                                           gcCCCCCCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCcga                                                                                       <  1:2737850/67‑1 (MQ=255)
                                                                                                                                             ccccccCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTAt  >  1:911501/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CTCGACCTGAGCCTGGATGTGAAGCTGTTCTACGACGAGATCACCGGCATGATCAGCGAGCCTCTGAGGAACAATCAGTACATCGCCAGCAACGCCAACAAGGCCCGTTTCAGCGGCAGCGAGGCGCAAATGGACTGGCGCCCCACCCTGCGCGACCGCCTGCGCCTGACCTATGCCTACGTCGATGCCTGGGCCAGCAACCCCGACGACCGTCGCCTCAGCGCCCGCAACAGCGGCTCGGCCGGGTGGATGCGCGAATGGGGTGCCGGCTGGTCGAGCGCGCTGTTCTAT  >  NC_002947/4713441‑4713731

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