Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP012306406,6770AG24.6% 507.1 / 24.5 199K294E (AAA→GAA) hemHferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (27/123);  new base G (49/0);  total (76/123)
Fisher's exact test for biased strand distribution p-value = 2.51e-27
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 233 total aligned reads displayed.

TGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAA  >  AP012306/406528‑406826
                                                                                                                                                     |                                                                                                                                                       
  ccTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGgaa                                                                                                                                                       >  2:717393/1‑150 (MQ=255)
  ccTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGgaa                                                                                                                                                       >  1:2496051/1‑150 (MQ=255)
  ccTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGgaa                                                                                                                                                       >  2:100390/1‑150 (MQ=255)
  ccTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGaaa                                                                                                                                                       <  1:2426511/150‑1 (MQ=255)
  ccTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGaaa                                                                                                                                                       <  1:965431/150‑1 (MQ=255)
     tataCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGaaaaaa                                                                                                                                                    <  2:549819/150‑1 (MQ=255)
      ataCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAat                                                                                                                                                   <  2:2006444/150‑1 (MQ=255)
      ataCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAat                                                                                                                                                   <  1:2586342/150‑1 (MQ=255)
       taCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGaaaaaaaa                                                                                                                                                  >  2:3167370/1‑148 (MQ=255)
       taCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAata                                                                                                                                                  <  1:4252222/150‑1 (MQ=255)
        aCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGgaaaaaaaa                                                                                                                                                 >  1:1577181/1‑147 (MQ=255)
        aCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAaaat                                                                                                                                                 >  2:3061732/1‑150 (MQ=255)
        aCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAatat                                                                                                                                                 <  1:69322/150‑1 (MQ=255)
          cgacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGgaaaaaaaaga                                                                                                                                               >  2:122876/1‑145 (MQ=255)
          cgacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGa                                                                                                                                               <  1:2752020/150‑1 (MQ=255)
           gacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGaa                                                                                                                                              <  1:1135244/150‑1 (MQ=255)
              gaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAatat                                                                                                                                           <  2:16194/150‑1 (MQ=255)
              gaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAatat                                                                                                                                           <  2:2630357/150‑1 (MQ=255)
                 aCGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑Att                                                                                                                                       <  2:1474571/150‑1 (MQ=255)
                   gCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTcc                                                                                                                                     <  2:2197272/150‑1 (MQ=255)
                     tGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAAAGAAGAT‑ATTctcg                                                                                                                                   >  2:4099303/1‑147 (MQ=255)
                      gAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATATAATTCCgc                                                                                                                                   >  1:3232606/1‑149 (MQ=255)
                        aaaTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATTAAAAT‑AATCCGg‑ggc                                                                                                                               >  2:3661979/1‑150 (MQ=255)
                         aaTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATAT‑ATTTCCG‑CGCt                                                                                                                              >  1:3085425/1‑150 (MQ=255)
                          aTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCtt                                                                                                                             <  1:3846873/150‑1 (MQ=255)
                            gCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATAT‑ATTCCCG‑CGCttta                                                                                                                           >  1:3768188/1‑148 (MQ=255)
                              tCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATg                                                                                                                         <  1:3571215/150‑1 (MQ=255)
                               cGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATTAATAT‑ATTTCCG‑CGCCTAAttc                                                                                                                        >  1:1008171/1‑148 (MQ=255)
                               cGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATAT‑ATTCCGG‑CGCTTTATGc                                                                                                                        >  2:3044278/1‑150 (MQ=255)
                                 gagaAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTTCGG‑CGCTTAATGCCa                                                                                                                      >  1:2838083/1‑150 (MQ=255)
                                  agaAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAAAGAAAAA‑ATTTCGG‑CCCTTAATGcccc                                                                                                                     >  2:37352/1‑148 (MQ=255)
                                    aaaaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGc                                                                                                                   <  2:1932535/150‑1 (MQ=255)
                                     aaaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATTCCACGcc                                                                                                                  >  1:2779160/1‑150 (MQ=255)
                                     aaaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGcc                                                                                                                  <  2:354641/150‑1 (MQ=255)
                                     aaaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGcc                                                                                                                  <  2:3745893/150‑1 (MQ=255)
                                      aaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATAT‑ATTTCGG‑CGCTTTATGCCACCCCg                                                                                                                 >  1:1759604/1‑150 (MQ=255)
                                        aGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATAT‑ATTTCCG‑CGCTTAATTCCACGcccgg                                                                                                               >  1:2149139/1‑147 (MQ=255)
                                           cGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATTAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGGAcc                                                                                                            >  1:1362606/1‑149 (MQ=255)
                                            gTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAAAGAATAA‑AATTCTG‑CGCCTAAAGACACGCCGGGAcca                                                                                                           >  2:379855/1‑148 (MQ=255)
                                            gTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACat                                                                                                           <  1:2030110/150‑1 (MQ=255)
                                             tAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATAT‑ATTCCGG‑CGCTTAATTCCACGCCGGAAAata                                                                                                          >  2:734639/1‑150 (MQ=255)
                                             tAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACata                                                                                                          <  1:1410909/150‑1 (MQ=255)
                                              aGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAAAT‑ATTTCGG‑CGCTTTATTCCACGCCGGGACatat                                                                                                         >  2:2317355/1‑150 (MQ=255)
                                               ggTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATc                                                                                                        <  2:2159880/150‑1 (MQ=255)
                                               ggTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATc                                                                                                        <  2:1693090/150‑1 (MQ=255)
                                                gTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAAAAT‑ATTCCGG‑CGCTTTAAGCCACGCCGGGACATATCg                                                                                                       >  1:540521/1‑150 (MQ=255)
                                                 tCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAACAAATAAATAT‑TTTCCCG‑CCCCTTAAGCCACCCCCGAAAatattgc                                                                                                      >  2:3546541/1‑147 (MQ=255)
                                                    tatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAat                                                                                                   <  2:1576759/150‑1 (MQ=255)
                                                    tatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAat                                                                                                   <  2:3407752/150‑1 (MQ=255)
                                                    tatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAat                                                                                                   <  2:2207336/150‑1 (MQ=255)
                                                      ttCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGGTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAatga                                                                                                 <  1:4190246/150‑1 (MQ=255)
                                                      ttCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATTAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGGACACATCGAAaaga                                                                                                 >  2:2943214/1‑150 (MQ=255)
                                                      ttCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAatga                                                                                                 <  2:3102937/150‑1 (MQ=255)
                                                       tCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAatgat                                                                                                <  2:3303093/150‑1 (MQ=255)
                                                        cAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATAT‑ATTTCCG‑CGCTTAATGCCACGCCCGAACATATCGAAaagatg                                                                                               >  1:1465424/1‑150 (MQ=255)
                                                         aGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATAT‑ATTTCCG‑CGCTTAATACCACCCCCGGACATATCGAAATGATgg                                                                                              >  2:37075/1‑150 (MQ=255)
                                                         aGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATgg                                                                                              <  2:1125291/150‑1 (MQ=255)
                                                         aGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATgg                                                                                              <  2:2239567/150‑1 (MQ=255)
                                                          ggTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGc                                                                                             <  1:2968498/150‑1 (MQ=255)
                                                              aTGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAAt                                                                                         <  2:2080950/150‑1 (MQ=255)
                                                               tgtgCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATc                                                                                        <  1:3366890/150‑1 (MQ=255)
                                                                 tgCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCtt                                                                                      <  1:4381298/150‑1 (MQ=255)
                                                                    ccGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATAT‑ATTTCCG‑CGCTTAATGCCACGCCCGGACATATCGAAATGATGGCTAATCttgtt                                                                                   >  1:866569/1‑150 (MQ=255)
                                                                    ccGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAAAT‑ATTTCGG‑CGCTTAATACCCCGCCGGGACAAATCGAAAAGATGGCTAAACttttt                                                                                   >  2:3728210/1‑150 (MQ=255)
                                                                    ccGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCttgtt                                                                                   <  1:251788/150‑1 (MQ=255)
                                                                     cGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTTCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCttgttg                                                                                  >  1:3900034/1‑150 (MQ=255)
                                                                        gCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCcg                                                                               <  1:256274/150‑1 (MQ=255)
                                                                        gCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCcg                                                                               <  2:2830731/150‑1 (MQ=255)
                                                                            tgctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATTAATAT‑ATTTCCG‑CGCTTAAAGCCACCCCCGGACATATACAAAAAATGGCCAAACTTTTTGCCCCGTa                                                                           >  2:4199528/1‑150 (MQ=255)
                                                                              ctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATAA‑ATTCCGG‑CGCCTAATGCCCCGCCGGAACATATCGAAAAGATGGCTAATCTTGTTGCCGCGCATc                                                                         >  2:1795781/1‑150 (MQ=255)
                                                                              ctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATc                                                                         <  1:1255175/150‑1 (MQ=255)
                                                                               tgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCg                                                                        <  1:1813819/150‑1 (MQ=255)
                                                                                 cGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCt                                                                      <  2:485639/150‑1 (MQ=255)
                                                                                   gATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATAT‑ATTTCGG‑CGCTTAATGCCCCGCCGGGACATATCGAAATGATGGCTAATCTTGTTTGCGCGTATCGCTaa                                                                    >  2:3479886/1‑150 (MQ=255)
                                                                                       gTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATAT‑ATTTCGG‑CGCTTAATGCCACGCCGGGACATATCCAAAAGATGGCTAATCTTTTTTCCCCGGATCGCTTAAACt                                                                >  2:3923151/1‑150 (MQ=255)
                                                                                       gTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATAT‑ATTTCGG‑CGCCTTATTCCACCCCGGGACAAATCGAAAAGATGGGGAATCTTTTTTCCGCGTATCGCTAAAGCt                                                                >  2:945338/1‑150 (MQ=255)
                                                                                         cTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGa                                                              <  1:791355/150‑1 (MQ=255)
                                                                                          tGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAg                                                             <  2:2896318/150‑1 (MQ=255)
                                                                                            gagaCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCg                                                           <  1:4013033/150‑1 (MQ=255)
                                                                                                cGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTaa                                                       <  2:2334202/150‑1 (MQ=255)
                                                                                                    ggAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGaa                                                   <  1:287067/150‑1 (MQ=255)
                                                                                                                cGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAAAA‑ATTTCGG‑CGCTTAATGCCACGCCGGGACATATCGAAAAGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAAATGAGCGCCGTc                                       >  2:2431432/1‑150 (MQ=255)
                                                                                                                     aaaaCCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGaa                                  <  2:604233/150‑1 (MQ=255)
                                                                                                                       aaCCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAAAT‑ATTTCGG‑CGCTTAATGCCCCGCCGGGAAAAATCGAAAAGATGGCTAATTTTTTTGCCGCGTATCGCTAAAACTGGGCGGGAAAGAAATGAGGGCCGTCGCGaaaa                                >  2:1951709/1‑148 (MQ=255)
                                                                                                                        aCCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTTATGCCACGCCGGGACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAACTGAGCGCCGTCGCGAagag                               >  1:1559891/1‑150 (MQ=255)
                                                                                                                        aCCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAagag                               <  2:1888695/150‑1 (MQ=255)
                                                                                                                            tGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtc                           <  2:2856058/150‑1 (MQ=255)
                                                                                                                               ggTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgtc                        <  2:3896640/150‑1 (MQ=255)
                                                                                                                               ggTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgtc                        <  1:3174224/150‑1 (MQ=255)
                                                                                                                                gTCTTCCTCGGTGCCGGCGGGGAAAAAAAATAATAA‑ATTCCGGCCCCTTAAT‑CCACGCCGGGACATATTGAAAAGATGGCTTAACTTGTTGCCGCCCATTCCTAAAGCTGAGCGGTAAAAAAACGAGCGCCGGCCCGAACAGCCccgtcg                       >  2:2551363/1‑150 (MQ=255)
                                                                                                                                  cTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTcggc                     <  2:3946323/150‑1 (MQ=255)
                                                                                                                                  cTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTcggc                     <  2:356537/150‑1 (MQ=255)
                                                                                                                                  cTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTcggc                     <  1:2232464/150‑1 (MQ=255)
                                                                                                                                    tCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTcggcgg                   <  2:4119301/150‑1 (MQ=255)
                                                                                                                                         ggTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtc              <  2:3427332/150‑1 (MQ=255)
                                                                                                                                              cggcggGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTcatcat         <  1:1090560/150‑1 (MQ=255)
                                                                                                                                                gcggGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATcc       <  2:872360/150‑1 (MQ=255)
                                                                                                                                                  gggAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGt     >  1:3555841/1‑150 (MQ=255)
                                                                                                                                                     aaaaaaTATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGaa  >  2:556878/1‑150 (MQ=255)
                                                                                                                                                     aaaaaaTATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGaa  <  2:93386/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATAT‑ATTCCGG‑CGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAA  >  AP012306/406528‑406826

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