Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,385,129 C→T 5.3% E111K (GAG→AAG)  PP_2986 ← oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,385,1290CT 5.3% 263.9 / 2.3 90E111K (GAG→AAG) PP_2986oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (42/42);  new base T (2/3);  total (45/45)
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

CGAGGATCTTGCTGCGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTCGACCTCGGCCAT  >  NC_002947/3384985‑3385276
                                                                                                                                                |                                                                                                                                                   
cGAGGATCTTGCTGCGGTAGTAGTCATCGGCGTCGGCGACCACTCTCCGCCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCGATGGACGCCGTATGGATGAACCATAGTATGCCGGCCCCGTAGGCCTGCTCAAGCa                                                                                                                                                <  2:514511/150‑1 (MQ=255)
 gAGGATCTTGCTGCGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAg                                                                                                                                               >  1:1258286/1‑150 (MQ=255)
   ggATCTTGCTGCGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTTGAGCAGgc                                                                                                                                             <  1:1876912/150‑1 (MQ=255)
            tgcGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTAc                                                                                                                                    >  2:336527/1‑150 (MQ=255)
             gcGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTAcc                                                                                                                                   <  1:594503/150‑1 (MQ=255)
              cGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCAATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTAccc                                                                                                                                  <  2:534454/150‑1 (MQ=255)
                   gtagtCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGacac                                                                                                                             <  1:1847843/150‑1 (MQ=255)
                     agtCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCTGGTGCGCCGTTGAGCCCAGCAATGGAAGACGTATGGATGAACCGTCGTCTGCAGTCCCCGTAGTCCGGCTCAAGCAAGCGCCGCGTACCCGCCGCGt                                                                                                                           <  2:884029/150‑1 (MQ=255)
                       tCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGttt                                                                                                                         >  2:2313341/1‑150 (MQ=255)
                        cATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTa                                                                                                                        <  1:1436652/150‑1 (MQ=255)
                             gcgtcggcgtcgGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGt                                                                                                                   <  2:2297565/150‑1 (MQ=255)
                             gcgtcggcgtcgGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGt                                                                                                                   <  1:1431797/150‑1 (MQ=255)
                               gtcggcgtcggcCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCTTCTTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTc                                                                                                                 <  1:2246246/150‑1 (MQ=255)
                                  ggcgtcggcCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAg                                                                                                              >  2:1486399/1‑150 (MQ=255)
                                    cggcggcCCCCAGGCAGTTCTCCTCGATGAACGTGCCAGGTGCGCCGTTGAGCACCGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGtt                                                                                                            <  2:1711089/147‑1 (MQ=255)
                                      tcggcCCTCAGGCAGACCGGCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCt                                                                                                          >  1:1289070/1‑150 (MQ=255)
                                          cccTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTcc                                                                                                      >  2:1111771/1‑150 (MQ=255)
                                           ccTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCa                                                                                                     >  2:833225/1‑150 (MQ=255)
                                             tCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGt                                                                                                   <  2:2286234/150‑1 (MQ=255)
                                             tCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGAACCGTTTCTCTGTTCGAGTTCTTGCCAGt                                                                                                   <  2:1841723/150‑1 (MQ=255)
                                               aGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTgg                                                                                                 >  2:653353/1‑150 (MQ=255)
                                                 gCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTACTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCt                                                                                               >  2:972426/1‑150 (MQ=255)
                                                       ctcctcGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTTGAGCAGgcctac                                                                                                                                         <  1:1601080/102‑5 (MQ=255)
                                                       ctcctcGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTTGAGCAGgcctac                                                                                                                                         >  2:1601080/1‑98 (MQ=255)
                                                        tcctcGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGccc                                                                                        >  2:970217/1‑150 (MQ=255)
                                                         cctcGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCt                                                                                       >  1:1238550/1‑150 (MQ=255)
                                                          ctcGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCtt                                                                                      >  2:513002/1‑150 (MQ=255)
                                                           tcGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTg                                                                                     <  2:1675263/150‑1 (MQ=255)
                                                           tcGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTg                                                                                     >  1:495166/1‑150 (MQ=255)
                                                           tcGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTg                                                                                     >  2:1778339/1‑150 (MQ=255)
                                                             ggtGGCCGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGTCGTACCGATCTGCCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTa                                                                                   <  1:760052/148‑1 (MQ=255)
                                                             gATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTa                                                                                   <  1:336527/150‑1 (MQ=255)
                                                             gATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTa                                                                                   <  1:513002/150‑1 (MQ=255)
                                                              aTGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAg                                                                                  >  1:587899/1‑150 (MQ=255)
                                                                    gTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATcg                                                                            >  2:302266/1‑150 (MQ=255)
                                                                    gTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCATTGTAATTATcg                                                                            >  2:1210063/1‑150 (MQ=255)
                                                                     tGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATcgc                                                                           <  2:495166/150‑1 (MQ=255)
                                                                     tGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATcgc                                                                           <  2:1238550/150‑1 (MQ=255)
                                                                      gCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATcgcg                                                                          >  1:1595441/1‑150 (MQ=255)
                                                                         aGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGaaa                                                                       <  2:1595441/150‑1 (MQ=255)
                                                                         aGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGaaa                                                                       <  2:587899/150‑1 (MQ=255)
                                                                          ggTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAg                                                                                                                                               >  1:855194/1‑77 (MQ=255)
                                                                          ggTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAg                                                                                                                                               <  2:855194/77‑1 (MQ=255)
                                                                                ccGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTc                                                                                                           >  2:208037/1‑107 (MQ=255)
                                                                                ccGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTc                                                                                                           <  1:208037/107‑1 (MQ=255)
                                                                                  gTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGc                                                              <  1:2009006/150‑1 (MQ=255)
                                                                                    tGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCtg                                                            >  1:1818010/1‑150 (MQ=255)
                                                                                       gCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCGGCCCTTGTAGTTATCGCGAAAGAACGACGCtgtgt                                                         <  1:193186/150‑1 (MQ=255)
                                                                                       gCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCtgtgt                                                         >  2:1615322/1‑150 (MQ=255)
                                                                                           aGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGaaa                                                     >  2:1050108/1‑150 (MQ=255)
                                                                                                  gAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGt                                                                                                   <  1:964619/97‑1 (MQ=255)
                                                                                                  gAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGt                                                                                                   >  2:964619/1‑97 (MQ=255)
                                                                                                  gAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCt                                                             <  2:1970162/135‑1 (MQ=255)
                                                                                                  gAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCt                                                             >  1:1970162/1‑135 (MQ=255)
                                                                                                    aGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATc                                            <  2:1258286/150‑1 (MQ=255)
                                                                                                     gACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCGTTCTGTAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCg                                           >  2:1667221/1‑150 (MQ=255)
                                                                                                     gACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATc                                            <  2:1274648/149‑1 (MQ=255)
                                                                                                     gACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATc                                            >  1:1274648/1‑149 (MQ=255)
                                                                                                           tGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCt                                     >  2:141761/1‑150 (MQ=255)
                                                                                                           tGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCt                                     <  2:1818010/150‑1 (MQ=255)
                                                                                                            ggATGAACCGTCGTATGCCGGCCCCTTAGGCCTGCTCAAGCAGTCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTGTCGCAACCtt                                    <  2:971743/150‑1 (MQ=255)
                                                                                                                gAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAg                                <  2:1407717/150‑1 (MQ=255)
                                                                                                                  aCCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCg                              >  1:106943/1‑150 (MQ=255)
                                                                                                                  aCCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCg                              >  2:1116231/1‑150 (MQ=255)
                                                                                                                   ccGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGa                             >  2:1620515/1‑150 (MQ=255)
                                                                                                                   ccGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGa                             >  1:30338/1‑150 (MQ=255)
                                                                                                                    cgtcgtATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGAt                            <  1:302266/150‑1 (MQ=255)
                                                                                                                        gtATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATcgc                                                                           <  1:2161289/99‑1 (MQ=255)
                                                                                                                        gtATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATcgc                                                                           >  2:2161289/1‑99 (MQ=255)
                                                                                                                         tATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGAtgctgc                       <  1:849950/150‑1 (MQ=255)
                                                                                                                         tATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGAtgctgc                       <  1:833225/150‑1 (MQ=255)
                                                                                                                            gCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGaa                    >  1:2261652/1‑150 (MQ=255)
                                                                                                                            gCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGaa                    <  1:1050108/150‑1 (MQ=255)
                                                                                                                              cGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATg                  >  1:2093951/1‑150 (MQ=255)
                                                                                                                              cGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATg                  <  1:2313341/150‑1 (MQ=255)
                                                                                                                               ggCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCt                                                                                       <  1:41783/80‑1 (MQ=255)
                                                                                                                               ggCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCt                                                                                       >  2:41783/1‑80 (MQ=255)
                                                                                                                                 ccccGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGt               >  2:1383219/1‑150 (MQ=255)
                                                                                                                                     gTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTCGAc           <  1:653353/150‑1 (MQ=255)
                                                                                                                                     gTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTCGAc           <  1:1486399/150‑1 (MQ=255)
                                                                                                                                         gCCTGCTTGAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAAt                   >  2:1299029/1‑138 (MQ=255)
                                                                                                                                         gCCTGCTTGAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAAt                   <  1:1299029/138‑1 (MQ=255)
                                                                                                                                         gCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCtgt                                                           >  2:2083918/1‑98 (MQ=255)
                                                                                                                                         gCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCtgt                                                           <  1:2083918/98‑1 (MQ=255)
                                                                                                                                          ccTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTCGACCTCgg      <  1:141761/150‑1 (MQ=255)
                                                                                                                                           cTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTCGACCTCGGc     >  2:1279179/1‑150 (MQ=255)
                                                                                                                                           cTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGAGTCGACCTCGGc     >  1:415172/1‑150 (MQ=255)
                                                                                                                                            tGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTCGACCTCGGcc    >  1:1538359/1‑150 (MQ=255)
                                                                                                                                            tGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTCGACCTCGGcc    >  2:990765/1‑150 (MQ=255)
                                                                                                                                              cTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAACCACCGTATCCCAACCTTGCAGCGATGCTGCAAAGGCGTCTACCCCGGCCAt  >  2:126198/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CGAGGATCTTGCTGCGGTAGTAGTCATCGGCGTCGGCGTCGGCCCTCAGGCAGTTCTCCTCGATGGACGTGCCAGGTGCGCCGTTGAGCACAGCAATGGAAGACGTATGGATGAACCGTCGTATGCCGGCCCCGTAGGCCTGCTCAAGCAGGCGCCGCGTACCCGACACGTTTATCTGTTCGAGTTCTTTCCAGTGGCTGCCGCCCTTGTAGTTATCGCGAAAGAACGACGCTGTGTGAAACACCGTATCGCAACCTTGCAGCGATGCTGCGAATGCGTCGACCTCGGCCAT  >  NC_002947/3384985‑3385276

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