Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,943,063 T→A 6.1% K292* (AAG→TAG)  PP_3478 ← secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,943,0630TA 6.1% 269.3 / 4.9 82K292* (AAG→TAG) PP_3478secretion protein
Reads supporting (aligned to +/- strand):  ref base T (37/40);  new base A (3/2);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTGGCGCAGCTGGTTTTGC  >  NC_002947/3942917‑3943211
                                                                                                                                                  |                                                                                                                                                    
cgccACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGAc                                                                                                                                                   >  2:1243914/3‑150 (MQ=255)
cTCCACTTCCCGCACCACCTCGGGGTTGGACTGGTCCGGCGACTGAAGCAGCTTCCCGGCCCTGCCCACCGCGTCCGGCTTGTCGCGCAGGGGCAGGGTGTTTAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGAc                                                                                                                                                   <  2:2079205/150‑1 (MQ=255)
 tCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGAcc                                                                                                                                                  >  1:1193218/1‑150 (MQ=255)
  ccACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTACGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCa                                                                                                                                                 >  1:869209/1‑150 (MQ=255)
   cACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCTGCTGCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGCAACGGCATTCAAAGGCAAATCAACAAACGCGTCTCGGTTCTTGAGCAGGGTCTAGatcgg                                                                                                                                                >  2:2329337/1‑146 (MQ=255)
          cgCACCCCCTCGGAGTTGGACTTGTCCTGCGACTGCAGCAGCTCCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAACATGGTCAAGACCATGCTCaag                                                                                                                                         <  1:564995/149‑2 (MQ=255)
          aGCACCACCTCGGGGTTGGACTTGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAAc                                                                                                                                         <  1:1243914/150‑1 (MQ=255)
              ccacgTCGGGGTTGGACTCGTCCCGCTACTGCAGCAGGTTCTCGGACATGCGCACCGCGTCCGACGTGTCGCGTATGGTCAGGGTCTTGAGGCTATAGTCAACGAACACGTCTCGGGTCTTGAGCATGGTCCTGACCATGTTCAACGCgg                                                                                                                                     <  2:114758/150‑1 (MQ=255)
                   tCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTg                                                                                                                                <  2:798573/150‑1 (MQ=255)
                   tCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTg                                                                                                                                <  2:2171185/150‑1 (MQ=255)
                       ggTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGt                                                                                                                            <  1:700961/150‑1 (MQ=255)
                           ggACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGAt                                                                                                                        <  2:856848/150‑1 (MQ=255)
                            gACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTAAGTAGTTTGTTGACCATGTTCAACGCGGTGTTGGCGTCGATg                                                                                                                       >  1:22193/1‑150 (MQ=255)
                               tGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTg                                                                                                                    <  1:2190512/150‑1 (MQ=255)
                                  tCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTc                                                                                                                 >  1:1113776/1‑150 (MQ=255)
                                    cTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCaa                                                                                                               >  1:960339/1‑149 (MQ=255)
                                    cTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAg                                                                                                               >  2:1335781/1‑150 (MQ=255)
                                    cTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAg                                                                                                               >  2:660418/1‑150 (MQ=255)
                                      gCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGt                                                                                                             <  1:660418/150‑1 (MQ=255)
                                       cGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATg                                                                                                                       <  2:1988813/139‑1 (MQ=255)
                                       cGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATg                                                                                                                       >  1:1988813/1‑139 (MQ=255)
                                             cagcagCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGt                                                                                                      <  1:1992770/150‑1 (MQ=255)
                                                  gcTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCa                                                                                                 <  1:968206/150‑1 (MQ=255)
                                                   cTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCAt                                                                                                >  1:233524/1‑150 (MQ=255)
                                                    ttCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATg                                                                                               >  2:2071295/1‑150 (MQ=255)
                                                      ctcGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGCGTGTAGAACCGCGCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGAc                                                                                             <  1:1740574/150‑1 (MQ=255)
                                                       tcGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGtt                                                                                                                                 <  2:1285540/113‑1 (MQ=255)
                                                       tcGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGtt                                                                                                                                 >  1:1285540/1‑113 (MQ=255)
                                                         ggCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAg                                                                                          >  2:1894470/1‑150 (MQ=255)
                                                         ggCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAg                                                                                          >  1:2250114/1‑150 (MQ=255)
                                                         ggCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGATCAGGTAGAACATGCAAATGACCAg                                                                                          >  2:1942572/1‑150 (MQ=255)
                                                         ggCCATGCGCACCGCGTACGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAg                                                                                          <  2:869209/150‑1 (MQ=255)
                                                         ggCCATGCCCATCGCGACCGACGTGTCGCGCATGGTCAGGTTGTTGAGGCGCTCGTAAACGAACACGTCGCGGGTCTTGTGCATGGTCTTGACCATCTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGACCGTGCGCATGACCAg                                                                                          <  2:744316/150‑1 (MQ=255)
                                                          gCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAAcgc                                                                                                                                       >  2:915630/1‑104 (MQ=255)
                                                          gCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAAcgc                                                                                                                                       <  1:915630/104‑1 (MQ=255)
                                                           ccATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGtt                                                                                        <  2:130014/150‑1 (MQ=255)
                                                           ccATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGtt                                                                                        >  2:2257315/1‑150 (MQ=255)
                                                           ccATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGtt                                                                                        >  2:929175/1‑150 (MQ=255)
                                                                cgcACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTgg                                                                                   <  1:1942572/150‑1 (MQ=255)
                                                                   aCCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAg                                                                                >  1:697278/1‑150 (MQ=255)
                                                                     cgcgTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTg                                                                                                                                >  2:383166/1‑100 (MQ=255)
                                                                     cgcgTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTg                                                                                                                                <  1:383166/100‑1 (MQ=255)
                                                                       cgTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCt                                                                            <  1:2071295/150‑1 (MQ=255)
                                                                       cgTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCt                                                                            <  1:1821546/150‑1 (MQ=255)
                                                                         tCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTg                                                                          <  2:233524/150‑1 (MQ=255)
                                                                         tCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTg                                                                          >  2:651271/1‑150 (MQ=255)
                                                                                 gtCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTg                                                                  >  1:2303254/1‑150 (MQ=255)
                                                                                     cgcATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTgcggc                                                              <  2:195789/150‑1 (MQ=255)
                                                                                      gcATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGtt                                                                                                                                             >  2:2287376/1‑70 (MQ=255)
                                                                                      gcATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGtt                                                                                                                                             <  1:2287376/70‑1 (MQ=255)
                                                                                      ccATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTgcggcg                                                             >  2:1657746/2‑150 (MQ=255)
                                                                                          ggTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATc                                                         >  1:280453/1‑150 (MQ=255)
                                                                                           gTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGAGTTGTCTTCGATGCTGGTTAGGTAGAACGTGCGCATGATCGGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCg                                                        >  1:160568/1‑150 (MQ=255)
                                                                                            tCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTg                                                                                    >  2:334781/1‑121 (MQ=255)
                                                                                            tCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTg                                                                                    <  1:334781/121‑1 (MQ=255)
                                                                                            tCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCgg                                                       <  2:2303254/150‑1 (MQ=255)
                                                                                             cAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCggg                                                      <  1:929175/150‑1 (MQ=255)
                                                                                             cAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCggg                                                      <  1:1657746/150‑1 (MQ=255)
                                                                                                  tgtTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGa                                                 <  1:2398045/150‑1 (MQ=255)
                                                                                                    ttGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCggg                                                      >  1:2041313/1‑143 (MQ=255)
                                                                                                    ttGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCggg                                                      <  2:2041313/143‑1 (MQ=255)
                                                                                                      gAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAg                                                                                                               <  2:1588605/84‑1 (MQ=255)
                                                                                                      gAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAg                                                                                                               >  1:1588605/1‑84 (MQ=255)
                                                                                                             tCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGgtg                                      <  1:1894470/150‑1 (MQ=255)
                                                                                                               gTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCggg                                                      <  1:2264128/132‑1 (MQ=255)
                                                                                                               gTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCggg                                                      >  2:2264128/1‑132 (MQ=255)
                                                                                                                  aaCGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTc                                 <  1:488965/150‑1 (MQ=255)
                                                                                                                   aCGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCg                                >  1:1172643/1‑150 (MQ=255)
                                                                                                                     gAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGtt                              <  2:1172643/150‑1 (MQ=255)
                                                                                                                     gAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGAGGCGAATCGGGGTAGACCATGACGGTGTTGTCGtt                              >  2:562437/1‑150 (MQ=255)
                                                                                                                       acacGTCGCGGGTCTAGAGGAGGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGa                            >  2:1193061/1‑150 (MQ=255)
                                                                                                                         acGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACg                          <  2:1113776/150‑1 (MQ=255)
                                                                                                                            tCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGAtc                       <  1:1196768/150‑1 (MQ=255)
                                                                                                                              gcgGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGAtctt                     >  1:758015/1‑150 (MQ=255)
                                                                                                                                   tCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTGgc                >  2:1445935/1‑150 (MQ=255)
                                                                                                                                     ttGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCa                                                                                                                >  2:885255/1‑52 (MQ=255)
                                                                                                                                     ttGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCa                                                                                                                <  1:885255/52‑1 (MQ=255)
                                                                                                                                     ttGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTGgcgc              >  2:13579/1‑150 (MQ=255)
                                                                                                                                       gaacATGGTCAAGACCATGCTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGAtc                       >  2:2361547/4‑139 (MQ=255)
                                                                                                                                       gaacATGGTCAAGACCATGCTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGAtc                       <  1:2361547/136‑1 (MQ=255)
                                                                                                                                        aacATGGTCAAGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGCGCGCATGACCAGTTCCTGGTAGTCCTGGGTCCTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTGCTTGACGATCTTCTCGCGCAGc           >  2:541492/3‑150 (MQ=255)
                                                                                                                                             ggTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTGGCGCAGCTGGtt      <  1:2151399/150‑1 (MQ=255)
                                                                                                                                                 ttGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTGGCGCAGCTGGTTTTGc  <  1:1335781/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTCCACTTCCAGCACCACCTCGGGGTTGGACTGGTCCTGCGACTGCAGCAGCTTCTCGGCCATGCGCACCGCGTCCGACGTGTCGCGCATGGTCAGGGTGTTGAGGCGCTCGTCAACGAACACGTCGCGGGTCTTGAGCATGGTCTTGACCATGTTCAACGCGGTGTTGGCGTCGATGCTGGTCAGGTAGAACGTGCGCATGACCAGTTCCTGGTAGTCCTTGGTCTTCTGCGGCGAATCGGGGTAGACCATGACGGTGTTGTCGTTGACGATCTTCTGGCGCAGCTGGTTTTGC  >  NC_002947/3942917‑3943211

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