Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,841,653 A→T 6.4% Y69* (TAT→TAA PP_5529 ← epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,841,6530AT 6.4% 613.8 / 22.4 186Y69* (TAT→TAAPP_5529epimerase
Reads supporting (aligned to +/- strand):  ref base A (106/68);  new base T (3/9);  total (109/77)
Fisher's exact test for biased strand distribution p-value = 2.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

Only 100 of 186 total aligned reads displayed.

CAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGCCTCCTGCGCGTAAC  >  NC_002947/2841505‑2841801
                                                                                                                                                    |                                                                                                                                                    
cAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCAt                                                                                                                                                     >  1:3180504/1‑150 (MQ=255)
  aTGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTATGCCAGCCCGGCTTAc                                                                                                                                                   <  1:5191481/150‑1 (MQ=255)
   tGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATAcc                                                                                                                                                  >  1:936538/1‑150 (MQ=255)
     aCATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAg                                                                                                                                                >  2:975732/1‑150 (MQ=255)
       aTGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGc                                                                                                                                              >  2:2173832/1‑150 (MQ=255)
         gCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTc                                                                                                                                            <  1:4790221/150‑1 (MQ=255)
           acaAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGa                                                                                                                                          <  1:1037684/150‑1 (MQ=255)
           acaAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGa                                                                                                                                          <  2:2701203/150‑1 (MQ=255)
           acaAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGa                                                                                                                                          <  1:4465157/150‑1 (MQ=255)
           acaAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGa                                                                                                                                          >  1:2302949/1‑150 (MQ=255)
             aaGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGc                                                                                                                                        >  2:1322661/1‑150 (MQ=255)
               gCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCgg                                                                                                                                      <  2:5288490/150‑1 (MQ=255)
                ccGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGc                                                                                                                                     >  1:4812968/1‑150 (MQ=255)
                 cGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCt                                                                                                                                    <  1:3872643/150‑1 (MQ=255)
                 cGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCt                                                                                                                                    >  1:2117417/1‑150 (MQ=255)
                 cGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCt                                                                                                                                    >  1:916768/1‑150 (MQ=255)
                     atatACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCa                                                                                                                                >  2:4802475/1‑150 (MQ=255)
                      tataCGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCat                                                                                                                               >  1:4390714/1‑150 (MQ=255)
                            cTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTc                                                                                                                         <  2:2448720/150‑1 (MQ=255)
                            cTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTc                                                                                                                         >  1:1705217/1‑150 (MQ=255)
                            cTCATAACGACGTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTc                                                                                                                         <  2:5393955/150‑1 (MQ=255)
                             tCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATTTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCAGCACCGGAGAGACGATGCTTTAGCGCCGACTGGTAAGCCGGGATGGCATACAAGGCTCGAGCGGCTGTCATATCTTCg                                                                                                                        >  1:1060299/1‑150 (MQ=255)
                               aTAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGaa                                                                                                                      >  1:1549511/1‑150 (MQ=255)
                               aTAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGaa                                                                                                                      >  2:2650096/1‑150 (MQ=255)
                               aTAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGaa                                                                                                                      >  1:1549244/1‑150 (MQ=255)
                                       ttGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGaa                                                                                                              >  2:2607478/1‑150 (MQ=255)
                                          tgCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGa                                                                                                           <  1:2607478/150‑1 (MQ=255)
                                          tgCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGa                                                                                                           >  1:996288/1‑150 (MQ=255)
                                           gCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGAc                                                                                                          >  1:1009548/1‑150 (MQ=255)
                                            cACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGAcc                                                                                                         >  1:3195366/1‑150 (MQ=255)
                                            cACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGAcc                                                                                                         >  2:4450338/1‑150 (MQ=255)
                                             aCGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCg                                                                                                        >  2:3799289/1‑150 (MQ=255)
                                                      aGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCtcga                                                                                               <  2:712303/150‑1 (MQ=255)
                                                       gTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCTGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCtcgat                                                                                              <  1:21001/150‑1 (MQ=255)
                                                       gTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCtcgat                                                                                              >  1:859988/1‑150 (MQ=255)
                                                       gTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCtcgat                                                                                              >  2:2501873/1‑150 (MQ=255)
                                                       gTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCtcgat                                                                                              >  2:5097952/1‑150 (MQ=255)
                                                           ggTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTcgatcgat                                                                                          <  1:4585821/150‑1 (MQ=255)
                                                                tgtgTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGAcc                                                                                     >  2:2120078/1‑150 (MQ=255)
                                                                tgtgTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGAcc                                                                                     >  1:71964/1‑150 (MQ=255)
                                                                    tCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTcc                                                                                 >  1:4416666/1‑150 (MQ=255)
                                                                       gTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAGACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAac                                                                              >  1:2444009/1‑150 (MQ=255)
                                                                       gTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAac                                                                              >  1:5271994/1‑150 (MQ=255)
                                                                       gTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAac                                                                              >  1:3026465/1‑150 (MQ=255)
                                                                         gCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAacac                                                                            <  1:3700577/150‑1 (MQ=255)
                                                                         gCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAacac                                                                            <  1:3480092/150‑1 (MQ=255)
                                                                           cTTCGATTAGGAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCt                                                                          >  2:2069651/1‑150 (MQ=255)
                                                                            ttCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCTATCGATGGACCTTCCAACACCTc                                                                         >  2:4009044/1‑150 (MQ=255)
                                                                            ttCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTc                                                                         <  1:1143597/150‑1 (MQ=255)
                                                                            ttCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTc                                                                         >  1:3835553/1‑150 (MQ=255)
                                                                               gATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAg                                                                      <  1:3301785/150‑1 (MQ=255)
                                                                               gATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAg                                                                      >  2:3076155/1‑150 (MQ=255)
                                                                               gATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAg                                                                      >  1:2212958/1‑150 (MQ=255)
                                                                               gATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAg                                                                      >  2:1847485/1‑150 (MQ=255)
                                                                               caatAGAAAAACTCCGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCCATCGATGGACCTTCCAACACCTCCAg                                                                      <  2:5065040/147‑1 (MQ=255)
                                                                                 ttAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGtt                                                                    <  2:3292595/150‑1 (MQ=255)
                                                                                   aGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGa                                                                  >  1:423473/1‑150 (MQ=255)
                                                                                    gAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGAc                                                                 >  2:1053758/1‑150 (MQ=255)
                                                                                    gAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGAc                                                                 <  2:5154731/150‑1 (MQ=255)
                                                                                    gAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGAc                                                                 >  1:3619641/1‑150 (MQ=255)
                                                                                     aaaaaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGTCTGGCATAACAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATTACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGaac                                                                >  1:5363330/1‑148 (MQ=255)
                                                                                     aaaaaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGAcc                                                                >  1:3185967/1‑150 (MQ=255)
                                                                                     aaaaaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGAcc                                                                >  1:3186005/1‑150 (MQ=255)
                                                                                       aaaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGt                                                              >  2:2074901/1‑150 (MQ=255)
                                                                                       aaaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGt                                                              >  1:3335685/1‑150 (MQ=255)
                                                                                        aaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTa                                                             <  1:3058760/150‑1 (MQ=255)
                                                                                        aaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTa                                                             <  1:620605/150‑1 (MQ=255)
                                                                                          cTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGa                                                           <  2:1548451/150‑1 (MQ=255)
                                                                                           tCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGaa                                                          >  2:1224062/1‑150 (MQ=255)
                                                                                           tCGGAACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGCAGGACCTTCCAACACCTCCAGTTGACCGTAGaa                                                          >  1:5221631/1‑150 (MQ=255)
                                                                                            cGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGaaa                                                         <  2:418315/150‑1 (MQ=255)
                                                                                              gTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGc                                                       <  1:4231221/150‑1 (MQ=255)
                                                                                                 cTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAg                                                    >  1:2061789/1‑150 (MQ=255)
                                                                                                  tCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCa                                                                                                                                                 <  1:4095225/56‑1 (MQ=255)
                                                                                                  tCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGa                                                   >  1:1288371/1‑150 (MQ=255)
                                                                                                  tCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGa                                                   >  2:375737/1‑150 (MQ=255)
                                                                                                    cGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGGTGACCGTAGAAAGCCAGACg                                                 >  2:2931588/1‑150 (MQ=255)
                                                                                                     gCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACgt                                                >  1:3681419/1‑150 (MQ=255)
                                                                                                     gCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACgt                                                <  2:4416330/150‑1 (MQ=255)
                                                                                                     gCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCAGATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACgt                                                >  1:1588229/1‑150 (MQ=255)
                                                                                                       aCCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACgtgt                                              <  2:2909785/150‑1 (MQ=255)
                                                                                                               gaCGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGAcgtcgt                                      <  2:3579128/150‑1 (MQ=255)
                                                                                                               gaCGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGAcgtcgt                                      <  1:1155770/150‑1 (MQ=255)
                                                                                                                aCGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGACACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTg                                     >  2:2393391/1‑150 (MQ=255)
                                                                                                                 cGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGc                                    <  1:3799289/150‑1 (MQ=255)
                                                                                                                       tttAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTtct                              >  2:1989219/1‑150 (MQ=255)
                                                                                                                       tttAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGAGGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTtct                              >  2:1526012/1‑150 (MQ=255)
                                                                                                                           gcgcCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTg                          >  2:1128362/1‑150 (MQ=255)
                                                                                                                           gcgcCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTg                          >  2:3868322/1‑150 (MQ=255)
                                                                                                                           gcgcCGCCTGGTAAGACGGGCTGGCAAACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTACCAACACCACCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCgag                          >  2:3311986/1‑148 (MQ=255)
                                                                                                                                  cTGGTATGCCAGCCCGGCTTACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAgg                   >  2:1663423/1‑150 (MQ=255)
                                                                                                                                   tGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGgc                  <  2:916768/150‑1 (MQ=255)
                                                                                                                                     tgatcacGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGgcgc                >  2:5163902/7‑150 (MQ=255)
                                                                                                                                     gTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGgcgc                >  1:4235338/1‑150 (MQ=255)
                                                                                                                                          ccGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGcctcct           >  2:3040743/1‑150 (MQ=255)
                                                                                                                                          ccGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGcctcct           >  1:1410721/1‑150 (MQ=255)
                                                                                                                                              gCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGCCTCCTgcgc       <  1:3524202/150‑1 (MQ=255)
                                                                                                                                              gCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGCCTCCTgcgc       <  2:4166833/150‑1 (MQ=255)
                                                                                                                                                tGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGCCTCCTGCGCGt     >  1:4991116/1‑150 (MQ=255)
                                                                                                                                                   cATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGCCTCCTGCGCGTata  <  2:4390714/150‑3 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGCCTCCTGCGCGTAAC  >  NC_002947/2841505‑2841801

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