Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,772,157 A→C 5.6% I269S (ATC→AGC)  betT‑I ← choline transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,772,1570AC 5.6% 225.9 / 4.2 72I269S (ATC→AGC) betT‑Icholine transporter
Reads supporting (aligned to +/- strand):  ref base A (35/33);  new base C (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

CCGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACTT  >  NC_002947/5772010‑5772304
                                                                                                                                                   |                                                                                                                                                   
ccGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCTCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCatgat                                                                                                                                                    >  7:179523/1‑149 (MQ=255)
     ttCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGTACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGtt                                                                                                                                                 >  1:201146/1‑147 (MQ=255)
       cGCCTTCCTACACATAGAGGGCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAg                                                                                                                                              <  6:409864/148‑1 (MQ=255)
       cGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCAACTGAACATCATGCTGTTCAg                                                                                                                                              <  1:155474/148‑1 (MQ=255)
        gCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACATCATGCTGTTCAgtg                                                                                                                                            <  6:360533/149‑3 (MQ=255)
         ccTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGtt                                                                                                                                           <  1:157231/149‑1 (MQ=255)
              aTACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGa                                                                                                                                        >  6:150212/1‑147 (MQ=255)
               tACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAg                                                                                                                                     >  2:40003/1‑149 (MQ=255)
                acacATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAgg                                                                                                                                    >  1:165158/1‑149 (MQ=255)
                 cacaTAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAggc                                                                                                                                   <  4:227012/149‑1 (MQ=255)
                  acaTAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGtt                                                                                                                                           >  1:173914/1‑140 (MQ=255)
                    aTAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAggcggc                                                                                                                                >  1:15291/1‑149 (MQ=255)
                    aTAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAggcggc                                                                                                                                >  6:101672/1‑149 (MQ=255)
                          gTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGCTCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGc                                                                                                                          >  1:311746/1‑149 (MQ=255)
                          gTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGc                                                                                                                          >  1:427116/1‑149 (MQ=255)
                           tCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGcc                                                                                                                         <  5:101672/149‑1 (MQ=255)
                           tCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGcc                                                                                                                         <  8:147001/149‑1 (MQ=255)
                             gAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCg                                                                                                                        <  2:427116/148‑1 (MQ=255)
                                  ttttCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTtctc                                                                                                                   >  7:187001/1‑148 (MQ=255)
                                  ttttCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTcc                                                                                                                  <  1:22424/149‑1 (MQ=255)
                                       aGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGcc                                                                                                             >  2:386356/1‑149 (MQ=255)
                                              tGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATAcc                                                                                                      >  6:199809/1‑149 (MQ=255)
                                                 cGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCt                                                                                                   >  2:328433/1‑149 (MQ=255)
                                                   ttCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGc                                                                                                 <  2:201146/149‑1 (MQ=255)
                                                   ttCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGc                                                                                                 <  5:249326/149‑1 (MQ=255)
                                                        gTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACgg                                                                                            >  4:327746/1‑149 (MQ=255)
                                                         tAGTCGCCAACGTGCTGCACGACGCCGTTGAGCAGGTGCAGGGGTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGt                                                                                           <  1:356046/149‑1 (MQ=255)
                                                           gTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTa                                                                                          <  3:50204/148‑1 (MQ=255)
                                                           gTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTa                                                                                          <  5:150212/148‑1 (MQ=255)
                                                              gCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCa                                                                                      >  2:223738/1‑149 (MQ=255)
                                                                 aaCGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGc                                                                                   <  2:165158/149‑1 (MQ=255)
                                                                       cTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGc                                                                                        <  2:35891/138‑1 (MQ=255)
                                                                       cTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGc                                                                                        >  1:35891/1‑138 (MQ=255)
                                                                          cACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCAt                                                                            <  8:179523/147‑1 (MQ=255)
                                                                          cACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGa                                                                          >  5:64153/1‑149 (MQ=255)
                                                                             gAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGAc                                                                         >  1:380907/1‑147 (MQ=255)
                                                                              aaGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAg                                                                      >  7:58554/1‑149 (MQ=255)
                                                                                 ccGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGAt                                                                   <  6:204541/149‑1 (MQ=255)
                                                                                   gTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATca                                                                 >  3:170010/1‑149 (MQ=255)
                                                                                       aGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATcaccac                                                             >  3:806/1‑149 (MQ=255)
                                                                                               cAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGAc                                                                         >  7:8396/1‑129 (MQ=255)
                                                                                               cAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGAc                                                                         <  8:8396/129‑1 (MQ=255)
                                                                                                   gTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGTAGAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACGACCAGCAGGTTGGt                                                 <  2:403241/149‑1 (MQ=255)
                                                                                                   gTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGt                                                 >  2:75641/1‑149 (MQ=255)
                                                                                                    tttCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGa                                                                    >  7:146831/1‑129 (MQ=255)
                                                                                                    tttCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGa                                                                    <  8:146831/129‑1 (MQ=255)
                                                                                                    tttCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTc                                                >  8:36047/1‑149 (MQ=255)
                                                                                                            cccAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACTGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCtgtg                                        <  2:15291/149‑1 (MQ=255)
                                                                                                             ccAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTgg                                       >  6:279189/1‑149 (MQ=255)
                                                                                                                 cagcaCGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCa                                   >  5:292938/1‑149 (MQ=255)
                                                                                                                  agcaCGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCAt                                  <  1:386356/149‑1 (MQ=255)
                                                                                                                   gcaCGAAGCTCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATg                                 <  2:311746/149‑1 (MQ=255)
                                                                                                                   gcaCGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATg                                 <  8:58554/149‑1 (MQ=255)
                                                                                                                    caCGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGc                                >  5:177185/1‑149 (MQ=255)
                                                                                                                     aCGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGcc                               <  1:75641/149‑1 (MQ=255)
                                                                                                                       gAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATcaccac                                                             <  6:74134/117‑1 (MQ=255)
                                                                                                                       gAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATcaccac                                                             >  5:74134/1‑117 (MQ=255)
                                                                                                                         aGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAAc                           <  2:380907/149‑1 (MQ=255)
                                                                                                                          gATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATg                                 >  6:238970/1‑142 (MQ=255)
                                                                                                                          gATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATg                                 <  5:238970/142‑1 (MQ=255)
                                                                                                                          gATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGaa                            >  7:167357/1‑147 (MQ=255)
                                                                                                                          gATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAAc                           <  5:150869/148‑1 (MQ=255)
                                                                                                                          gATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAAc                           >  6:150869/1‑148 (MQ=255)
                                                                                                                            tCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAg                         >  3:361729/1‑148 (MQ=255)
                                                                                                                              aaCAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGt                       <  1:223738/148‑1 (MQ=255)
                                                                                                                               aCAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTAt                     >  1:205385/1‑149 (MQ=255)
                                                                                                                                 aGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTc                   <  4:806/149‑1 (MQ=255)
                                                                                                                                   gCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCa                 <  1:40003/149‑1 (MQ=255)
                                                                                                                                       cTGAACATCATGCTGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGcc             >  3:115681/1‑149 (MQ=255)
                                                                                                                                         gAACATCATGCTGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCg            >  6:320841/1‑148 (MQ=255)
                                                                                                                                           aCAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGa          <  6:64153/148‑1 (MQ=255)
                                                                                                                                             aGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGa       <  3:327746/149‑1 (MQ=255)
                                                                                                                                             aGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGa       >  5:179801/1‑149 (MQ=255)
                                                                                                                                                 tgatgTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACt   >  5:324219/1‑149 (MQ=255)
                                                                                                                                                   atgTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACtt  <  1:6871/148‑1 (MQ=255)
                                                                                                                                                   atgTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACtt  <  5:199809/148‑1 (MQ=255)
                                                                                                                                                   atgTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACtt  <  1:328433/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CCGGCTTCGCCTTCATACACATAGAGGTCGAAGGTTTTCAGCACGATGCCGTTCAGGTAGTCGCCAACGTTCTGCACGAAGCCGTTGAGCAGGTGCAGGGTTTCGCCGCCCAGCAGCACGAAGATCAACAGGCCACTGAACAGCATGATGTTCAGGTTGGACAGGCGGCGGATGCCGTTCTCCACGCCCGATACCGCTGCCACGGTAGCCACGCCGGCCATGACCAGGATCACCACCAGCAGGTTGGTCTTGCTGTGGTCCATGCCGAACAGGTATTCCAGGCCGGACGATACTT  >  NC_002947/5772010‑5772304

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: ‑