Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP012306406,6770AG24.0% 552.1 / 24.7 221K294E (AAA→GAA) hemHferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (43/125);  new base G (53/0);  total (96/125)
Fisher's exact test for biased strand distribution p-value = 8.22e-25
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 251 total aligned reads displayed.

TGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAA  >  AP012306/406528‑406826
                                                                                                                                                     |                                                                                                                                                      
  ccTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGgaa                                                                                                                                                      >  1:4557276/1‑150 (MQ=255)
   cTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGgaaa                                                                                                                                                     >  2:3873067/1‑150 (MQ=255)
    ttATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGcaaaa                                                                                                                                                    >  1:1985772/1‑150 (MQ=255)
    ttATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGaaaaa                                                                                                                                                    <  2:3310164/150‑1 (MQ=255)
     tataCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGgaaaaa                                                                                                                                                   >  1:3488103/1‑150 (MQ=255)
      ataCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGgaaaaaa                                                                                                                                                  >  1:3303518/1‑149 (MQ=255)
         ccGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATg                                                                                                                                               <  1:4534824/150‑1 (MQ=255)
          cgacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGTCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGa                                                                                                                                              >  1:1022891/1‑150 (MQ=255)
          cgacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAatatta                                                                                                                                              >  1:3452513/1‑148 (MQ=255)
          cgacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGa                                                                                                                                              >  1:668604/1‑150 (MQ=255)
          cgacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGa                                                                                                                                              <  2:316498/150‑1 (MQ=255)
           gacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATaaa                                                                                                                                             >  2:1043709/1‑150 (MQ=255)
           gacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGaa                                                                                                                                             >  2:291092/1‑150 (MQ=255)
             cgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAata                                                                                                                                           <  1:1274476/150‑1 (MQ=255)
              gaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAatat                                                                                                                                          >  2:2891085/1‑150 (MQ=255)
              gaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAatat                                                                                                                                          >  1:3137534/1‑150 (MQ=255)
               aaaCGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAatata                                                                                                                                         >  1:3232479/1‑150 (MQ=255)
                    cTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAAAGAAAATAtttcc                                                                                                                                    >  2:639618/1‑147 (MQ=255)
                     tGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCgg                                                                                                                                   <  1:4314010/150‑1 (MQ=255)
                      gAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGTAAAAATTCCGgc                                                                                                                                  >  1:4049931/1‑150 (MQ=255)
                           tGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTa                                                                                                                             <  1:1006089/150‑1 (MQ=255)
                               cGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCGGCGCTTAA‑TGc                                                                                                                        >  1:1392431/1‑150 (MQ=255)
                               cGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGc                                                                                                                        <  1:3068661/150‑1 (MQ=255)
                                    aaaaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATATATTTCCGCGCTTAA‑TGCCACGc                                                                                                                   >  1:4427311/1‑150 (MQ=255)
                                     aaaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCCGCGCTTAA‑TTCCCCGcc                                                                                                                  >  1:3618589/1‑150 (MQ=255)
                                     aaaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGCAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGcc                                                                                                                  <  1:4269825/150‑1 (MQ=255)
                                      aaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAAAGAAAAAATTTCCGCGCTTAA‑TTCCAccccc                                                                                                                 >  2:36769/1‑149 (MQ=255)
                                      aaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCg                                                                                                                 <  1:2579585/150‑1 (MQ=255)
                                       aaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAATATATTTCCGCGCTTAA‑TGCCACGcccg                                                                                                                >  1:2681645/1‑148 (MQ=255)
                                       aaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATAAATATAATACGGCGCTTAA‑TGGCACGCCgg                                                                                                                >  2:2589599/1‑150 (MQ=255)
                                        aGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATATATTTCCGCGCTTAA‑TGCCCCGCCggg                                                                                                               >  2:3520051/1‑149 (MQ=255)
                                         ggCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCTGGCGGGGAAAAAAATGAATATATTTCGGCGCCTAA‑TGGCACGCCggga                                                                                                              >  2:2369863/1‑148 (MQ=255)
                                         ggCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGGAAAAAATAAAATTATTACGGCGCCTAA‑TTCCACGcccgga                                                                                                              >  2:2138767/1‑146 (MQ=255)
                                         ggCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGaa                                                                                                              <  2:4526737/150‑1 (MQ=255)
                                            gTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAAAAATTCCGGCGCTTAATTGCCCCCCCGGacaa                                                                                                            >  2:2666463/1‑147 (MQ=255)
                                            gTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACat                                                                                                           <  2:869128/150‑1 (MQ=255)
                                             tAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACata                                                                                                          <  1:1375455/150‑1 (MQ=255)
                                              aGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACatat                                                                                                         <  2:1572548/150‑1 (MQ=255)
                                                gTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCg                                                                                                       <  1:1749382/150‑1 (MQ=255)
                                                  cATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGaa                                                                                                     <  1:3032725/150‑1 (MQ=255)
                                                   atatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATTAATATATTTCCGCGCCTAA‑TTCCACCCCCGGACAAAACGaaa                                                                                                    >  1:1103611/1‑150 (MQ=255)
                                                   atatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATAAAAAAATTTCCGCGCTTAA‑TACCACCCCCGGACAAATCCaaa                                                                                                    >  2:1250042/1‑150 (MQ=255)
                                                    tatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCCCGCCCGGACATATCGAAat                                                                                                   >  1:1316372/1‑150 (MQ=255)
                                                    tatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAat                                                                                                   <  1:2339032/150‑1 (MQ=255)
                                                    tatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAAATAATTTATTCCCGCGCCTTA‑TTCCCCCCCGGGAAATATctgaat                                                                                                   >  2:4238859/1‑145 (MQ=255)
                                                     atTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAAATAATAAATTTCGGCGGTTAA‑TACCACGCCCGGACAAATCGaaaag                                                                                                  >  2:705875/1‑148 (MQ=255)
                                                     atTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCCGCGCCTAA‑TTCCACGCCGGGACATATCGAAata                                                                                                  >  1:3435146/1‑149 (MQ=255)
                                                        cAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAAAGAAAAAATTTCCGCGCCTAA‑TGCCACGCCGGGACATAACGAAaagatg                                                                                               >  2:4349513/1‑150 (MQ=255)
                                                        cAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAatgatg                                                                                               <  1:1218859/150‑1 (MQ=255)
                                                        cAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAatgatg                                                                                               <  2:2104343/150‑1 (MQ=255)
                                                        cAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAatgatg                                                                                               <  2:1287384/150‑1 (MQ=255)
                                                         aGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAAAATATTTCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAAAGATgg                                                                                              >  1:231281/1‑150 (MQ=255)
                                                           gTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTTCCGCGCTTAA‑TGCCACGCCGGGACATATCGAAAAGAAGGCt                                                                                            >  1:639667/1‑150 (MQ=255)
                                                            tGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTa                                                                                           <  1:3846626/150‑1 (MQ=255)
                                                               tgtgCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATTAATATATTTCCGCGCTTTA‑TGCCACGCCCGAACATATCGAAAAGATTGCTAATc                                                                                        >  1:2066399/1‑150 (MQ=255)
                                                               tgtgCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATc                                                                                        <  2:4296209/150‑1 (MQ=255)
                                                                 tgCACGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCtt                                                                                      <  1:2187525/150‑1 (MQ=255)
                                                                   cccGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCttgt                                                                                    <  2:812864/150‑1 (MQ=255)
                                                                     cGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCttgttg                                                                                  <  2:4570318/150‑1 (MQ=255)
                                                                       ggCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAAATATTTCTGGGCTTTA‑TTCCACACCCGGACATATCGAAAATATGGGTTAACttgttttc                                                                                >  2:3438324/1‑147 (MQ=255)
                                                                        gCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCcg                                                                               <  1:1882205/150‑1 (MQ=255)
                                                                            tgctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATTAATATATTCCGGAGCGTTA‑TTCCCCGCCGGAACATAATGAAATGATGGCTAATCTTGTTGTCGCGTa                                                                           >  2:2400399/1‑150 (MQ=255)
                                                                            tgctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAAATATTTCCGCGCTTAA‑ATCCCCCCCGGGACATAACGAAAAAATGGGTAAACTTTTTGTcgcgga                                                                           >  2:3571389/1‑148 (MQ=255)
                                                                             gctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTAt                                                                          <  1:2824969/150‑1 (MQ=255)
                                                                               tgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCg                                                                        <  1:3586522/150‑1 (MQ=255)
                                                                                    aTTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTaaa                                                                   <  1:4166890/150‑1 (MQ=255)
                                                                                    aTTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTaaa                                                                   <  2:1998650/150‑1 (MQ=255)
                                                                                     ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAg                                                                  <  2:3079042/150‑1 (MQ=255)
                                                                                       gTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATTAAAATATTCCCGCGCTTAA‑TTCCCCCCCGGGACAAATCGAAAAAATGGGCAAACTTTTTGCCCCGTATCGCTAAAACt                                                                >  2:1752962/1‑150 (MQ=255)
                                                                                       gTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCt                                                                <  1:440677/150‑1 (MQ=255)
                                                                                        tCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTg                                                               <  1:2452571/150‑1 (MQ=255)
                                                                                           ggAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGc                                                            <  1:2762253/150‑1 (MQ=255)
                                                                                                cGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATATATTTCCGCGCTTAA‑TGCCACGCCCGAACATATCGGAATGATGGGTAATCTTGTTGCCGCGGATCGCTAAAGCTGAGCGGGaa                                                       >  2:1685251/1‑150 (MQ=255)
                                                                                                cGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTaa                                                       <  2:2425075/150‑1 (MQ=255)
                                                                                                 gCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATATATTCCGGCGCTTAA‑TTCCACGCCGGGACATATCGAAATGATGGCTAATCTTTTTTCCGCGTATCGCTAAAACTGAGCGGGaaa                                                      >  2:635310/1‑150 (MQ=255)
                                                                                                   tGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGa                                                    <  2:2288952/150‑1 (MQ=255)
                                                                                                      aaGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACt                                                 <  2:2526857/150‑1 (MQ=255)
                                                                                                         agaTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAg                                              <  2:2696457/150‑1 (MQ=255)
                                                                                                          gaTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAgc                                             <  1:487830/150‑1 (MQ=255)
                                                                                                                cGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTc                                       <  1:847931/150‑1 (MQ=255)
                                                                                                                      aaaCCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAAAAAATTTCGGCGCCTTA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATATTTTTTGCGCGTATCGCTAAAACTGAGCGGGAAAAAAATGAGCGCCCTCGCGAag                                 >  2:3766018/1‑150 (MQ=255)
                                                                                                                        aCCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAagag                               <  1:2048711/150‑1 (MQ=255)
                                                                                                                            tGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtc                           <  2:554111/150‑1 (MQ=255)
                                                                                                                              aGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgt                         <  1:71155/150‑1 (MQ=255)
                                                                                                                              aGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAAAATATTCCGGCGCTTTA‑TGCCACCCCGGGACAAATCGAAATGATGGCTAATCTTTTTGCCCCGTATCGCTAAAGCTGAGCGGGAAAAAACTGAGCGCCGTCGCGAAGAGCCtcgt                         >  2:1375605/1‑150 (MQ=255)
                                                                                                                                 tCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTcgg                      <  1:524376/150‑1 (MQ=255)
                                                                                                                                       tCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGttt                <  1:3860260/150‑1 (MQ=255)
                                                                                                                                        cGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGtttt               <  1:402359/150‑1 (MQ=255)
                                                                                                                                        cGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGtttt               <  2:3335072/150‑1 (MQ=255)
                                                                                                                                            gCCGGCGGGAAAAAAAATGAATATATTTCGGCGCTTAA‑TGCCACCCCGGGACATATCGAAAAGATGGCTAATCTTTTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGgttttttt           >  2:4476886/1‑146 (MQ=255)
                                                                                                                                             ccGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcatca          <  1:4189544/150‑1 (MQ=255)
                                                                                                                                               ggcggGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCatcatc        <  1:2980251/150‑1 (MQ=255)
                                                                                                                                               ggcggGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCatcatc        <  1:2204640/150‑1 (MQ=255)
                                                                                                                                                gcggGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATcc       <  2:4134691/150‑1 (MQ=255)
                                                                                                                                                gcggGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATcc       <  2:3731829/150‑1 (MQ=255)
                                                                                                                                                 cggGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCg      >  1:1597301/1‑150 (MQ=255)
                                                                                                                                                 cggGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATAATCCg      >  2:2227974/1‑150 (MQ=255)
                                                                                                                                                   ggAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTg    <  2:3945576/150‑1 (MQ=255)
                                                                                                                                                    gAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGa   >  2:1068225/1‑150 (MQ=255)
                                                                                                                                                     aaaaaaTATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGaa  <  1:905671/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAA‑TGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAA  >  AP012306/406528‑406826

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑