Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP012306406,6770AG23.8% 602.3 / 24.4 239K294E (AAA→GAA) hemHferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (39/143);  new base G (57/0);  total (96/143)
Fisher's exact test for biased strand distribution p-value = 1.90e-29
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 266 total aligned reads displayed.

TGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAA  >  AP012306/406528‑406826
                                                                                                                                                     |                                                                                                                                                      
tGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGa                                                                                                                                                        <  2:942589/150‑1 (MQ=255)
tGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGa                                                                                                                                                        <  1:3256473/150‑1 (MQ=255)
    ttATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGgaaaa                                                                                                                                                    >  1:3175160/1‑150 (MQ=255)
      ataCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAat                                                                                                                                                  <  1:133/150‑1 (MQ=255)
      ataCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAat                                                                                                                                                  <  1:2125055/150‑1 (MQ=255)
      ataCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCCGCGGGAAAAAat                                                                                                                                                  <  2:4207843/150‑1 (MQ=255)
        aCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAatat                                                                                                                                                <  2:2226966/150‑1 (MQ=255)
        aCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAatat                                                                                                                                                <  1:1507582/150‑1 (MQ=255)
         ccGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATg                                                                                                                                               >  1:3701810/1‑150 (MQ=255)
           gacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATTaa                                                                                                                                             >  1:2655559/1‑150 (MQ=255)
           gacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGaa                                                                                                                                             <  1:4207462/150‑1 (MQ=255)
            acgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAat                                                                                                                                            <  2:2521789/150‑1 (MQ=255)
            acgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAat                                                                                                                                            <  1:2557507/150‑1 (MQ=255)
            acgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATTAat                                                                                                                                            >  1:946224/1‑150 (MQ=255)
              gaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAata‑t                                                                                                                                         <  1:3878813/150‑1 (MQ=255)
              gaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAAAGAata‑t                                                                                                                                         >  2:1587809/1‑150 (MQ=255)
               aaaCGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAaaa‑ta                                                                                                                                        >  1:4014672/1‑150 (MQ=255)
                aaCGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAaaa‑tat                                                                                                                                       >  1:1970244/1‑150 (MQ=255)
                 aCGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAATAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGCAAA‑TAtt                                                                                                                                      >  1:4071540/1‑150 (MQ=255)
                 aCGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGCATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATA‑TAtt                                                                                                                                      >  2:3258612/1‑150 (MQ=255)
                  cGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATA‑TAttt                                                                                                                                     >  2:603536/1‑149 (MQ=255)
                    cTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATA‑TATTctg                                                                                                                                   >  2:24945/1‑148 (MQ=255)
                     tGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATA‑TATTTCgg                                                                                                                                  >  2:535847/1‑150 (MQ=255)
                     tGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTcccg                                                                                                                                  >  1:23279/1‑148 (MQ=255)
                      gAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGgc                                                                                                                                 <  1:1210690/150‑1 (MQ=255)
                       aaaaTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATA‑TATTTCGgcg                                                                                                                                >  1:2140580/1‑150 (MQ=255)
                       aaaaTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTTCGgcg                                                                                                                                >  1:1958401/1‑150 (MQ=255)
                        aaaTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATTAAAA‑TATTCCGgcgc                                                                                                                               >  1:4184559/1‑150 (MQ=255)
                            gCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATA‑TATTCCGGCGCTTaa                                                                                                                           >  1:1430438/1‑150 (MQ=255)
                            gCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAGAAAATATGAATA‑TATTCCGGCGCTTaa                                                                                                                           <  2:643424/150‑1 (MQ=255)
                            gCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTaa                                                                                                                           <  2:378021/150‑1 (MQ=255)
                             cTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATA‑TATTTCGCCGCTTaaa                                                                                                                          >  1:3022585/1‑149 (MQ=255)
                              tCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAAA‑TATTCCGGCGCTTAATg                                                                                                                         >  1:3573290/1‑150 (MQ=255)
                              tCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATg                                                                                                                         <  2:179443/150‑1 (MQ=255)
                              tCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATg                                                                                                                         <  1:3739209/150‑1 (MQ=255)
                                ggAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGGAAAAAAAGAAAA‑TATTTCGGCCGTTAATGcc                                                                                                                       >  2:3008834/1‑150 (MQ=255)
                                 gagaAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAAATAATATTTTTCCGGGGCtaatggcc                                                                                                                       >  2:2430172/1‑143 (MQ=255)
                                    aaaaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATA‑TATTTCGGCGCTTAATGcccccc                                                                                                                   >  2:2826093/1‑146 (MQ=255)
                                    aaaaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAAAA‑TATTTCCGCGCCTAATTCCAccc                                                                                                                   >  2:1456503/1‑148 (MQ=255)
                                     aaaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGcc                                                                                                                  <  1:4135844/150‑1 (MQ=255)
                                        aGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGa                                                                                                               <  2:77971/150‑1 (MQ=255)
                                          gCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAAc                                                                                                             <  2:4120423/150‑1 (MQ=255)
                                           cGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACa                                                                                                            <  2:3348317/150‑1 (MQ=255)
                                            gTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATTAATA‑TATTTCGGCGCGTTATGCTACGCCGGAAcca                                                                                                           >  1:970994/1‑148 (MQ=255)
                                            gTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACat                                                                                                           <  1:3113307/150‑1 (MQ=255)
                                            gTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACat                                                                                                           <  2:876439/150‑1 (MQ=255)
                                             tAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTTTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACata                                                                                                          <  2:395919/150‑1 (MQ=255)
                                             tAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACata                                                                                                          <  1:2193343/150‑1 (MQ=255)
                                               ggTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATc                                                                                                        <  2:3650493/150‑1 (MQ=255)
                                               ggTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATc                                                                                                        <  1:1711578/150‑1 (MQ=255)
                                                gTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATTAATA‑TATTCCCGCGCCTTATTCCACGCCGGAACATATCg                                                                                                       >  1:1003134/1‑150 (MQ=255)
                                                gTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCg                                                                                                       <  1:2786055/150‑1 (MQ=255)
                                                    tatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAat                                                                                                   <  1:1055827/150‑1 (MQ=255)
                                                    tatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAat                                                                                                   <  1:3951416/150‑1 (MQ=255)
                                                      ttCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAatga                                                                                                 <  1:878707/150‑1 (MQ=255)
                                                       tCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAatgat                                                                                                <  1:720271/150‑1 (MQ=255)
                                                         aGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATgg                                                                                              <  2:3082010/150‑1 (MQ=255)
                                                          ggTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATA‑TATTTCGGCGCTTAATTCCCCGCCGGAACATATCGAAATGATGGc                                                                                             >  1:2835391/1‑150 (MQ=255)
                                                          ggTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAAA‑TATTCCGGCGCTTTATGCCACGCCCGGACATATTGGAAAGAAGGc                                                                                             >  2:803218/1‑150 (MQ=255)
                                                            tGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATA‑TATTCCCGCGCTTAATTCCACCCCGGGACATATCTAAAAGATGGGTa                                                                                           >  1:517570/1‑150 (MQ=255)
                                                              aTGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGGACATATCGAAAAGATGGCTAAt                                                                                         >  1:2275783/1‑150 (MQ=255)
                                                                   cccGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCttgt                                                                                    <  2:3412582/150‑1 (MQ=255)
                                                                    ccGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCttgtt                                                                                   <  2:1110965/150‑1 (MQ=255)
                                                                       ggCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGcc                                                                                <  1:449213/150‑1 (MQ=255)
                                                                         cTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAAA‑TATTTCGGCGCTTAATGCCACGCCGGGACATATCGAAATGATGGGTAATCTTTTTGCcgc                                                                              >  2:464233/1‑150 (MQ=255)
                                                                         cTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGATTA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCcgc                                                                              <  2:2162789/150‑1 (MQ=255)
                                                                         cTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCcgc                                                                              <  1:361913/150‑1 (MQ=255)
                                                                          tttGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATTAAAA‑TATTTCCGCGCTTAATGCCACCCCCGAACAAATTGAAATGATGGCTTAACTTTTTTCcgcg                                                                             >  2:123237/1‑150 (MQ=255)
                                                                          tttGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATA‑TATTTCGGCGGTTAATTCCACGCCGGGACATATCGAAAAGATGGCTAATCTTTTTGCcgcg                                                                             >  2:1152186/1‑150 (MQ=255)
                                                                            tgctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTa                                                                           <  1:3869381/150‑1 (MQ=255)
                                                                             gctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTAt                                                                          <  1:2630088/150‑1 (MQ=255)
                                                                              ctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATc                                                                         <  1:3977675/150‑1 (MQ=255)
                                                                               tgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCg                                                                        <  1:2867982/150‑1 (MQ=255)
                                                                                   gATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAAA‑TATTTCGGCGCTTAATGCCACCCCGGGACATAACGAAAAGATGGCTAAACTTGTTGCCGCGTATCGCTaa                                                                    >  2:677027/1‑150 (MQ=255)
                                                                                     ttGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAg                                                                  <  1:1130155/150‑1 (MQ=255)
                                                                                      tGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGc                                                                 <  1:3462511/150‑1 (MQ=255)
                                                                                        tCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTg                                                               <  1:3794460/150‑1 (MQ=255)
                                                                                         cTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATA‑TATTCCCGCGCTTAATTCCACGCCGGGACAAAACCAAATGATGGCTAATCTTTTTGGCCCGTATCCCTAAAACttg                                                              >  2:2644866/1‑148 (MQ=255)
                                                                                          tGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCTGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGGATCGCTAAAGCTGAg                                                             >  1:1707702/1‑150 (MQ=255)
                                                                                              gaCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTAGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGt                                                         <  1:4095343/150‑1 (MQ=255)
                                                                                                 gCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTaaa                                                      <  1:3128812/150‑1 (MQ=255)
                                                                                                  cTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAg                                                     <  1:157248/150‑1 (MQ=255)
                                                                                                     gAAGAGATTGACGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAAc                                                  <  1:4238896/150‑1 (MQ=255)
                                                                                                         agaTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATTAAAA‑TTTTTCGGCGCTTAAATCCACGGCGGGAAAAATCGAAAAGATGGGGAATCTTGTTGCCGCGGATTGCTAAAACTGAACGGGAAAGaaattgg                                              >  2:43728/1‑145 (MQ=255)
                                                                                                                    cAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCCCGCCGGGACATATCGAAAAGATGGGTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAAATGAGCGCCGTCGCGa                                   >  2:1702788/1‑150 (MQ=255)
                                                                                                                       atccGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAaga                                <  1:3860386/148‑1 (MQ=255)
                                                                                                                        aCCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATTAATA‑TATTCCGGCGCTTAATGCCACGCCGGGACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAACTGAGCGCCGTCGCGAagag                               >  1:2347217/1‑150 (MQ=255)
                                                                                                                          cGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGcc                             <  2:2889220/150‑1 (MQ=255)
                                                                                                                           gTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCt                            <  1:2382567/150‑1 (MQ=255)
                                                                                                                            tGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtc                           <  1:3363214/150‑1 (MQ=255)
                                                                                                                            tGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtc                           <  1:3489865/150‑1 (MQ=255)
                                                                                                                               ggTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTGTCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCGGTCGCGAAGAGCCtcgtc                        <  2:1815230/150‑1 (MQ=255)
                                                                                                                                     ccTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGt                  <  2:548484/150‑1 (MQ=255)
                                                                                                                                     ccTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGt                  <  2:2286257/150‑1 (MQ=255)
                                                                                                                                        cGGTGCCGGCGGGAAAAAAAAAAAATA‑AATTCCGGCGCTTAATGCCACGCCGGGACATATCGAAATGATGGCTAATCTTTTTGCCGCGGATCGCTAAAGCTGAGCGGGAAAAAACTGAGCGCCGTCGCGGAGAGCCTCGTcggcggggtt               >  2:3953848/1‑146 (MQ=255)
                                                                                                                                             ccGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcatca          <  2:1897711/150‑1 (MQ=255)
                                                                                                                                                 cggGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTTATTATCCg      >  2:2037108/1‑150 (MQ=255)
                                                                                                                                                  gggAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGt     <  2:3578494/150‑1 (MQ=255)
                                                                                                                                                  gggAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGt     >  1:835232/1‑150 (MQ=255)
                                                                                                                                                   ggAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTg    >  1:764331/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATA‑TATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAA  >  AP012306/406528‑406826

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