Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,533,185 C→G 5.1% T399S (ACC→AGC) ‡ secA → protein translocase subunit SecA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,533,1850CG 5.1% 257.3 / 2.9 79T399S (ACC→AGC) ‡secAprotein translocase subunit SecA
Reads supporting (aligned to +/- strand):  ref base C (37/38);  new base G (2/2);  total (39/40)
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.96e-01

CGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATCTACGGCCTGAACGTGATGGTCATTCCGCCGAACAAGCCACTGGCGCGT  >  NC_002947/1533037‑1533333
                                                                                                                                                    |                                                                                                                                                     
cGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGAc                                                                                                                                                       <  6:216748/149‑1 (MQ=255)
   aaCTTCGAGTACATCGTCCAGGACGGCTAGGTTCTGCTGATCGACGATCACACCTGCCGTAGCATGACGGGCCGTCGCGTGTCCGAGGGCGTGCACCAGGCCATCTAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCt                                                                                                                                                    <  2:120823/149‑1 (MQ=255)
    aCGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTg                                                                                                                                                   <  6:20053/149‑1 (MQ=255)
    aCGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTg                                                                                                                                                   <  2:347408/149‑1 (MQ=255)
       tCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGGGTCTGGCt                                                                                                                                                >  3:352644/1‑149 (MQ=255)
         gAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGGGTCTGGCtc                                                                                                                                               <  4:352644/148‑2 (MQ=255)
           gTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGTAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGa                                                                                                                                            <  1:100820/149‑1 (MQ=255)
             aCATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGac                                                                                                                                           >  3:120712/1‑148 (MQ=255)
              cATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTGCGAGGGCCTGCACCAGGCCATCGAAGCAAAAGAGAACCTGAACATCCAGGCCGTGAGCCAGACCCTGGCTTCGacc                                                                                                                                          <  2:184187/148‑1 (MQ=255)
                tCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGaccacc                                                                                                                                       >  7:127891/1‑149 (MQ=255)
                  gTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCtt                                                                                                                                    >  4:76630/1‑149 (MQ=255)
                   tCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGGGTCTGGCtctcg                                                                                                                                            >  1:353318/1‑137 (MQ=255)
                   tCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGGGTCTGGCtctcg                                                                                                                                            <  2:353318/141‑5 (MQ=255)
                   tCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTc                                                                                                                                   <  6:3750/149‑1 (MQ=255)
                     cAGGACGGCCAGGTTCTGCTGAGCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTcc                                                                                                                                  <  7:381732/148‑1 (MQ=255)
                        gACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGaa                                                                                                                              <  5:182594/149‑1 (MQ=255)
                        gACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGaa                                                                                                                              >  3:185918/1‑149 (MQ=255)
                             ccAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAactact                                                                                                                         <  3:91587/149‑1 (MQ=255)
                             ccAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAactact                                                                                                                         >  6:91153/1‑149 (MQ=255)
                               aGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGACCCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTc                                                                                                                       <  6:166962/149‑1 (MQ=255)
                               aGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTc                                                                                                                       <  3:43610/149‑1 (MQ=255)
                               aGGTGCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTc                                                                                                                       <  2:271880/149‑1 (MQ=255)
                                 gTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCg                                                                                                                     >  7:124198/1‑149 (MQ=255)
                                  ttCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGc                                                                                                                    <  6:246241/149‑1 (MQ=255)
                                        tgATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTa                                                                                                               >  6:188573/1‑148 (MQ=255)
                                        tgATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTa                                                                                                               >  4:103763/1‑148 (MQ=255)
                                            cgacgaACACACCGGCCGTACCATGCCGCGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACAcc                                                                                                           >  1:230553/1‑148 (MQ=255)
                                             gacgaACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGc                                                                                                                    >  5:347486/1‑138 (MQ=255)
                                             gacgaACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGc                                                                                                                    <  6:347486/138‑1 (MQ=255)
                                              acgaACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAg                                                                                                        <  8:124198/149‑1 (MQ=255)
                                              acgaACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTACGCCTGTACACCAAg                                                                                                        >  7:324892/1‑149 (MQ=255)
                                                    acacCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCt                                                                                                                  >  6:367432/1‑133 (MQ=255)
                                                    acacCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCt                                                                                                                  <  5:367432/133‑1 (MQ=255)
                                                      acCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGa                                                                                                                                            >  6:373620/1‑106 (MQ=255)
                                                      acCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGa                                                                                                                                            <  5:373620/106‑1 (MQ=255)
                                                      acCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCgg                                                                                                >  3:23707/1‑149 (MQ=255)
                                                       ccGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGc                                                                                               <  3:103763/149‑1 (MQ=255)
                                                                ccATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCgg                                                                                       >  5:88002/1‑148 (MQ=255)
                                                                 cATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTa                                                                                     >  6:182104/1‑149 (MQ=255)
                                                                         gCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCg                                                                                                 >  6:101345/1‑129 (MQ=255)
                                                                         gCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCg                                                                                                 <  5:101345/129‑1 (MQ=255)
                                                                          ccGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGaca                                                                            <  4:23707/149‑1 (MQ=255)
                                                                          ccGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGaca                                                                            >  5:132387/1‑149 (MQ=255)
                                                                           cGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGacac                                                                           <  3:76630/149‑1 (MQ=255)
                                                                             tagCCTGTCCGTTGTGCTGCTCCAGGCCATCGAAGCCAAACAGAACCTGAACATCAAGGCCGAGATCCAGACCCTGGCTTCGACCACCATCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGAATGACGGGTACCGCCGACACCg                                                                         <  8:324892/147‑1 (MQ=255)
                                                                             tCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACAcc                                                                          >  3:274128/1‑148 (MQ=255)
                                                                               gCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCg                                                                         >  5:211008/1‑147 (MQ=255)
                                                                                ccTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAg                                                                      <  7:73399/149‑1 (MQ=255)
                                                                                ccTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAg                                                                      <  6:211008/149‑1 (MQ=255)
                                                                                ccTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAg                                                                      <  6:88002/149‑1 (MQ=255)
                                                                                   gTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTa                                                                                                               >  3:105131/1‑105 (MQ=255)
                                                                                   gTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTa                                                                                                               >  3:110713/1‑105 (MQ=255)
                                                                                   gTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTa                                                                                                               <  4:105131/105‑1 (MQ=255)
                                                                                   gTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTa                                                                                                               <  4:110713/105‑1 (MQ=255)
                                                                                    ttcGATGGCCTGGGGCAGGCCATCGAAGCAAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACACTGGCTTCGACCACCTTCCAGAACTACATCCGCCAGTACACCAAACTCTCCGGCAAGACCGGTACCGCCGCCACCGAAGCAtt                                                                  >  8:353318/3‑149 (MQ=255)
                                                                                    tCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCAtt                                                                  >  3:151858/1‑149 (MQ=255)
                                                                                      cGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTc                                                                 <  4:151858/148‑1 (MQ=255)
                                                                                      cGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTc                                                                 >  2:68589/1‑148 (MQ=255)
                                                                                        aGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGaa                                                                       <  5:300106/140‑1 (MQ=255)
                                                                                        aGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGaa                                                                       >  6:300106/1‑140 (MQ=255)
                                                                                           gCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTc                                                           <  5:188573/149‑1 (MQ=255)
                                                                                           gCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTc                                                           <  4:185918/149‑1 (MQ=255)
                                                                                            ccTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTc                                                           >  8:121072/1‑148 (MQ=255)
                                                                                             cTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCGGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTc                                                           >  1:267547/1‑147 (MQ=255)
                                                                                                cACCAGGCCAGCGAAGCCAAAGTGAACCTGAAGATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCGTGTACACCAAGCTGTCCGGCATGACCGTTACCGCCGACACCGATGCATTCGAGTTCCAg                                                        <  4:119347/147‑1 (MQ=255)
                                                                                                                      agaACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATCTACGGCCTGAACGTGATg                                <  4:120712/149‑1 (MQ=255)
                                                                                                                        aaCCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATCTACGGCCTGAACGTGATGGt                              <  2:230553/149‑1 (MQ=255)
                                                                                                                          ccTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCGTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATCTACGGCCTGAACGTGATGGt                              >  1:77234/1‑148 (MQ=255)
                                                                                                                             gAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATCTACGGCCTGAACGTGATGGTCATTc                         >  1:152925/1‑149 (MQ=255)
                                                                                                                               aCATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCACGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATATACTGCCTGAACGTGATGGTCATTcc                        <  8:127891/148‑1 (MQ=255)
                                                                                                                                 aTCCATGCCGAGAGCCAGACCCTGGCTTCGACCAGCTTCGAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATCTACTGCCTGAACGTGATGGTCATTccgcc                     <  7:353318/149‑1 (MQ=255)
                                                                                                                                 aTCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATCTACGGCCTGAACGTGATGGTCATTccgcc                     >  4:270419/1‑149 (MQ=255)
                                                                                                                                              gCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACATCGAAGCATTCGAGTTCCAGTCGATCTACGGCCTGAACGTGATGGTCATTCCGCCGAACAAGCCACTg        >  4:120059/1‑149 (MQ=255)
                                                                                                                                              gCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATCTACGGCCTGAACGTGTTGGTCATTCCGCCGAACAAGCCACTg        >  3:264739/1‑149 (MQ=255)
                                                                                                                                               ccAGACCCTGGCTTCGACCACCTTCCTGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATCTACGGCCTGAACGTGATGGTCATTCCGCCGAACAAGCCACTgg       >  1:218316/1‑149 (MQ=255)
                                                                                                                                               ccAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATCTACGGCCTGAACGTGATGGTCATTCCGCCGAACAAGCCACTgg       >  3:391666/1‑149 (MQ=255)
                                                                                                                                               ccAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATCTACGGCCTGAACGTGATGGTCATTCCGCCGAACAAGCCACTgg       >  2:298718/1‑149 (MQ=255)
                                                                                                                                                cAGCCCCTGCCTTCGACCACCTTCCAGAACTACTTCCGCCTCTACACCAAGCTGTCCGGCATGACCCGTACCGCCGACACCCAAGCATTCGACTTCCAGACGATCTACGCCCTCAACCTGATGCTCAATCCGCCGAACAACCCACtctc      >  4:271146/1‑146 (MQ=255)
                                                                                                                                                  gACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATCTACGGCCTGAACGTGATGGTCATTCCGCCGAACAAGCCACTGgcgc    <  2:267547/149‑1 (MQ=255)
                                                                                                                                                  gACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATCTACGGCCTGAACGTGATGGTCATTCCGCCGAACAAGCCACTGgcgc    <  2:152925/149‑1 (MQ=255)
                                                                                                                                                    cccTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATCTACGGCCTGAACGTGATGGTCATTCCGCCGAACAAGCCACTGGCGCGt  <  7:267936/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGCAACGTCGAGTACATCGTCCAGGACGGCCAGGTTCTGCTGATCGACGAACACACCGGCCGTACCATGCCGGGCCGTCGCCTGTCCGAGGGCCTGCACCAGGCCATCGAAGCCAAAGAGAACCTGAACATCCAGGCCGAGAGCCAGACCCTGGCTTCGACCACCTTCCAGAACTACTTCCGCCTGTACACCAAGCTGTCCGGCATGACCGGTACCGCCGACACCGAAGCATTCGAGTTCCAGTCGATCTACGGCCTGAACGTGATGGTCATTCCGCCGAACAAGCCACTGGCGCGT  >  NC_002947/1533037‑1533333

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