Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,236,545 A→T 7.3% A76A (GCT→GCA)  catA‑I ← catechol 1,2‑dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,236,5450AT 7.3% 185.9 / 3.3 64A76A (GCT→GCA) catA‑Icatechol 1,2‑dioxygenase
Reads supporting (aligned to +/- strand):  ref base A (23/36);  new base T (3/2);  total (26/38)
Fisher's exact test for biased strand distribution p-value = 3.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GTCGTCCATGCGCGCTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGCGCGG  >  NC_002947/4236398‑4236691
                                                                                                                                                   |                                                                                                                                                  
gtcgtcCATGCGCGTTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGgccag                                                                                                                                                   <  8:306305/149‑1 (MQ=255)
gtcgtcCATGCGCGCTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGgccag                                                                                                                                                   >  6:263640/1‑149 (MQ=255)
 tcgtcCATGCGCGCTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCACCCTGCAGCAGGTCGAGGAAGTGCTCGATA‑CCAGgccagca                                                                                                                                                 <  5:119702/149‑1 (MQ=255)
   gtcCATGCGCGCTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGgccagcag                                                                                                                                                >  8:109312/1‑149 (MQ=255)
   gtcCATGCGCGCTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGgccagca                                                                                                                                                 <  7:105019/148‑1 (MQ=255)
       aTGCGCGGTTCGCTCTGGGCCAGCGGCGCCCCTGCTAGGTACAAGGGGCCTTGGATGGTGCGCGGGGTGCCGGCGCCAATGCCGGCTGCGGCATCGTTGGCATCCTGCAGCAGGTCGAGGATGTGCTCGATACCCAGgccagcagccag                                                                                                                                            <  6:379614/149‑1 (MQ=37)
        tGCGCGCTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGgccagcagccag                                                                                                                                            >  6:394855/1‑148 (MQ=255)
           gcgcTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGAGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAgg                                                                                                                                        >  8:282910/1‑149 (MQ=255)
            cgcTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGgc                                                                                                                                       <  5:358735/149‑1 (MQ=255)
                   cgcTGGGCCAGGGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGGTTCGGCATCCTTGGCATCCTGCAGCTGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTc                                                                                                                                <  4:373029/147‑1 (MQ=255)
                   cccTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTc                                                                                                                                >  4:138366/1‑149 (MQ=255)
                   cccTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTc                                                                                                                                >  8:193208/1‑149 (MQ=255)
                    ccTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCg                                                                                                                               >  1:403685/1‑149 (MQ=255)
                        ggCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTAc                                                                                                                           >  4:267085/1‑149 (MQ=255)
                         gCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACg                                                                                                                          <  7:99949/149‑1 (MQ=255)
                         gCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACg                                                                                                                          <  3:267085/149‑1 (MQ=255)
                          ccAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCTGCTGGCTGCTGgcctgg                                                                                                                                   >  4:222813/1‑138 (MQ=25)
                          ccAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCTGCTGGCTGCTGgcctgg                                                                                                                                   <  3:222813/139‑2 (MQ=25)
                          ccAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGgc                                                                                                                                       <  3:119976/135‑1 (MQ=255)
                          ccAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACgg                                                                                                                         >  7:96301/1‑149 (MQ=255)
                          ccAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTAGATACCCAGGCCAGCAGCCAGCAGgc                                                                                                                                       >  4:119976/1‑135 (MQ=255)
                                   ccccGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAgg                                                                                                                <  3:138366/149‑1 (MQ=255)
                                    ccccGCAACGTCCCACTGTCCTGCCATGGGGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCCTCCTGCAGCCGGTCGCGGAAGTGCTCGATCCCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGc                                                                                                               <  4:321415/149‑1 (MQ=255)
                                    cccGGCAACGTACAACGGGCCTTCTATGGTGCGCGGGGTGCCGCCGCCAAGGCAGGCTTCTGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGATCGATAACCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCTCCAAGGc                                                                                                               >  5:361888/1‑149 (MQ=37)
                                        gCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGtt                                                                                                           <  5:394855/149‑1 (MQ=255)
                                         cAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTg                                                                                                          <  5:190729/149‑1 (MQ=255)
                                          aaCGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGa                                                                                                         <  7:113946/149‑1 (MQ=255)
                                          aaCGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGa                                                                                                         <  2:116265/149‑1 (MQ=255)
                                           aCGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAg                                                                                                        <  5:256861/149‑1 (MQ=255)
                                                cAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAg                                                                                                    <  7:282910/148‑1 (MQ=255)
                                                   cGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGa                                                                                                <  7:63025/149‑1 (MQ=255)
                                                     ggCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACg                                                                                              >  5:184580/1‑149 (MQ=255)
                                                         ttCGATGTTGCGCGGGTTGCCGACGCCAAGGCCGGCTTCGGCATCCTTGGAATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCATCATGCCTGCCTCGTTACGGCCGCACAGGCGGTTGACGTAGTCGACGGCGt                                                                                          >  4:374096/1‑149 (MQ=37)
                                                           cGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGt                                                                                          <  2:403685/147‑1 (MQ=255)
                                                           cGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGc                                                                                        >  5:8357/1‑149 (MQ=255)
                                                            gATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGc                                                                                            >  1:292866/1‑144 (MQ=255)
                                                            gATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGc                                                                                            <  2:292866/144‑1 (MQ=255)
                                                            gATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGc                                                                                        <  5:263640/148‑1 (MQ=255)
                                                            gATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGc                                                                                        <  4:128362/148‑1 (MQ=255)
                                                            gATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGc                                                                                        >  3:125480/1‑148 (MQ=255)
                                                              tGGTGCGGGGGGTGCCGCCGCTAAGGCGGGCTTCGGAATCCTTGGCATCCTGGAGCAGTTCGAGGAAGTGCTCGATACCCATGCCTGCAGCCAGCATGCCTGCCTGGTTACTGCCGCCCTGGCGGTTGTTGTAGTCGACGGCGTTCCa                                                                                      <  6:361888/148‑1 (MQ=255)
                                                                gTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGaa                                                                                   <  3:138688/149‑1 (MQ=255)
                                                                 tGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAAc                                                                                  <  2:284351/149‑1 (MQ=255)
                                                                   cgcgGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTc                                                                                <  7:146726/149‑1 (MQ=255)
                                                                    gcgGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGTAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGCCGGCGTGCCAGCACTCg                                                                               >  8:264874/1‑149 (MQ=255)
                                                                    gcgGGGTGCCGCCGCAAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCg                                                                               <  3:374096/149‑1 (MQ=255)
                                                                        ggcgccgccgccAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCt                                                                           <  5:370848/146‑1 (MQ=255)
                                                                        ggTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCt                                                                           <  7:193208/149‑1 (MQ=255)
                                                                                    aaGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAAtctc                                                                 >  1:113559/1‑147 (MQ=255)
                                                                                               cGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGAt                                                     >  1:407518/1‑148 (MQ=255)
                                                                                                ggCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATc                                                    <  6:184580/148‑1 (MQ=255)
                                                                                                  cATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAgg                                                 <  8:292838/149‑1 (MQ=255)
                                                                                                      cTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAggcggg                                             <  8:96301/149‑1 (MQ=255)
                                                                                                      cTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAggcggg                                             >  1:305881/1‑149 (MQ=255)
                                                                                                          gCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGGCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGgt                                         <  7:109312/149‑1 (MQ=255)
                                                                                                             tgcTGCAGCTTGTCGAGGATGTGCTGGATACCCAGGCCAGCTGCCAGCAGGCCTGCCTGGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGTACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCTGTGTc                                      <  2:115304/147‑1 (MQ=25)
                                                                                                                 gcagcaGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTc                                                        >  1:239593/1‑127 (MQ=255)
                                                                                                                 gcagcaGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTc                                                        <  2:239593/127‑1 (MQ=255)
                                                                                                                  cagcagGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTgca                                 <  6:8357/149‑1 (MQ=255)
                                                                                                                         tCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCAcgcg                          >  5:315866/1‑149 (MQ=255)
                                                                                                                            aGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCAcgcg                          >  7:38448/1‑146 (MQ=255)
                                                                                                                                  tGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGAt                 <  1:276457/149‑1 (MQ=255)
                                                                                                                                     tctaTACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGAGTCTTGCAGCACGCGCAGAATGATCTg              >  7:265831/4‑149 (MQ=255)
                                                                                                                                       gATAGCCAGGCCAGCAGGCATCAGGCCTGCCTCGTTACTGCCTCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGc             <  4:172028/148‑1 (MQ=37)
                                                                                                                                       gATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCt            <  1:75510/149‑1 (MQ=255)
                                                                                                                                        cgaggcAGGCCTGCTGGCTGCTGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCtt           >  5:79380/6‑149 (MQ=25)
                                                                                                                                         tACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCATGTCTTCGATCAGGCGGTCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTg          <  8:370994/149‑1 (MQ=255)
                                                                                                                                           cccAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGaa        <  2:305881/149‑1 (MQ=255)
                                                                                                                                            ccAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAg       >  7:296022/1‑149 (MQ=255)
                                                                                                                                            ccAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAg       >  8:55734/1‑149 (MQ=255)
                                                                                                                                              aGGCCTGCTGGCTGCTGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAgcg     >  5:80612/1‑149 (MQ=25)
                                                                                                                                                gccagcagccagcaGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCtt                                                         >  4:118792/1‑95 (MQ=255)
                                                                                                                                                gccagcagccagcaGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCtt                                                         <  3:118792/95‑1 (MQ=255)
                                                                                                                                                gccagcagccagcaGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGcgcg   >  4:157776/1‑149 (MQ=255)
                                                                                                                                                 ccagcagccagcagGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGCGCgg  <  8:38448/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTCGTCCATGCGCGCTTCGCCCTGGGCCAGCGGCGCCCCGGCAACGTACAACGGGCCTTCGATGGTGCGCGGGGTGCCGCCGCCAAGGCCGGCTTCGGCATCCTTGGCATCCTGCAGCAGGTCGAGGAAGTGCTCGATACCCAGGCCAGCAGCCAGCAGGCCTGCCTCGTTACGGCCGCCCAGGCGGTTGAGGTAGTCGACGGCGTGCCAGAACTCGTCCTCGGTAATCTCCAGGTCTTCGATCAGGCGGGCGGTGTCTTGCAGCACGCGCAGAATGATCTGCTTGAAGCGCGG  >  NC_002947/4236398‑4236691

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