Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP012306406,6770AG22.2% 517.0 / 18.2 199K294E (AAA→GAA) hemHferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (37/117);  new base G (44/0);  total (82/117)
Fisher's exact test for biased strand distribution p-value = 6.57e-22
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 236 total aligned reads displayed.

TGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGA  >  AP012306/406528‑406825
                                                                                                                                                       |                                                                                                                                                     
     tataCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑gaaaaa                                                                                                                                                  >  2:2912113/1‑150 (MQ=255)
         ccGACGAAACGCTGAAAAAGCTCGGAAAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGGGGCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATg                                                                                                                                              <  2:832742/150‑1 (MQ=255)
          cgacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑aaaaaaaaaga                                                                                                                                             >  2:3688016/1‑145 (MQ=255)
           gacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑aaaaaaaaaaaa                                                                                                                                            >  2:1034507/1‑144 (MQ=255)
           gacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAAAATGaa                                                                                                                                            >  2:2408315/1‑150 (MQ=255)
            acgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAATAAGAat                                                                                                                                           >  1:2535222/1‑150 (MQ=255)
             cgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAATATTAata                                                                                                                                          >  1:1308826/1‑150 (MQ=255)
             cgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAATATGAata                                                                                                                                          >  1:956155/1‑150 (MQ=255)
             cgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAAAATGAata                                                                                                                                          >  1:1916125/1‑150 (MQ=255)
             cgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAata                                                                                                                                          >  1:3340821/1‑150 (MQ=255)
             cgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAata                                                                                                                                          <  2:1673604/150‑1 (MQ=255)
             cgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAata                                                                                                                                          <  2:910211/150‑1 (MQ=255)
              gaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAATATGaaaaa                                                                                                                                         >  2:3953013/1‑147 (MQ=255)
              gaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAAAATGAaaat                                                                                                                                         >  2:2649615/1‑150 (MQ=255)
              gaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AATAAATATGAatat                                                                                                                                         <  2:4276386/150‑1 (MQ=255)
              gaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAatat                                                                                                                                         <  2:3271817/150‑1 (MQ=255)
               aaaCGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAatat‑a                                                                                                                                       <  2:2580723/150‑1 (MQ=255)
               aaaCGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAatat‑a                                                                                                                                       <  1:1954676/150‑1 (MQ=255)
                aaCGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAatat‑at                                                                                                                                      <  2:2372001/150‑1 (MQ=255)
                    cTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCg                                                                                                                                  <  2:1394587/150‑1 (MQ=255)
                     tGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAATATAAATAT‑Atttcag                                                                                                                                 >  1:735845/1‑146 (MQ=255)
                       aaaaTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAATATGAATAT‑ATTCCcgcg                                                                                                                               >  1:817601/1‑150 (MQ=255)
                         aaTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCt                                                                                                                             <  1:1202983/150‑1 (MQ=255)
                          aTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCtt                                                                                                                            <  1:3564572/150‑1 (MQ=255)
                           tGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGG‑AAAAAAAATAAAAAA‑TTTACGGCGCtc                                                                                                                            >  2:2084508/1‑149 (MQ=255)
                           tGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAAAATTAATAT‑ATTTCCGCGCCTa                                                                                                                           >  1:101576/1‑150 (MQ=255)
                            gCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAAAATTAAAATAATTCGGGGGCTTa                                                                                                                           >  1:4354752/1‑150 (MQ=255)
                            gCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTaa                                                                                                                          <  2:2201971/150‑1 (MQ=255)
                             cTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAAt                                                                                                                         <  1:3194548/150‑1 (MQ=255)
                               cGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGc                                                                                                                       <  2:2664283/150‑1 (MQ=255)
                               cGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGc                                                                                                                       <  2:2707005/150‑1 (MQ=255)
                                ggAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAAAATGAATAT‑ATTCCGGCGCCTTATGcc                                                                                                                      >  1:2408376/1‑150 (MQ=255)
                                 gagaAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCa                                                                                                                     <  2:3792301/150‑1 (MQ=255)
                                      aaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAATATGAATAT‑ATTTCCGCGCTTAATTCCACGccc                                                                                                                >  1:2026499/1‑149 (MQ=255)
                                      aaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATTCCACGCCg                                                                                                                >  2:1066689/1‑150 (MQ=255)
                                       aaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCgg                                                                                                               <  1:1305808/150‑1 (MQ=255)
                                         ggCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAAAAAGAAAAA‑ATTTCCGCGCCTAATTCCCCCCCggga                                                                                                             >  2:1498631/1‑148 (MQ=255)
                                         ggCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTTCCGCGCCTAATGCCACGCCGGaa                                                                                                             >  1:3072387/1‑150 (MQ=255)
                                         ggCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGaa                                                                                                             <  1:2632427/150‑1 (MQ=255)
                                              aGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAATATTAATAT‑ATTTCCGCGCTTTATGCCACACCGGGACatat                                                                                                        >  1:3117713/1‑150 (MQ=255)
                                              aGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAAAATGAATAT‑ATTTCGGCGCTTAATGCCACGCCCGGACatat                                                                                                        >  1:994800/1‑150 (MQ=255)
                                              aGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAAAATGAATAA‑ATTCCGGCGCTTAATGCCACGCCGGAACatat                                                                                                        >  1:1542065/1‑150 (MQ=255)
                                              aGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACatat                                                                                                        <  1:399153/150‑1 (MQ=255)
                                               ggTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GGAAAATATGAATAT‑ATTCCGGCGCTTAATTCCCCGCCCGAACATATc                                                                                                       >  1:3439051/1‑150 (MQ=255)
                                               ggTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAAAATGAAAAT‑ATTTTGGGGCTTAATGCCACGCCGGGACATATc                                                                                                       >  2:3423731/1‑150 (MQ=255)
                                               ggTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGACTTCCTCGGTGCCGGC‑GGG‑AGAAAACAAGAAATA‑ATTCCCGCGCTTCATGCCCCCCCGGAAcccatc                                                                                                       >  2:1762316/1‑145 (MQ=255)
                                                  cATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGaa                                                                                                    <  1:3856180/150‑1 (MQ=255)
                                                   atatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGaaa                                                                                                   <  2:531680/150‑1 (MQ=255)
                                                    tatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTTCCGCGCTTAATGCCACGCCGGGACATATCGAAat                                                                                                  >  1:458180/1‑150 (MQ=255)
                                                     atTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAAAATGAAAAT‑ATTTCGGCGCTTAAAGCCCCGCCGGGACATATCGAAatt                                                                                                 >  2:334712/1‑149 (MQ=255)
                                                      ttCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATAAGAATAT‑ATTACCGCGCCTAATGCCACGCCCGGGCATATCTTAatga                                                                                                >  1:1638201/1‑150 (MQ=255)
                                                       tCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAAAATTAATAT‑ATTTCCGCGCTTAATGCCACCCCGGAACATATCCGAatgat                                                                                               >  1:3968497/1‑150 (MQ=255)
                                                         aGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCCGCGCCTCATGCCACGCCGGAACATATCGAAATGATgg                                                                                             >  1:2998614/1‑150 (MQ=255)
                                                            tGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTa                                                                                          <  1:3722852/150‑1 (MQ=255)
                                                              aTGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAAt                                                                                        <  2:1900407/150‑1 (MQ=255)
                                                                 tgCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCtt                                                                                     <  2:2822050/150‑1 (MQ=255)
                                                                  gCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCttg                                                                                    <  1:4210729/150‑1 (MQ=255)
                                                                  gCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCttg                                                                                    <  1:3348687/150‑1 (MQ=255)
                                                                   cccGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCttgt                                                                                   <  2:3099071/150‑1 (MQ=255)
                                                                    ccGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCttgtt                                                                                  <  1:3314802/150‑1 (MQ=255)
                                                                        gCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAATATGAATAT‑ATTTCGGCGCCTAATGCCACGCCGGGACATATCGAAATGATGGCTAATCTTGTTGCcg                                                                              >  2:632338/1‑150 (MQ=255)
                                                                           ttGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAAAATGAATAT‑ATTTCGGCGCTTTATGCCACGCCGGGACATATCGAAAAGATGGCTAATCTTTTTTCCGCGt                                                                           >  2:4201754/1‑150 (MQ=255)
                                                                             gctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAAAATGAATAT‑ATTTCCGCGGCTCATTACCCGCCGGGACATATCGAAATGATAGGTAAACTTTTTGCCGCGTAt                                                                         >  2:2122445/1‑150 (MQ=255)
                                                                               tgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCg                                                                       <  2:350883/150‑1 (MQ=255)
                                                                                gcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAAAATGAATAT‑ATTTCGGCGCTTAATGCCACGCCGGGAAAAAACGAAATGATGGCTAATCTTGTTGCCGCGGATCGc                                                                      >  2:6850/1‑150 (MQ=255)
                                                                                  ggATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTTCGGCGCTTAATTCCACGCCGGAACATATCGAAAAGATGGCTAATCTTGTTGTCGCGGATCGCTa                                                                    >  2:321833/1‑150 (MQ=255)
                                                                                         cTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGa                                                             <  2:562847/150‑1 (MQ=255)
                                                                                          tGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAAAATGAATAT‑ATTCCGGCGCTTAATGCCACCCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGaa                                                            >  1:1327880/1‑149 (MQ=255)
                                                                                           ggAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGc                                                           <  2:1602457/150‑1 (MQ=255)
                                                                                            gagaCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCg                                                          <  1:4150425/150‑1 (MQ=255)
                                                                                             agaCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑GAAAAAAATGAATAT‑ATTTCGGCGCTTAATTCCACCCCGGGACAAATCGAAATGATGGGTAAACTTTTTTCCGCGGATCGCCAAAACTGAGCgg                                                         >  2:2002789/1‑150 (MQ=255)
                                                                                                cGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACTTATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTaa                                                      <  1:269465/150‑1 (MQ=255)
                                                                                                 gCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTaaa                                                     <  1:1874543/150‑1 (MQ=255)
                                                                                                 gCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTaaa                                                     <  1:2773415/150‑1 (MQ=255)
                                                                                                    ggAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGaa                                                  <  1:2573982/150‑1 (MQ=255)
                                                                                                     gAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTTCGGCGCTTAATGCCACCCCCGGACATATCGGAATAATGGCTAATCTTGTTGCCGCGTATCGCTAAAACTGCGCGGTAAAGAAc                                                 >  2:2258515/1‑150 (MQ=255)
                                                                                                      aaGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAAAATGAATAT‑ATTCCGGCGCTTGATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGGATCGCTAAAGCTGAGCGGGAAAGaaat                                                >  2:3500628/1‑148 (MQ=255)
                                                                                                       agagATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTg                                               <  1:1980017/150‑1 (MQ=255)
                                                                                                         agaTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAAAATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAACTGAg                                             >  1:1211678/1‑150 (MQ=255)
                                                                                                          gaTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAgc                                            <  2:3047962/150‑1 (MQ=255)
                                                                                                            ttGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAgcgc                                          <  1:2503851/150‑1 (MQ=255)
                                                                                                              gCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCg                                        <  2:740726/150‑1 (MQ=255)
                                                                                                               ccGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGt                                       <  2:1179150/150‑1 (MQ=255)
                                                                                                                    cAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGa                                  <  2:3265724/150‑1 (MQ=255)
                                                                                                                            tGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtc                          <  1:1929564/150‑1 (MQ=255)
                                                                                                                               ggTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgtc                       <  2:1924534/150‑1 (MQ=255)
                                                                                                                                 tCTTCCTCGGTGCCGGC‑GGG‑AAAAAAAATGAATAT‑ATTTCGGCGCTTAATGCCCCGCCGGGACATATCGAAATGATGGCTAATCTTTTTGCCCCGTATCGCTAAAACTGAGCGGGAAAGAACTGAGCGCCGTCGTGAAAAGCCTCGTcgg                     >  2:4002728/1‑150 (MQ=255)
                                                                                                                                 tCTTCCTCGGTGCCGGC‑GGG‑AAAAAAAATGAAAAT‑ATTTCGGCGCTTAATTCCCCGCCGGGAAATATCGAAAAGATGGCTAATCTTTTTGCCGCGGATCGCCAAAACTGGGCGGGAAAAAAATGAGGGCCGTCGCGAAGAGccccttccg                     >  2:4117467/1‑143 (MQ=255)
                                                                                                                                     ccTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGt                 <  1:2911306/150‑1 (MQ=255)
                                                                                                                                     ccTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGt                 <  1:1554142/150‑1 (MQ=255)
                                                                                                                                       tCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGttt               <  1:3828061/150‑1 (MQ=255)
                                                                                                                                        cGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGtttt              <  2:460739/150‑1 (MQ=255)
                                                                                                                                          gTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtca            <  2:4039633/150‑1 (MQ=255)
                                                                                                                                            gCCGGC‑GGG‑GAAAAAAATGAATAT‑ATTTCGGCGCTTAATTCCCCCCCGGGACATATCGAAATGATGGCTAATCTTGTTGCCGCGGATCGCTAAAGCTGAGCGGGAAAGAACTGAGCGCCGTCGCGAAGAGGCTCGTCGGCGGgttttttc          >  2:2458561/1‑146 (MQ=255)
                                                                                                                                             ccGGG‑GGGAAAAAAATATAAATAT‑TTCCGGGCGTTAAATCCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcatc          <  1:2427875/150‑1 (MQ=255)
                                                                                                                                              cggc‑ggG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTcatcat        <  1:4213627/150‑1 (MQ=255)
                                                                                                                                              cggc‑ggG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGGTTTcatcat        >  1:2358703/1‑150 (MQ=255)
                                                                                                                                                gc‑ggG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATcc      <  1:3618298/150‑1 (MQ=255)
                                                                                                                                                gc‑ggG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATcc      >  1:1710261/1‑150 (MQ=255)
                                                                                                                                                   ggg‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGt    >  1:2134242/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
TGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGC‑GGG‑AAAAAATATGAATAT‑ATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGA  >  AP012306/406528‑406825

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