Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,533,770 G→A 19.4% R594H (CGT→CAT)  secA → protein translocase subunit SecA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,533,7700GA19.4% 209.1 / 38.8 98R594H (CGT→CAT) secAprotein translocase subunit SecA
Reads supporting (aligned to +/- strand):  ref base G (39/40);  new base A (9/10);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

AACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGG  >  NC_002947/1533625‑1533917
                                                                                                                                                 |                                                                                                                                                   
aaCTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTAt                                                                                                                                                  >  5:298242/1‑149 (MQ=255)
 aCTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATc                                                                                                                                                 >  8:127315/1‑149 (MQ=255)
  cTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCg                                                                                                                                                <  5:212414/149‑1 (MQ=255)
  cTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCg                                                                                                                                                <  3:119756/149‑1 (MQ=255)
   tGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGa                                                                                                                                               >  4:99539/1‑149 (MQ=255)
     ggAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGa                                                                                                                                               <  6:368037/147‑1 (MQ=255)
       aaGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAAc                                                                                                                                           <  2:52873/149‑1 (MQ=255)
        aGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAAcc                                                                                                                                          <  7:127315/149‑1 (MQ=255)
         gCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAAc                                                                                                                                           >  6:237541/1‑147 (MQ=255)
            gAAGTGGCAGCCCTTGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCt                                                                                                                                      <  6:298242/149‑1 (MQ=255)
            gAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCt                                                                                                                                      >  1:149152/1‑149 (MQ=255)
                  gCAGCCCTGGATAACCCGACAGCCGAGCAGATCGCCGAGATCAAGGCCGATTTGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGTCCTGCATGTTATCTCTTCCTAGCGCCACGAATCGGGCCGTATCGACAACCAGCTGCGCgg                                                                                                                                <  5:207265/149‑1 (MQ=255)
                  gCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCgg                                                                                                                                <  3:404776/149‑1 (MQ=255)
                  ccAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCgg                                                                                                                                >  4:58987/2‑149 (MQ=255)
                      cccTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGc                                                                                                                            >  1:186903/1‑149 (MQ=255)
                         tGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTcc                                                                                                                         >  7:182904/1‑149 (MQ=255)
                           gAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCgg                                                                                                                       >  2:286146/1‑149 (MQ=255)
                             aaaCCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGACGCGTCACCAGCAAGTGAGCGAGACGGGTGGCCTGCATGTAATCGCTTCCGCGCGCCACGAATCGCGCCGTATCGACAACCTGCTGCGCGGCCGCTCCGgtc                                                                                                                     <  5:427997/149‑1 (MQ=255)
                              aaCCCGACAGCCGAGCAGATAGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGgtcg                                                                                                                    >  6:443545/1‑149 (MQ=255)
                                cccGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGgtcgt                                                                                                                   >  8:141625/1‑148 (MQ=255)
                                cccGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGgtcgt                                                                                                                   >  8:20937/1‑148 (MQ=255)
                                 ccGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCa                                                                                                                 <  7:141625/149‑1 (MQ=255)
                                   gACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAgg                                                                                                               <  8:194512/149‑1 (MQ=255)
                                   gACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAgg                                                                                                               >  1:92073/1‑149 (MQ=255)
                                    aCAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAggg                                                                                                              <  3:99539/149‑1 (MQ=255)
                                        ccGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGAc                                                                                                          >  1:390110/1‑149 (MQ=255)
                                         cGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGAcc                                                                                                         <  7:20937/149‑1 (MQ=255)
                                           aGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCg                                                                                                       >  8:288816/1‑149 (MQ=255)
                                            gCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCgg                                                                                                      <  1:371043/149‑1 (MQ=255)
                                               gATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGtt                                                                                                   >  6:241523/1‑149 (MQ=255)
                                                 tCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTcc                                                                                                 <  4:197374/149‑1 (MQ=255)
                                                   gCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCa                                                                                                <  2:149152/148‑1 (MQ=255)
                                                    cccAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGc                                                                                              <  7:354627/149‑1 (MQ=255)
                                                    cccAGATCAAGGCAGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGc                                                                                              >  6:393773/1‑149 (MQ=255)
                                                       aGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGc                                                                                           <  6:201259/149‑1 (MQ=255)
                                                            aaGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGCCAGGGCGACCCGGGTTCCAGCCGCTTCTa                                                                                      <  8:71648/149‑1 (MQ=255)
                                                              ggCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTAcc                                                                                    >  1:363280/1‑149 (MQ=255)
                                                              ggCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCggg                                                                                                     <  4:169456/132‑1 (MQ=255)
                                                              ggCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCggg                                                                                                     >  3:169456/1‑132 (MQ=255)
                                                               gCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCt                                                                                   >  6:382464/1‑149 (MQ=255)
                                                                   aTTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCACCCGCTTCTACCTCTCg                                                                               >  4:279700/1‑149 (MQ=255)
                                                                      ggCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCt                                                                             <  5:33404/148‑1 (MQ=255)
                                                                      ggCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCt                                                                             <  5:312733/148‑1 (MQ=255)
                                                                       gCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTgg                                                                           <  8:281677/149‑1 (MQ=255)
                                                                       gCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTgg                                                                           <  3:303644/149‑1 (MQ=255)
                                                                       gCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTg                                                                            >  5:303502/1‑148 (MQ=255)
                                                                       gCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTgg                                                                           <  2:92073/149‑1 (MQ=255)
                                                                        cAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTgg                                                                           <  5:443545/148‑1 (MQ=255)
                                                                         aGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGaa                                                                         <  7:288816/149‑1 (MQ=255)
                                                                               gTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGc                                                                   <  1:286146/149‑1 (MQ=255)
                                                                                tCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGCAGACAGcc                                                                  >  4:60252/1‑149 (MQ=255)
                                                                                tCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCgg                                                                                                      >  7:81328/1‑113 (MQ=255)
                                                                                tCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCgg                                                                                                      <  8:81328/113‑1 (MQ=255)
                                                                                    cAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGAt                                                              <  5:393773/149‑1 (MQ=255)
                                                                                     aGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATg                                                             <  6:98951/149‑1 (MQ=255)
                                                                                        aaGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGGCCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATgcg                                                           >  6:289311/1‑148 (MQ=255)
                                                                                        aaGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATgcgc                                                          <  2:363280/149‑1 (MQ=255)
                                                                                          gTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCAt                                                        <  2:186903/149‑1 (MQ=255)
                                                                                             aTCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCtt                                                     >  3:2990/1‑149 (MQ=255)
                                                                                              tCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTc                                                    <  5:382464/149‑1 (MQ=255)
                                                                                              tCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTc                                                    >  2:312807/1‑149 (MQ=255)
                                                                                                gagaCGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGc                                                  >  8:184776/1‑149 (MQ=255)
                                                                                                 agccGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGc                                                  >  7:37359/4‑148 (MQ=255)
                                                                                                  gaCGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCct                                                <  4:169679/149‑1 (MQ=255)
                                                                                                      ggTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTc                                                    >  7:267031/1‑141 (MQ=255)
                                                                                                      ggTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTc                                                    <  8:267031/141‑1 (MQ=255)
                                                                                                       gTGGCCTGCATGTATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGAcc                                          >  4:172600/1‑149 (MQ=255)
                                                                                                          gCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCgg                                        <  5:237541/149‑1 (MQ=255)
                                                                                                            cTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGt                                      >  5:317903/1‑149 (MQ=255)
                                                                                                            cTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGt                                      >  2:346826/1‑149 (MQ=255)
                                                                                                            cTGCATGTAATCGCTTCCGAGCGCCACGAATCGCACCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGt                                      >  1:64653/1‑149 (MQ=255)
                                                                                                             tGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGt                                      >  7:303723/1‑148 (MQ=255)
                                                                                                               cATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTga                                    <  8:182904/148‑1 (MQ=255)
                                                                                                               cATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTg                                     >  6:407700/1‑147 (MQ=255)
                                                                                                                 tGTAATCGCTTGCGAGCGCCACGAATCGCGCCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTgaaga                                 <  2:390110/149‑1 (MQ=255)
                                                                                                                 tGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATg                                                             >  2:194306/1‑121 (MQ=255)
                                                                                                                 tGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATg                                                             <  1:194306/121‑1 (MQ=255)
                                                                                                                  gTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTgaagaa                                >  7:96451/1‑149 (MQ=255)
                                                                                                                  gTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTgaagaa                                >  1:106479/1‑149 (MQ=255)
                                                                                                                     aTCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACt                              >  6:211739/1‑148 (MQ=255)
                                                                                                                      tCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTGCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTGTGACCGGGTGAAGAACTTc                            <  6:329462/149‑1 (MQ=255)
                                                                                                                      tCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGCAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTc                            <  6:303502/149‑1 (MQ=255)
                                                                                                                         cTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCAt                          >  6:375969/1‑148 (MQ=255)
                                                                                                                          ttCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATg                         >  2:368282/1‑148 (MQ=255)
                                                                                                                            ccGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGaa                       <  5:241523/148‑1 (MQ=255)
                                                                                                                                   ccACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCGGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGAGTGAAGAACTTCATGAAGGCACTgg               >  3:291406/1‑149 (MQ=255)
                                                                                                                                   ccACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTgg               <  7:184776/149‑1 (MQ=255)
                                                                                                                                    cACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGt                                                                                                    >  4:245007/1‑63 (MQ=255)
                                                                                                                                    cACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGt                                                                                                    <  3:245007/63‑1 (MQ=255)
                                                                                                                                    cACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTggg              >  2:415380/1‑149 (MQ=255)
                                                                                                                                    cACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTggg              >  7:150086/1‑149 (MQ=255)
                                                                                                                                     aCGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGc             <  6:317903/149‑1 (MQ=255)
                                                                                                                                       gAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCTTCAGGGCGACCCGGGTTCCATCCGCTTCTACCTGTCGCTTGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCAt           <  8:4222/149‑1 (MQ=255)
                                                                                                                                        aaTCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCGGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGAGTGAAGAACTTCATGAAGGCACTGGGCATg          <  4:291406/149‑1 (MQ=255)
                                                                                                                                            gcgcCATATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGCGCGACCAGGGTTCCAGCCGCTTCTACCTGTCGCACGAAGACAGACTGATGCTCATCTTCGCCTCTGACCGGGAGAACAACTTCATGAAGGCACTGGGCATGCAga      >  8:243467/1‑148 (MQ=255)
                                                                                                                                               ccGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGCGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCgg   <  2:106479/149‑1 (MQ=255)
                                                                                                                                               ccGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTGCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCgg   <  7:187594/149‑1 (MQ=255)
                                                                                                                                                cGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCggg  >  3:1255/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AACTGGGAAGCCGAAGTGGCAGCCCTGGAAAACCCGACAGCCGAGCAGATCGCCCAGATCAAGGCCGATTGGCAGAAGCGTCACCAGCAAGTGATCGAGACGGGTGGCCTGCATGTAATCGCTTCCGAGCGCCACGAATCGCGCCGTATCGACAACCAGCTGCGCGGCCGCTCCGGTCGTCAGGGCGACCCGGGTTCCAGCCGCTTCTACCTGTCGCTGGAAGACAGCCTGATGCGCATCTTCGCCTCTGACCGGGTGAAGAACTTCATGAAGGCACTGGGCATGCAGTCGGG  >  NC_002947/1533625‑1533917

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