Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP012306406,6770AG23.6% 486.4 / 28.7 195K294E (AAA→GAA) hemHferrochelatase
Reads supporting (aligned to +/- strand):  ref base A (37/112);  new base G (45/1);  total (82/113)
Fisher's exact test for biased strand distribution p-value = 2.66e-20
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 212 total aligned reads displayed.

GCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAA  >  AP012306/406529‑406826
                                                                                                                                                    |                                                                                                                                                      
   ttATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGaaaaa                                                                                                                                                    <  1:4067004/150‑1 (MQ=255)
    tataCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGaaaaaa                                                                                                                                                   >  2:1389868/1‑150 (MQ=255)
    tataCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGaaaaaa                                                                                                                                                   <  1:1580190/150‑1 (MQ=255)
     ataCCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAat                                                                                                                                                  <  1:3078528/150‑1 (MQ=255)
        ccGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAatatt                                                                                                                                               >  1:878274/1‑149 (MQ=255)
         cgacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGa                                                                                                                                              <  2:2414872/150‑1 (MQ=255)
          gacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGaa                                                                                                                                             <  1:353261/150‑1 (MQ=255)
          gacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGaa                                                                                                                                             <  1:1083911/150‑1 (MQ=255)
          gacgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGaa                                                                                                                                             >  1:2752613/1‑150 (MQ=255)
           acgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGaaa                                                                                                                                            >  2:1332942/1‑149 (MQ=255)
           acgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAat                                                                                                                                            >  1:3725502/1‑150 (MQ=255)
            cgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATTaaaa                                                                                                                                           >  2:51105/1‑148 (MQ=255)
            cgaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAata                                                                                                                                           <  1:2183226/150‑1 (MQ=255)
             gaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAatat                                                                                                                                          >  1:563558/1‑150 (MQ=255)
             gaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAatat                                                                                                                                          <  2:3969816/150‑1 (MQ=255)
             gaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAatat                                                                                                                                          <  1:2278209/150‑1 (MQ=255)
             gaAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAatat                                                                                                                                          <  1:3131881/150‑1 (MQ=255)
              aaaCGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAAAGAatata                                                                                                                                         >  2:1927666/1‑150 (MQ=255)
               aaCGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAatatat                                                                                                                                        <  1:1720219/150‑1 (MQ=255)
                aCGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATAtt                                                                                                                                       <  2:426128/150‑1 (MQ=255)
                  gCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGGAAAATATGAATATAtttc                                                                                                                                     >  1:1960764/1‑148 (MQ=255)
                  gCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAAATAtttc                                                                                                                                     >  2:2134642/1‑148 (MQ=255)
                  gCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTcc                                                                                                                                     <  1:2916102/150‑1 (MQ=255)
                  gCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTcc                                                                                                                                     <  1:3663339/150‑1 (MQ=255)
                   cTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCg                                                                                                                                    <  1:1279232/150‑1 (MQ=255)
                    tGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCgg                                                                                                                                   <  2:2039957/150‑1 (MQ=255)
                       aaaTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGgcg‑c                                                                                                                               <  2:2618240/150‑1 (MQ=255)
                        aaTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATATATTCCGgcg‑cc                                                                                                                              >  2:1564716/1‑149 (MQ=255)
                        aaTGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑Ct                                                                                                                              <  2:1604125/150‑1 (MQ=255)
                           gCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTTCGGCG‑CTTaa                                                                                                                           >  1:1608758/1‑150 (MQ=255)
                            cTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAAt                                                                                                                          <  1:51499/150‑1 (MQ=255)
                            cTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAAt                                                                                                                          <  1:4140929/150‑1 (MQ=255)
                              cGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAATATATTCCTTCGCCTTAATg                                                                                                                         >  2:3262835/1‑150 (MQ=255)
                                gagaAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATTAATATATTCCGGCG‑CTTAATGCCa                                                                                                                      >  1:233792/1‑150 (MQ=255)
                                gagaAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCa                                                                                                                      <  1:1397738/150‑1 (MQ=255)
                                gagaAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCa                                                                                                                      <  1:1640414/150‑1 (MQ=255)
                                gagaAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCa                                                                                                                      <  2:2710888/150‑1 (MQ=255)
                                gagaAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAAATAAAATATTCCGGCG‑CTTAAAGCCa                                                                                                                      >  2:1645891/1‑150 (MQ=255)
                                    aaaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATTAATATATTTCCGCG‑CTTTATTGCACGcc                                                                                                                  >  1:3692427/1‑150 (MQ=255)
                                    aaaaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAAATATTTCGGCG‑CGTAATTCCACGcc                                                                                                                  >  2:3161878/1‑150 (MQ=255)
                                      aaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCGGCG‑CTTCATGCCACGCCgg                                                                                                                >  1:185366/1‑150 (MQ=255)
                                      aaGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCgg                                                                                                                <  1:1997231/150‑1 (MQ=255)
                                        ggCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTTCCGCG‑CCTAAAGCCACGCCGGaa                                                                                                              >  1:2313657/1‑150 (MQ=255)
                                        ggCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGaa                                                                                                              <  1:4060751/150‑1 (MQ=255)
                                          cGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACa                                                                                                            <  2:2312872/150‑1 (MQ=255)
                                           gTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATTAATATATTCCCGCG‑CTTAATTGCCCCCCGGAAcca                                                                                                           >  1:3410148/1‑148 (MQ=255)
                                            tAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACata                                                                                                          <  1:3031096/150‑1 (MQ=255)
                                               gTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCg                                                                                                       <  2:3058819/150‑1 (MQ=255)
                                                tCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGGACATATCGa                                                                                                      >  1:774776/1‑150 (MQ=255)
                                                tCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGa                                                                                                      <  1:3369259/150‑1 (MQ=255)
                                                tCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGa                                                                                                      <  1:379040/150‑1 (MQ=255)
                                                   tatTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAat                                                                                                   <  2:815364/150‑1 (MQ=255)
                                                         ggTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATTAATATATTCCGGCG‑CTTTATAGCCCGCCCGGACATATCGAAATGATGGc                                                                                             >  1:1242068/1‑150 (MQ=255)
                                                          gTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATTAATATATTTCGGCG‑CCTTATGCCACGCCGGGACATATCGGAATGAGGGCt                                                                                            >  2:2080796/1‑150 (MQ=255)
                                                             aTGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAAt                                                                                         <  2:2544974/150‑1 (MQ=255)
                                                              tgtgCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATc                                                                                        <  1:577702/150‑1 (MQ=255)
                                                                tgCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCtt                                                                                      <  2:729225/150‑1 (MQ=255)
                                                                 gCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCttg                                                                                     <  2:1249011/150‑1 (MQ=255)
                                                                    cGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCttgttg                                                                                  <  1:1235309/150‑1 (MQ=255)
                                                                        cTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATAAATTCCGGCG‑CTTAATTCCACCCCCGGACATATCGAAATGATGGCTAATCTTGTTGCcgc                                                                              >  2:2500797/1‑150 (MQ=255)
                                                                          ttGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGt                                                                            <  1:3254757/150‑1 (MQ=255)
                                                                             ctgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATTAATATATTTCGGCG‑CTTAATGCCACGCCCGGAAAAATCGAAATGATGGCTAATCTTTTTGTCGCGTATc                                                                         >  2:1431942/1‑150 (MQ=255)
                                                                              tgcGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCg                                                                        <  2:1747667/150‑1 (MQ=255)
                                                                                cGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATATATTCCGGCG‑CTTCATTCCCCCCCCGAACATAACGAAATGATGGCTAATTTTTTTTCCGCGGAACGCt                                                                      >  2:2170078/1‑150 (MQ=255)
                                                                                  gATTGTCTGGAGGCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAAATATTTCGGCG‑CTTAATTCCACGCCGGGACATATCGAAAAGATGGCTAATCTTTTTGCCGCGGATCGCCaa                                                                    >  2:4085932/1‑150 (MQ=255)
                                                                                   aTTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTaaa                                                                   <  2:1264313/150‑1 (MQ=255)
                                                                                      gTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGGTAATCTTTTTGCCGCGTATCGCTAAAACt                                                                >  2:2361022/1‑150 (MQ=255)
                                                                                      gTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAAAGAAAAAATTTCGGCG‑CTTAATTCCCCGCCGGGACATATCGAAATGAAGGGTAATCTTTTTTTCGCGGATTGCTAAAACt                                                                >  2:660676/1‑150 (MQ=255)
                                                                                      gTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCt                                                                <  2:1297384/150‑1 (MQ=255)
                                                                                       tCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTg                                                               <  2:3730268/150‑1 (MQ=255)
                                                                                         tGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAg                                                             <  2:1080839/150‑1 (MQ=255)
                                                                                         tGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAg                                                             <  1:237311/150‑1 (MQ=255)
                                                                                           gagaCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCg                                                           <  1:3288208/150‑1 (MQ=255)
                                                                                             gaCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATTAAAATATTCCGGCG‑CGTAATGCCACGCCGGGAAATATCGAAAAGATGGAGAATCTTGTTGCCGCGTATCGGTAAAACTGAggggg                                                         >  2:3502036/1‑149 (MQ=255)
                                                                                              aCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATAAATATATTTCGGCG‑CTTTATGCCCCGCCGGAAAATAAAGAAATGATAGCTAAACTTTTTGCCGCGTATCGCTTAAGCTGAGAGGTa                                                        >  2:1721134/1‑150 (MQ=255)
                                                                                                   ggAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAAAATATTTCGGCG‑CTTAATGCCACGCCGGGAAATATCGAAATGATGGCTAAACTTGTTGCCGCGGATCGCTAAAGCTGAGCGGGAAAGaa                                                   >  2:3759734/1‑150 (MQ=255)
                                                                                                       gagaTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGa                                               <  1:3993496/150‑1 (MQ=255)
                                                                                                       gagaTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGa                                               <  2:2699269/150‑1 (MQ=255)
                                                                                                         gaTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATATATTTCGGCG‑CTTAATGCCACGCCGGAACATATCGAAAAGATGGGTAATCTTGTTTCCGCGTATCGCTAAAACTGAGCGGGAAAAAAACGAgc                                             >  2:38563/1‑150 (MQ=255)
                                                                                                          aTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATATATTCCGGCG‑CTTAAATCCACGCCGGAACATATTGAAAAGATGGCTAAATTTGTTGGCGCGGTTCGCTAAAGATGAGGGGGAAAGAACTCAggg                                            >  2:2816711/1‑148 (MQ=255)
                                                                                                          aTTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAAAATAATTCGGCG‑CTTAATACCACGCCGGGACATATCCAAATGATGGCTAATCTTGTTGCCGCGGATCCCTAAAGCTTGGCGGTTAAAAACTGAgag                                            >  2:1785570/1‑148 (MQ=255)
                                                                                                            tGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGcc                                          <  1:1424923/150‑1 (MQ=255)
                                                                                                             gCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAATATATTTCGGCG‑CCTTATGCCCCGCCGGAACATATCGAAATGATGGCTAATCTTTTTGCCCCGGATCGCTAAAGCTGAGCGGTAAAAAAATGAGGGccc                                         >  2:1237591/1‑149 (MQ=255)
                                                                                                                         cGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGcc                             <  1:145995/150‑1 (MQ=255)
                                                                                                                           tGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCAAGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtc                           <  1:1159065/150‑1 (MQ=255)
                                                                                                                           tGAGGTCTTCCTCGGTGCCGGCGGGAAAAAAAATTAATATATTTCGGCG‑CTTTATGCCACGCCGGGACATATCGAAATGATGGCTAATCTTTTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAACTGAGCGCCGTCGCGGAGAGCCtc                           >  2:2395851/1‑150 (MQ=255)
                                                                                                                              ggTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgtc                        <  2:2328816/150‑1 (MQ=255)
                                                                                                                               gTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGGACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgtcg                       >  1:3215973/1‑150 (MQ=255)
                                                                                                                               gTCTTCCTCGGTGCCGGCGGGAAAAAAAATGAATATATTTCGGCG‑CTTTATTCCACCCCGGGACATATCGAAATGATGGCTAATCTTTTTGCCGCGTATCGCCAAAACTGAGCGGGAAAAAAATGGGCGCCGTCGCGAAGAGCCtcgtcg                       >  2:1501630/1‑150 (MQ=255)
                                                                                                                                tCTTCCTCGGTGCCGGCGGGGAAAAAAATGAATATATTCCGGCG‑CTTAATTCCCCCCCGGGACAAATCGAAAAGATGGCTAATCTTTTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAACTGAGCGCCGGCGCCAAGAGCCCCGTcgg                      >  2:1589205/1‑150 (MQ=255)
                                                                                                                                tCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTcgg                      <  2:3520457/150‑1 (MQ=255)
                                                                                                                                     cTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGtt                 <  1:1029487/150‑1 (MQ=255)
                                                                                                                                      tCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGttt                <  1:1949330/150‑1 (MQ=255)
                                                                                                                                       cGGTGCCGGCGGGAAAAAAAAAAAAAATATTCCGGCG‑CTTAATGCCACGCCGGGACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAACTGAGCGGGAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGgttt               >  2:1115479/1‑150 (MQ=255)
                                                                                                                                        ggTGCCGGCGGGAAAAAATATGAATATATTTCGGCG‑CTTAATGCCACCCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGGATCGCTAAAGGTGGGCGGGAAAGAACTGAGGGCCGTCGCGAAAAGCCTCGTCGGCGGgtttt              >  2:3391169/1‑149 (MQ=255)
                                                                                                                                             cggcggGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGGACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGGAAAGAACTGAGCGCCGTCGCGAAAAGCCTCGTCGGCGGGTTTtatcat         >  2:132042/1‑150 (MQ=255)
                                                                                                                                               gcggGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATcc       <  1:840812/150‑1 (MQ=255)
                                                                                                                                                cggGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCg      >  1:2242049/1‑150 (MQ=255)
                                                                                                                                                 gggAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGt     >  2:1862578/1‑150 (MQ=255)
                                                                                                                                                  ggAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTg    >  2:2854374/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCG‑CTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAA  >  AP012306/406529‑406826

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