Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,252,194 C→A 6.0% intergenic (‑162/+78) PP_3726 ← / ← rocE enoyl‑CoA hydratase/isomerase family protein/amino acid permease RocE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,252,1940CA 6.0% 210.2 / 2.5 67intergenic (‑162/+78)PP_3726/rocEenoyl‑CoA hydratase/isomerase family protein/amino acid permease RocE
Reads supporting (aligned to +/- strand):  ref base C (37/26);  new base A (2/2);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

TTGTGATGGTGTTGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCA  >  NC_002947/4252050‑4252338
                                                                                                                                                   |                                                                                                                                                
ttgtGATGGTGTTGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCTCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCAc                                                                                                                                              >  8:355213/1‑149 (MQ=255)
        gagtTGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGAGTCTGGCGTTAAACAA‑‑GGCCGTTCGTCATCACGGCTGTATTGCCTGTGAGATCGCGAACCGCGTAgg                                                                                                                                      >  8:93437/3‑149 (MQ=255)
          gtTGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACAGAA‑‑GGGCGTTCGTCATGACGGCGTTATGGCCTGTGAGAGGGCGCACGGCGTAgagg                                                                                                                                    >  4:333944/1‑146 (MQ=255)
          gtTGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACcgcg                                                                                                                                          <  3:333944/143‑1 (MQ=255)
            tGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAggcggc                                                                                                                                  >  1:113764/1‑149 (MQ=255)
              gtggaCTTTATGAGGGAGATACTAACTGAGTGTTCAGCAAGGCGTGATGGGCTGTTACTAGTTCGGACAACTCGGTGGCTGGCGTTACACAA‑‑GGGCGTTCGTCATTACGGCTTTATTGCCTGTGAGATCTCGCACCGCGTAGGCGgcgc                                                                                                                                <  8:85632/147‑1 (MQ=255)
              gaggaCTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCTCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGgcgc                                                                                                                                <  7:355213/149‑1 (MQ=255)
              gaggaCTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGgcgc                                                                                                                                >  3:24991/1‑149 (MQ=255)
                    ttGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGAcc                                                                                                                          >  6:327416/1‑149 (MQ=255)
                    ttGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGAcc                                                                                                                          <  5:233777/149‑1 (MQ=255)
                      gatgaGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCTCTCGACCt                                                                                                                         >  5:279777/1‑148 (MQ=255)
                        tgaGCGACATACTAAACGCGAGTTCACCAACGCGGGATTGGCGGTTACAACATCGCCCCACTCGGGGTCTGGCGTTACACAACTAGGCGTTCGGCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTc                                                                                                                        <  1:299125/149‑1 (MQ=255)
                          aGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTcaca                                                                                                                     <  3:95286/148‑1 (MQ=255)
                              agaTACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGgcgc                                                                                                                <  5:147199/149‑1 (MQ=255)
                                   cTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGaaa                                                                                                           <  5:164447/149‑1 (MQ=255)
                                              ttCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGAc                                                                                                 <  2:113764/148‑1 (MQ=255)
                                                  gCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATTGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAaca                                                                                            >  4:65901/1‑149 (MQ=255)
                                                    aaGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACAcc                                                                                          >  6:161303/1‑149 (MQ=255)
                                                     aGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGCCCTCCCAGGCGCTGAAAAAGCAGCGACCAACACCt                                                                                         >  4:19566/1‑149 (MQ=255)
                                                      gtcGGGATGTGCTGTTACAACTTCGCACATCTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCTTGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTg                                                                                        <  7:93437/147‑1 (MQ=255)
                                                      ggCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTg                                                                                        <  1:340963/149‑1 (MQ=255)
                                                      ggCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTg                                                                                        >  7:346857/1‑149 (MQ=255)
                                                         gggATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCcg                                                                                     >  3:121404/1‑149 (MQ=255)
                                                          ggATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGTCGGCTTTATGGCCTGTGAGATCGCGCACCGCTTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCcga                                                                                    <  4:24991/149‑1 (MQ=255)
                                                              gggCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCcgacgac                                                                                >  6:245355/1‑149 (MQ=255)
                                                                    ttACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTg                                                                          >  1:60612/1‑149 (MQ=255)
                                                                      aCAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGTGACCAACACCTGCCGACGACGAGCTgcg                                                                        >  7:64074/1‑149 (MQ=255)
                                                                        aaCTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTgcggc                                                                      >  6:342312/1‑149 (MQ=255)
                                                                        aaCTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTgcggc                                                                      >  6:336511/1‑149 (MQ=255)
                                                                        aaCTTAGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTAAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTgcggc                                                                      >  5:261256/1‑149 (MQ=255)
                                                                         aCTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGc                                                                                       >  5:51298/1‑131 (MQ=255)
                                                                         aCTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGc                                                                                       <  6:51298/131‑1 (MQ=255)
                                                                                     tCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTgg                                                         >  4:47128/1‑149 (MQ=255)
                                                                                       ggTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGgtg                                                       >  2:258561/1‑149 (MQ=255)
                                                                                          gtCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGc                                                    >  7:254366/1‑149 (MQ=255)
                                                                                             tGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGATTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTgg                                                 <  5:245355/149‑1 (MQ=255)
                                                                                               gCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTggg                                                >  1:66907/1‑149 (MQ=255)
                                                                                               gCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGt                                               <  2:60612/149‑1 (MQ=255)
                                                                                               gCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGt                                               <  5:336511/149‑1 (MQ=255)
                                                                                               gCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGt                                               <  5:342312/149‑1 (MQ=255)
                                                                                                      acaA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAgg                                        >  6:87795/1‑149 (MQ=255)
                                                                                                      acaA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAgg                                        >  5:16704/1‑149 (MQ=255)
                                                                                                        aa‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCTGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTggg                                                <  6:64337/139‑1 (MQ=255)
                                                                                                        aa‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTggg                                                >  5:64337/1‑139 (MQ=255)
                                                                                                         a‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGccc                                     <  5:161303/149‑1 (MQ=255)
                                                                                                                gTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGc                                <  1:258561/149‑1 (MQ=255)
                                                                                                                    gTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGa                             <  3:19566/148‑1 (MQ=255)
                                                                                                                       agtaCGTCTTTATGGCCTTTGAGATCGCGCACTGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGAcggc                         <  6:261256/146‑1 (MQ=255)
                                                                                                                        tGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGAcggcg                        >  3:303062/1‑149 (MQ=255)
                                                                                                                         gatggCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGAcggcgg                       >  6:164330/4‑149 (MQ=255)
                                                                                                                           cGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGTCCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGtt                     <  3:47128/149‑1 (MQ=255)
                                                                                                                            ggCTTTATGGCCTGTGAGATCGCGCACCGAGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGAGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTCCCGACGGCGGtt                     >  2:262351/1‑148 (MQ=255)
                                                                                                                             gCTTTTTGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTGGACCTCACAGGGGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACACCTGGTGTGCTGGGTGCGAATGCCCCCTGCCGACGGCGGTTGa                   <  6:279777/149‑1 (MQ=255)
                                                                                                                                ttATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCACCACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGCCGGCGGTTGACg                 >  7:291131/1‑148 (MQ=255)
                                                                                                                                ttATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACg                 >  6:73143/1‑148 (MQ=255)
                                                                                                                                   agcgCCTGTGAGGTCGAGCGCCGCCTACGCGGTGCGCGATCCCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACgaag              <  5:321479/146‑1 (MQ=37)
                                                                                                                                      gCCTGTGAGATTGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACgaagaa            <  3:65901/148‑1 (MQ=255)
                                                                                                                                        cTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACgaagaa            <  4:322741/146‑1 (MQ=255)
                                                                                                                                        cTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGaagaagag         >  3:322741/1‑147 (MQ=255)
                                                                                                                                         tgtgAGGTCGAGCGCCGCCTACGCGGTGCGCGATCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTAc        >  1:308822/1‑149 (MQ=38)
                                                                                                                                          gtgAGGTCGAGCGCCGCCTACGCGGTGCGCGATCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTAcg       >  2:121592/1‑149 (MQ=38)
                                                                                                                                          gtgAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTAcg       >  3:189800/1‑149 (MQ=255)
                                                                                                                                          gtgAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTAcg       >  3:177973/1‑149 (MQ=255)
                                                                                                                                             aggtCGAGCGCCGCCTACGCGGTGCGCGATCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTAcgcc     <  2:308822/145‑1 (MQ=38)
                                                                                                                                              gaTCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTAcgccg    >  2:291308/1‑148 (MQ=255)
                                                                                                                                               aTCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGCCCAACACCTGCCGACGACGAGTTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCa  >  5:48008/1‑149 (MQ=255)
                                                                                                                                               aTCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCa  >  1:211155/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TTGTGATGGTGTTGGAGGACTTGATGAGGGAGATACTAAACGAGTGTTCAGCAAGGCGGGATGGGCTGTTACAACTTCGCACAACTCGGTGTCTGGCGTTACACAA‑‑GGGCGTTCGTCATGACGGCTTTATGGCCTGTGAGATCGCGCACCGCGTAGGCGGCGCTCGACCTCACAGGCGCTGAAAATGCAGCGACCAACACCTGCCGACGACGAGCTGCGGCGTCAAACCCCTGGTGTGCTGGGTGCGAAGGCCCCCTGCCGACGGCGGTTGACGAAGAAGTACGCCGCA  >  NC_002947/4252050‑4252338

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