Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,010,502 T→A 6.4% E325D (GAA→GAT) ‡ PP_1792 ← group 2 family glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,010,5020TA 6.4% 274.7 / 6.9 89E325D (GAA→GAT) ‡PP_1792group 2 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (38/45);  new base A (4/2);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAGGCCGGGTCTGAACCTGGATGGGAT  >  NC_002947/2010353‑2010651
                                                                                                                                                     |                                                                                                                                                     
tGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAgttgt                                                                                                                                                       <  2:1006893/150‑1 (MQ=255)
  gCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTc                                                                                                                                                     >  2:1915099/1‑150 (MQ=255)
   ccTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCa                                                                                                                                                    >  2:654979/1‑150 (MQ=255)
     ttGTGGCCCGTCGCGGCCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCaa                                                                                                                                                   <  2:1133816/149‑1 (MQ=255)
     ttGTGGCCCGTCGCGGCCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCaa                                                                                                                                                   >  1:1133816/1‑149 (MQ=255)
        tgGCCCGTCGCGGTCAGCATGGACGTCGGCCCAGTCGCAACCTTCCACCAGCAAATTGCGCGGACCTGTTGGCGCGCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGc                                                                                                                                               <  1:43406/150‑1 (MQ=255)
         ggCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCg                                                                                                                                              >  1:899818/1‑150 (MQ=255)
          gCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCgg                                                                                                                                             >  2:2331188/1‑150 (MQ=255)
               tCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGga                                                                                                                                          >  2:1630496/1‑148 (MQ=255)
               tCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGga                                                                                                                                          <  1:1630496/148‑1 (MQ=255)
                   ggTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGt                                                                                                                                    <  1:1000972/150‑1 (MQ=255)
                       acaGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGc                                                                                                                                <  2:1858885/150‑1 (MQ=255)
                         aGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCg                                                                                                                              >  2:86605/1‑150 (MQ=255)
                          gTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCgg                                                                                                                             >  2:439399/1‑150 (MQ=255)
                            cGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCgggg                                                                                                                           <  2:1054786/150‑1 (MQ=255)
                                aTAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCCGCGGTGGTTCAACTGCGGGGAGAAGGTCGGCCGGggcggc                                                                                                                       >  1:500756/1‑150 (MQ=255)
                                     cccAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCg                                                                                                                  <  1:654979/150‑1 (MQ=255)
                                       cAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGcc                                                                                                                >  1:1913677/1‑150 (MQ=255)
                                          gcgcAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGt                                                                                                             <  2:955674/150‑1 (MQ=255)
                                            gcAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTg                                                                                                           >  1:2237190/1‑150 (MQ=255)
                                            gcAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGATGCCAGCAATCGATCAAATCACTGCAGCAGATCTGCCGGGGCAGCCGTCGCCAGTTg                                                                                                           >  2:1453941/1‑150 (MQ=255)
                                              aaCCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGAt                                                                                                         <  2:1436509/150‑1 (MQ=255)
                                              aaCCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGAt                                                                                                         >  2:1683516/1‑150 (MQ=255)
                                              aaCCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACGCCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGAt                                                                                                         >  1:1204816/1‑150 (MQ=255)
                                               cccTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATg                                                                                                        <  2:563158/149‑1 (MQ=255)
                                               aCCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATg                                                                                                        <  1:1384267/150‑1 (MQ=255)
                                                   tcctccAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGaca                                                                                                    <  2:1850468/150‑1 (MQ=255)
                                                   tccgccAGCAAATCGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCCGCAGCGATCTGGACAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGaca                                                                                                    <  2:1300609/150‑1 (MQ=255)
                                                    cctccAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTTTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGCGCTTTCCGGCCGGGGCAGCCTTCGCCAGTTGATGGacac                                                                                                   <  1:789763/150‑1 (MQ=255)
                                                        cAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTc                                                                                               <  1:586235/150‑1 (MQ=255)
                                                          gCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAg                                                                                             >  2:419314/1‑150 (MQ=255)
                                                          gCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTCTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAg                                                                                             >  2:1875865/1‑150 (MQ=255)
                                                           cAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGc                                                                                            <  2:1787748/150‑1 (MQ=255)
                                                           cAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGc                                                                                            <  2:899818/150‑1 (MQ=255)
                                                           cAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGc                                                                                            <  1:1601787/150‑1 (MQ=255)
                                                                 gCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATc                                                                                      >  2:396172/1‑150 (MQ=255)
                                                                   gAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGc                                                                                    >  2:2219098/1‑150 (MQ=255)
                                                                   gAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGc                                                                                    >  1:573510/1‑150 (MQ=255)
                                                                      ggctcgcCATGCATCGGGGCGCAAGTCTGTCCTCTCCGGCGATATCGATAATGCTGTAGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCa                                                                                 <  2:500756/146‑1 (MQ=255)
                                                                           gcCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTa                                                                            >  2:148331/1‑150 (MQ=255)
                                                                              aTGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCg                                                                         >  1:1114141/1‑150 (MQ=255)
                                                                                  cTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGc                                                                     <  1:86605/150‑1 (MQ=255)
                                                                                     ggggCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGa                                                                  >  1:2431260/1‑150 (MQ=255)
                                                                                      gggCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGAt                                                                 <  2:2431260/150‑1 (MQ=255)
                                                                                       ggCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGAtt                                                                <  2:338268/150‑1 (MQ=255)
                                                                                           ccAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGAtttttt                                                            >  1:1246824/1‑150 (MQ=255)
                                                                                            cAGTTTCTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTc                                                           <  1:1875865/150‑1 (MQ=255)
                                                                                            aaGTTCGTCAACTGCAGCGATCTGGATAAGTCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTc                                                           <  2:1961091/149‑1 (MQ=255)
                                                                                              gTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTACCGAGGCTGATTTTTTCga                                                         <  1:419314/150‑1 (MQ=255)
                                                                                               tttGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCgag                                                        <  1:1683516/150‑1 (MQ=255)
                                                                                                ttGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCgaga                                                       >  1:1245692/1‑150 (MQ=255)
                                                                                                ttGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCgaga                                                       <  2:1814438/150‑1 (MQ=255)
                                                                                                  gTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCAAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGt                                                     <  1:2333377/150‑1 (MQ=255)
                                                                                                     aCCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCCGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGc                                                  <  1:148331/150‑1 (MQ=255)
                                                                                                       cTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGt                                                >  1:1266681/1‑150 (MQ=255)
                                                                                                        tGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTc                                               <  1:1915099/150‑1 (MQ=255)
                                                                                                          cAGCTATCTGGCTAATGCTGTAGCCGCCGGCTGAAAGCAAGGGTCCATCTGCGGGGGGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCAc                                             <  2:1204816/150‑1 (MQ=255)
                                                                                                                tCTGTATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGc                                                                                            <  1:2379230/97‑1 (MQ=255)
                                                                                                                tCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGc                                                                                            >  2:2379230/1‑97 (MQ=255)
                                                                                                                    gATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAg                                   <  2:1913677/150‑1 (MQ=255)
                                                                                                                     aTAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCTGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGc                                  <  2:1245692/150‑1 (MQ=255)
                                                                                                                      tAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCga                                                         <  2:628868/126‑1 (MQ=255)
                                                                                                                      tAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCga                                                         >  1:628868/1‑126 (MQ=255)
                                                                                                                          gCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTc                                                           <  1:1465299/120‑1 (MQ=255)
                                                                                                                          gCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTc                                                           >  2:1465299/1‑120 (MQ=255)
                                                                                                                          gCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTaa                             >  2:2421281/1‑150 (MQ=255)
                                                                                                                            tGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGAtttttt                                                            <  1:1604636/117‑1 (MQ=255)
                                                                                                                            ggTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGAtttttt                                                            >  2:1604636/2‑117 (MQ=255)
                                                                                                                              tCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATAGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAgg                         >  1:388268/1‑150 (MQ=255)
                                                                                                                                    ccGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAGGCCGGGt                   <  2:2237190/150‑1 (MQ=255)
                                                                                                                                      ggCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTa                                    <  1:1632770/131‑1 (MQ=255)
                                                                                                                                      ggCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTa                                    >  2:1632770/1‑131 (MQ=255)
                                                                                                                                      ggCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAGGCCGGGTCt                 <  1:2631/150‑1 (MQ=255)
                                                                                                                                         tGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCGAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAGGCCGGGTCTGaa              >  2:2280639/1‑150 (MQ=255)
                                                                                                                                         tGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAGGCCGGGTCTGaa              <  2:1246824/150‑1 (MQ=255)
                                                                                                                                           cccgcagTTGAACAACTGCTGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAgg                                                                               >  2:975010/4‑83 (MQ=38)
                                                                                                                                           cccgcagTTGAACAACTGCTGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAgg                                                                               <  1:975010/80‑1 (MQ=38)
                                                                                                                                           cccgcagTTGAACAACTGCTGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCa                                        <  1:1250911/119‑1 (MQ=255)
                                                                                                                                           cccgcagTTGAACAACTGCTGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCa                                        >  2:1250911/4‑122 (MQ=255)
                                                                                                                                           cccgcagTTGAACAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAGGCCGGGTCTGAAcc            >  1:2356986/4‑150 (MQ=255)
                                                                                                                                            cagcagTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACt                                      <  1:833837/123‑1 (MQ=255)
                                                                                                                                            cagcagTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACt                                      >  2:833837/1‑123 (MQ=255)
                                                                                                                                            cagcagTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAgg                         <  1:1409223/136‑1 (MQ=255)
                                                                                                                                            cagcagTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAgg                         >  2:1409223/1‑136 (MQ=255)
                                                                                                                                            cagcagTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAGGCCGGGTCTGAACCt           >  2:105063/1‑150 (MQ=255)
                                                                                                                                            cagcagTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAGGCCGGGTCTGAACCt           >  2:160222/1‑150 (MQ=255)
                                                                                                                                            cagcagTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAGGCCGGGTCTGAACCt           >  2:1339118/1‑150 (MQ=255)
                                                                                                                                             agcagTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCga                                                         >  2:1360055/1‑103 (MQ=255)
                                                                                                                                             agcagTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCga                                                         <  1:1360055/103‑1 (MQ=255)
                                                                                                                                                agTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAGGCCGGGTCTGAACCTggat       <  1:396172/150‑1 (MQ=255)
                                                                                                                                                    gttCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATAGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAGGCCGGGTCTGAACCTggatggga   <  2:388268/150‑1 (MQ=255)
                                                                                                                                                     ttCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGTGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAGGCCGGGTCTGAACCTggatgggat  >  2:914499/1‑150 (MQ=255)
                                                                                                                                                     ttCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAGGCCGGGTCTGAACCTggatgggat  >  1:1878140/1‑150 (MQ=255)
                                                                                                                                                     ttCAACTGCGGGGAGAAGGTCCGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTCCAGCGTCCTCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAAAAGGCCGGGTCTGAACCTggatgggat  <  2:1213577/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGCCTTGTGGCCCGTCGCGGTCACAGTCGGCATAAGCCCAAGCGCAACCTTCCTCCAGCAAATTGCGAAGCTGCGCCATGGCTCGGGGCGCCAGTTTGTCACCTGCAGCGATCTGGATAATGCTGTCGCAGCCGGCTGCCAGCAGTTGTTCAACTGCGGGGAGAAGGTCGGCCGGGGCAGCCTTCGCCAGTTGATGGACACTTTCAGCGTCATCGCCCAGGCTATCGAGGCTGATTTTTTCGAGAGTGGCGTCACCTTCACTTAGCAGTAACAGGCCGGGTCTGAACCTGGATGGGAT  >  NC_002947/2010353‑2010651

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: