Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,257,874 A→G 6.1% intergenic (‑73/+71) PP_4634 ← / ← PP_4635 alpha/beta family hydrolase/trans‑2‑enoyl‑CoA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,257,8740AG 6.1% 234.3 / 4.4 83intergenic (‑73/+71)PP_4634/PP_4635alpha/beta family hydrolase/trans‑2‑enoyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base A (40/37);  new base G (2/3);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

AGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCG  >  NC_002947/5257726‑5258020
                                                                                                                                                    |                                                                                                                                                        
agcagAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGaa                                                                                                                                                          >  5:173182/1‑149 (MQ=255)
    gAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCg                                                                                                                                                      >  2:343389/1‑149 (MQ=255)
      aGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCt                                                                                                                                                    <  6:272343/149‑1 (MQ=255)
       gccgccGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGATAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTc                                                                                                                                                   >  8:206062/1‑149 (MQ=255)
         cgccgAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGCAAGCCCGCTc‑‑cc                                                                                                                                               <  8:115045/149‑2 (MQ=255)
          gccgAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑t                                                                                                                                             >  2:191609/1‑148 (MQ=255)
           ccgAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TAc                                                                                                                                           <  2:238525/149‑1 (MQ=255)
            cgAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCTc‑‑cca                                                                                                                                              >  3:348415/1‑145 (MQ=255)
              aGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAgg                                                                                                                                        <  2:314870/149‑1 (MQ=255)
                gCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAgggg                                                                                                                                      <  5:22485/149‑1 (MQ=255)
                gCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAgggg                                                                                                                                      <  3:200699/149‑1 (MQ=255)
                 cTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGATAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAggggg                                                                                                                                     <  7:206062/149‑1 (MQ=255)
                  tGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑Ta                                                                                                                                            >  1:239693/1‑141 (MQ=255)
                  tGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTGCTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑Ta                                                                                                                                            <  2:239693/141‑1 (MQ=255)
                    tATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGgag                                                                                                                                   >  4:148543/1‑148 (MQ=255)
                    tATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGga                                                                                                                                    <  7:259373/147‑1 (MQ=255)
                      tCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAAc                                                                                                                                <  4:406477/149‑1 (MQ=255)
                       ccGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGACTCTTCGCGGGTGAACCCGCTC‑‑C‑‑AACAGGGGGAGAAcc                                                                                                                               >  8:51487/1‑149 (MQ=255)
                        cGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑t                                                                                                                            >  5:83150/1‑149 (MQ=255)
                        cGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCAACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑t                                                                                                                            >  1:306526/1‑149 (MQ=255)
                             gCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTgg                                                                                                                        >  2:153618/1‑148 (MQ=255)
                                gcgGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGc                                                                                                                    <  8:399010/149‑1 (MQ=255)
                                  ggggTGGCGAAATGGGCTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGGGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCgg                                                                                                                  <  1:287036/149‑1 (MQ=255)
                                   gggTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCgg                                                                                                                  <  1:374940/148‑1 (MQ=255)
                                     gTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCt                                                                                                               <  6:306163/149‑1 (MQ=255)
                                        gCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCCAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTga                                                                                                            >  1:126410/1‑148 (MQ=255)
                                        gCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTg                                                                                                             <  4:154131/148‑1 (MQ=255)
                                        gCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGc                                                                                                            >  6:229003/1‑149 (MQ=255)
                                           aaaTGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCcg                                                                                                         >  2:224413/1‑149 (MQ=255)
                                            aaTGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCtt                                                                                                              >  5:26867/1‑143 (MQ=255)
                                            aaTGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCtt                                                                                                              <  6:26867/143‑1 (MQ=255)
                                                  gTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAGCTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAagag                                                                                                  >  6:326808/1‑149 (MQ=255)
                                                   tAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCTC‑‑C‑‑CACAGGTTCTCCCCC‑‑TGTAGGAGCGGGTTCACCCGCGAagag                                                                                                  <  3:303817/148‑1 (MQ=14)
                                                   tAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCTC‑‑C‑‑CACAGGTTCTCCCCC‑‑TGTAGGAGCGGGTTCACCCGCGAagag                                                                                                  >  4:303817/1‑148 (MQ=14)
                                                     gCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAg                                                                                                                     >  6:273805/1‑127 (MQ=255)
                                                     gCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAg                                                                                                                     <  5:273805/127‑1 (MQ=255)
                                                       gggTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACCGGCCTCTTCGCGGGCAAGCCCGCTC‑‑C‑‑CACAGGTTCTCCCCC‑‑TGTAGGAGCGGGTTCACCCGCGAAGAGGccag                                                                                             <  4:348415/149‑3 (MQ=14)
                                                              aTGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCTCCAAAACTAGCAAAGATTTTGTAGCGTGTACTGGCCTGTGCGCGGGTGAACCCGCTC‑‑C‑‑TACTGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGc                                                                                      <  4:407597/149‑1 (MQ=39)
                                                                      aCTCAAGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAGGAGGCCGGTACAGGCTTACCCGa                                                                              <  3:95461/149‑1 (MQ=255)
                                                                       cTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTATCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTCCTC‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTAccc                                                                                <  1:343239/148‑1 (MQ=255)
                                                                        tCAAGGCCCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGCCCTCATCGCGGGTGAACCCCCTC‑‑C‑‑TACACGGGGAGAACC‑‑TGTGGGAGCGCGCTTGCCCGCGAAGAGGCAGGTACAGGCTTACCCGAt                                                                             >  6:322057/1‑148 (MQ=255)
                                                                          aTGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTa                                                                           >  1:203313/1‑148 (MQ=255)
                                                                              cgcATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCt                                                                      >  5:375399/1‑149 (MQ=255)
                                                                              cgcATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCt                                                                      <  5:229003/149‑1 (MQ=255)
                                                                               gcATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTc                                                                     >  4:87436/1‑149 (MQ=255)
                                                                                   gTCCTGTTGTTGTGCGCCAAAAGTCGCAAAGATTTTGTAGCCTTTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGAc                                                                  <  7:51487/148‑1 (MQ=255)
                                                                                       tgttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGt                                                             >  8:183983/1‑149 (MQ=255)
                                                                                       tgttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑A‑‑TACAGGGGGAGAACC‑‑TGTGTGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACACGATTACAGCTCGACGCAg                                                              >  7:279139/1‑148 (MQ=255)
                                                                                       tgttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTAATGGCCTCTTACCGGGTGCAACCCCGC‑‑C‑‑TACAGGGGGAGAAAA‑‑TGTGGGATCGGCCATGCCACACAAGAGCCCGGAACATGCATACCCGATTACATCTATACGCAGt                                                             >  8:215018/1‑149 (MQ=255)
                                                                                         ttgttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTc                                                            <  6:375399/148‑1 (MQ=255)
                                                                                            ttgtGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGaa                                                         <  3:253494/148‑1 (MQ=255)
                                                                                              gtgCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACg                                                                 <  5:211068/138‑1 (MQ=255)
                                                                                              gtgCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACg                                                                 >  6:211068/1‑138 (MQ=255)
                                                                                               tgCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACtt                                                      >  2:144340/1‑148 (MQ=255)
                                                                                                gcgcCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑CGGTACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACtt                                                      >  1:125529/1‑149 (MQ=255)
                                                                                                 cgcCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTca                                                    <  7:183983/148‑1 (MQ=255)
                                                                                                  gcCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTcac                                                   >  6:99141/1‑148 (MQ=255)
                                                                                                   ccAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCCGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACAt                                                 >  4:363119/1‑149 (MQ=255)
                                                                                                     aaaaCTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑AACAGGGGGAGCACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGGACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCg                                               >  7:141517/1‑149 (MQ=255)
                                                                                                        ccttgCAATGATTTTGTAGCCTTTACTGGCCGATTCGCGGGTGAACCCGCTC‑‑C‑‑TACCGGGGGATAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCTATTACATCTCGACGCAGTCGAACTTCACATCGGTg                                            <  7:215018/145‑1 (MQ=38)
                                                                                                          tAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGc                                          >  7:397636/1‑149 (MQ=255)
                                                                                                             cAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCAc                                       <  1:224413/149‑1 (MQ=255)
                                                                                                              aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACCTGTGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCa                                        >  1:126038/1‑149 (MQ=255)
                                                                                                              aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACg                                      >  7:159754/1‑149 (MQ=255)
                                                                                                              aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACg                                      >  1:415501/1‑149 (MQ=255)
                                                                                                              aaaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACCCAGTCGAACATCACATCGGTGGCCAcc                                      >  4:224884/1‑148 (MQ=255)
                                                                                                               aaGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGCACTTCACATCGGTGGCCACGt                                     >  1:205446/1‑149 (MQ=255)
                                                                                                                 gATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCt                                   <  5:312260/149‑1 (MQ=255)
                                                                                                                 gATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCt                                   <  3:148543/149‑1 (MQ=255)
                                                                                                                  aTTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCtt                                  >  3:337220/1‑149 (MQ=255)
                                                                                                                   ttttGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCttt                                 <  1:191609/149‑1 (MQ=255)
                                                                                                                   ttttGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCttt                                 >  4:138215/1‑149 (MQ=255)
                                                                                                                       gTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTc                              <  6:83150/148‑1 (MQ=255)
                                                                                                                          gCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTa                           <  3:138215/148‑1 (MQ=255)
                                                                                                                          gCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAg                          >  4:270990/1‑149 (MQ=255)
                                                                                                                          gCCTCTACTGTCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCAGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTa                           <  8:279139/148‑1 (MQ=255)
                                                                                                                               tACTGGCCTCTTCGGGGGTGATCCCGGTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCTGGAAGAGGCCGGTAGAGGCTTACCCGATTAGAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCaa                      <  5:322057/148‑1 (MQ=255)
                                                                                                                               tACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAAc                     <  1:153618/149‑1 (MQ=255)
                                                                                                                                 cTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACg                    <  3:87436/148‑1 (MQ=255)
                                                                                                                                      ctctTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcg              >  1:227049/1‑149 (MQ=255)
                                                                                                                                       tctTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcgc             >  6:311440/1‑149 (MQ=255)
                                                                                                                                         ttCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcgcgc           >  3:186527/1‑149 (MQ=255)
                                                                                                                                           cgcgGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAgcgcgc           <  7:32180/147‑1 (MQ=255)
                                                                                                                                                 tGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGcc   <  2:227049/149‑1 (MQ=255)
                                                                                                                                                   aaCCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCg  >  6:27991/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                        
AGCAGAAGCCGCCGAGGCTGTATCCGTCGGCCGCGGGGTGGCGAAATGGGGTAGCGGGTTGGATGAGGCTACTCATGGCGCATGTCCTGTTGTTGTGCGCCAAAACTAGCAAAGATTTTGTAGCCTGTACTGGCCTCTTCGCGGGTGAACCCGCTC‑‑C‑‑TACAGGGGGAGAACC‑‑TGTGGGAGCGGGCTTGCCCGCGAAGAGGCCGGTACAGGCTTACCCGATTACAGCTCGACGCAGTCGAACTTCACATCGGTGGCCACGTCTTTGTCGTAGTCAACGTCAGCGCGCTCGAAGCCG  >  NC_002947/5257726‑5258020

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