Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,841,641 A→T 7.1% A73A (GCT→GCA PP_5529 ← epimerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,841,6410AT 7.1% 582.8 / 25.0 183A73A (GCT→GCAPP_5529epimerase
Reads supporting (aligned to +/- strand):  ref base A (105/65);  new base T (4/9);  total (109/74)
Fisher's exact test for biased strand distribution p-value = 3.91e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 184 total aligned reads displayed.

GATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGC  >  NC_002947/2841492‑2841787
                                                                                                                                                     |                                                                                                                                                  
     ttCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCgg                                                                                                                                               >  2:1352833/1‑150 (MQ=255)
      tCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCAGCCTGGTATGccagc                                                                                                                                              >  1:505068/1‑147 (MQ=255)
        cAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTCAGCCGGGCt                                                                                                                                            >  1:4335574/1‑150 (MQ=255)
          gTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTgg                                                                                                                                          >  2:2624576/1‑150 (MQ=255)
             cAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCAt                                                                                                                                       >  1:3180504/1‑150 (MQ=255)
               aTGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTATGCCAGCCCGGCTTAc                                                                                                                                     <  1:5191481/150‑1 (MQ=255)
                  aCATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAg                                                                                                                                  >  2:975732/1‑150 (MQ=255)
                    aTGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTATGCCAGCCCGGCTTACCAGGc                                                                                                                                <  2:1906481/150‑1 (MQ=255)
                    aTGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGc                                                                                                                                >  2:2173832/1‑150 (MQ=255)
                     tGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTATGCCAGCCCGGCTTGCCAGGcg                                                                                                                               <  2:2024492/150‑2 (MQ=255)
                     tGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTATGCCAGCCCGGCTTACCAGGcg                                                                                                                               <  2:373363/150‑2 (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:2404946/1‑150 (MQ=255)
                         caAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAg                                                                                                                           <  1:3533268/150‑1 (MQ=255)
                            gCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCgg                                                                                                                        <  2:5288490/150‑1 (MQ=255)
                              cGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCt                                                                                                                      >  1:2484112/1‑150 (MQ=255)
                              cGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCt                                                                                                                      <  1:3872643/150‑1 (MQ=255)
                              cGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCt                                                                                                                      >  1:2117417/1‑150 (MQ=255)
                                  atatACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCa                                                                                                                  >  2:4802475/1‑150 (MQ=255)
                                   tataCGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCat                                                                                                                 >  1:4390714/1‑150 (MQ=255)
                                     taCGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCatat                                                                                                               >  1:4194663/1‑150 (MQ=255)
                                        gCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCtt                                                                                                            <  2:5158759/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:1549244/1‑150 (MQ=255)
                                            aTAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGaa                                                                                                        >  1:1549511/1‑150 (MQ=255)
                                            aTAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGaa                                                                                                        >  2:2650096/1‑150 (MQ=255)
                                                 gACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTa                                                                                                   >  2:1186033/1‑150 (MQ=255)
                                                    ttGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGaa                                                                                                <  1:3853129/150‑1 (MQ=255)
                                                        gCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGAc                                                                                            >  2:2553611/1‑150 (MQ=255)
                                                        gCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGAc                                                                                            >  1:1009548/1‑150 (MQ=255)
                                                         cACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGAcc                                                                                           >  2:4450338/1‑150 (MQ=255)
                                                         cACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGAcc                                                                                           >  1:3195287/1‑150 (MQ=255)
                                                          aCGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCg                                                                                          >  2:3799289/1‑150 (MQ=255)
                                                                 tCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCtc                                                                                   >  1:712303/1‑150 (MQ=255)
                                                                   aGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCtcga                                                                                 <  2:712303/150‑1 (MQ=255)
                                                                    gTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCTGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCtcgat                                                                                <  1:21001/150‑1 (MQ=255)
                                                                    gTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCtcgat                                                                                >  2:5097952/1‑150 (MQ=255)
                                                                    gTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCtcgat                                                                                >  1:859988/1‑150 (MQ=255)
                                                                        ggTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTcgatcgat                                                                            <  1:4585821/150‑1 (MQ=255)
                                                                        ggTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTcgatcgat                                                                            <  1:3269481/150‑1 (MQ=255)
                                                                        ggTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTcgatcgat                                                                            >  2:2467602/1‑150 (MQ=255)
                                                                           gATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGa                                                                         <  1:3768105/150‑1 (MQ=255)
                                                                             tgtgTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGAcc                                                                       >  2:2120078/1‑150 (MQ=255)
                                                                                 tCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTcc                                                                   >  1:4416666/1‑150 (MQ=255)
                                                                                   ggTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCaa                                                                 <  1:2173832/150‑1 (MQ=255)
                                                                                    gTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAac                                                                >  1:891768/1‑150 (MQ=255)
                                                                                     tGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAaca                                                               >  1:5319870/1‑150 (MQ=255)
                                                                                      gcattCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAacac                                                              <  1:929020/147‑1 (MQ=255)
                                                                                      gCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAacac                                                              <  1:3480092/150‑1 (MQ=255)
                                                                                         ttCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCTATCGATGGACCTTCCAACACCTc                                                           >  2:4009044/1‑150 (MQ=255)
                                                                                         ttCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTc                                                           >  1:3835498/1‑150 (MQ=255)
                                                                                         ttCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTc                                                           >  2:3360042/1‑150 (MQ=255)
                                                                                         ttCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTc                                                           <  1:1143597/150‑1 (MQ=255)
                                                                                            gATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAg                                                        >  1:2212958/1‑150 (MQ=255)
                                                                                            gATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAg                                                        >  2:3076155/1‑150 (MQ=255)
                                                                                            gATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAg                                                        >  2:1847485/1‑150 (MQ=255)
                                                                                            gATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAg                                                        <  1:3301919/150‑1 (MQ=255)
                                                                                            caatAGAAAAACTCCGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCCATCGATGGACCTTCCAACACCTCCAg                                                        <  2:5065040/147‑1 (MQ=255)
                                                                                             aTTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGt                                                       >  1:560662/1‑150 (MQ=255)
                                                                                              ttAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGtt                                                      <  2:3292595/150‑1 (MQ=255)
                                                                                                aGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGa                                                    >  1:423473/1‑150 (MQ=255)
                                                                                                 gAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGAc                                                   <  2:5154731/150‑1 (MQ=255)
                                                                                                  aaaaaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGTCCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGAcc                                                  >  1:4909290/1‑150 (MQ=255)
                                                                                                  aaaaaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGAcc                                                  >  1:3185967/1‑150 (MQ=255)
                                                                                                    aaaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGt                                                >  2:2074901/1‑150 (MQ=255)
                                                                                                     aaCTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTa                                               <  1:3058760/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:1101367/150‑1 (MQ=255)
                                                                                                         cGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGaaa                                           >  2:3017944/1‑150 (MQ=255)
                                                                                                         cGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGaaa                                           <  2:418315/150‑1 (MQ=255)
                                                                                                              cTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAg                                      >  1:302983/1‑150 (MQ=255)
                                                                                                              cTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAg                                      >  1:2061789/1‑150 (MQ=255)
                                                                                                               tCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCa                                                                                                                                   <  1:4095225/56‑1 (MQ=255)
                                                                                                               tCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGa                                     >  2:375737/1‑150 (MQ=255)
                                                                                                               tCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGa                                     >  1:1288371/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:2910048/150‑1 (MQ=255)
                                                                                                                    aCCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACgtgt                                <  1:3360042/150‑1 (MQ=255)
                                                                                                                    aCCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACgtgt                                <  2:2909785/150‑1 (MQ=255)
                                                                                                                            gaCGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGAcgtcgt                        <  1:1155770/150‑1 (MQ=255)
                                                                                                                             aCGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGACACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTg                       >  2:2393391/1‑150 (MQ=255)
                                                                                                                              cGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGc                      <  1:953392/150‑1 (MQ=255)
                                                                                                                              cGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGc                      <  2:2404946/150‑1 (MQ=255)
                                                                                                                                    tttAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGAGGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTtct                >  2:1526012/1‑150 (MQ=255)
                                                                                                                                     ttAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTtctc               >  2:4401009/1‑150 (MQ=255)
                                                                                                                                        gcgcCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTg            >  2:3868322/1‑150 (MQ=255)
                                                                                                                                               cTGGTATGCCAGCCCGGCTTACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAgg     >  2:1663423/1‑150 (MQ=255)
                                                                                                                                                tGGTATGCCAGCCCGGCTTACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTtctc               >  2:1045536/1‑139 (MQ=255)
                                                                                                                                                tGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGgc    >  1:2646757/1‑150 (MQ=255)
                                                                                                                                                tGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGgc    <  2:916768/150‑1 (MQ=255)
                                                                                                                                                  gTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGgcgc  >  1:4235338/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GATTGTTCCAGTGCAATGACATGCACAAGCCGAGATATACGCTCATAACGACTTGTGCACGCGTGTCAGTTAGGTGATGTGTCGGTGCCTTCGATTAGAAAAACTCGGTACTCCGCACCGGAGAGACGATGCTTTAGCGCCGCCTGGTAAGCCGGGCTGGCATACCAGGCTCGAGCGGCTGTCATATCTTCGAAACGTAGAATGACCGCTCCCTCGATCGATGGACCTTCCAACACCTCCAGTTGACCGTAGAAAGCCAGACGTGTGACGTCGTGCCCTTCTCGTGCGGCAGGCGC  >  NC_002947/2841492‑2841787

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