Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0123062,326,9790AC21.1% 686.7 / 15.3 273intergenic (‑5/+68)gshA/yqaAgamma‑glutamate‑cysteine ligase/conserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (127/86);  new base C (0/57);  total (127/146)
Fisher's exact test for biased strand distribution p-value = 2.45e-19
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 334 total aligned reads displayed.

CTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACAAAAAAAA  >  AP012306/2326830‑2327127
                                                                                                                                                     |                                                                                                                                                    
 tAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGAc                                                                                                                                                     <  2:1727411/150‑1 (MQ=255)
    tGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTc                                                                                                                                                  <  1:3943207/150‑1 (MQ=255)
     gCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTcc                                                                                                                                                 >  2:159334/1‑150 (MQ=255)
      cTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTccc                                                                                                                                                >  1:2607876/1‑150 (MQ=255)
      cTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTccc                                                                                                                                                <  1:1754442/150‑1 (MQ=255)
       ttCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGACGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCg                                                                                                                                               >  1:4676203/1‑150 (MQ=255)
         cAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGCCCTCCCGcc                                                                                                                                             <  2:621867/150‑1 (MQ=255)
          aGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCt                                                                                                                                            <  2:4231554/150‑1 (MQ=255)
            gATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGt                                                                                                                                          >  1:509499/1‑150 (MQ=255)
              tGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTTTCa                                                                                                                                        <  2:3036877/150‑1 (MQ=255)
                   tgttgtTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGCCCTCCCGCCTGTCAAAATc                                                                                                                                   <  1:3459447/150‑1 (MQ=255)
                   tgttgtTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATc                                                                                                                                   <  1:1654075/150‑1 (MQ=255)
                    gttgttGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCg                                                                                                                                  >  1:4038321/1‑150 (MQ=255)
                            cAGTGTGCCATCTGCTTTAACACGCAAGGTTGCGCGCTCCAGCCCACGCTATACCCCATTAAACGCTGGAGGATTTTTTCCCTGCCAGGCCAGGGCCTGTAATACCTCCGGGATCAAATTGACCCCCCGCCTGTCAAAATCGTTTTAAtt                                                                                                                          <  2:3185610/150‑1 (MQ=255)
                            cAGTGTGCCATCAGCATTAACACGCAAATTTTCGCGCCCCAGCCCCCGCTGTACCCCTTTTAACGCCGGAGGATTTTTTCCCAGCCAGGCCAGCGCTTGTGATACGTCCGGGATCAAATTGCCCTCCCGCCTGTCAAAATCGTTTTAAtt                                                                                                                          <  2:2812643/150‑1 (MQ=255)
                            cAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAAtt                                                                                                                          <  1:884934/150‑1 (MQ=255)
                            aaGTGTGCCATCAGCATCAACCCGCAAAGTTTCGCGCCCCACCCCACGCTGTATCCCTTTTAACGCCTGGGGATTTTTTTCCACCCGGCCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAAtt                                                                                                                          <  2:4327772/149‑1 (MQ=255)
                             aGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTa                                                                                                                         <  1:2663259/150‑1 (MQ=255)
                             aGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTa                                                                                                                         >  1:2425345/1‑150 (MQ=255)
                             aGTGTGCCATCAGCATTAACACGAAAAGTTTCGCGCCCCAGCCCACGCTGTATCCCCTTTACCGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATAAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTa                                                                                                                         <  2:1924617/150‑1 (MQ=255)
                              gtgtGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAg                                                                                                                        <  1:3082414/150‑1 (MQ=255)
                                  gCCATCAGCATTAACCCGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTCCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATAAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATa                                                                                                                    <  2:4501360/150‑1 (MQ=255)
                                  gCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATa                                                                                                                    >  1:2827625/1‑150 (MQ=255)
                                  gCCATCAGCATTAACACGCAAAGTTCCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATa                                                                                                                    <  1:316256/150‑1 (MQ=255)
                                   ccATCAGCATTACCCCGAAAATTTTCCCGCCCCAGCCCCCGCTGTACCCCCTTAAACGCCTGAGAATTTTTTCCCACCCAGGCCGGCGCCTGTGATCCGTCCGGGATCAAATTGCCCTCCCGCCTGTCAAAATCGTTTTAATTAGCATaa                                                                                                                   <  2:1169469/150‑1 (MQ=255)
                                   ccATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGCCAAAATCGTTTTAATTAGCATaa                                                                                                                   <  1:1492681/150‑1 (MQ=255)
                                    cATCAGAATTACCCCGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGCCCCCCCGCCTGTCAAAATCGTTTTAATTAGCATAAc                                                                                                                  <  2:1087703/150‑1 (MQ=255)
                                     aTCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACt                                                                                                                 <  1:1563999/150‑1 (MQ=255)
                                     aTCAGCATTAACACGCAAAGTTCCGCCCCCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATTTTTTCCCAGCCAGGCCAGCGCCTGTTATACGTCCGGGATAAAATTGCCCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACt                                                                                                                 <  2:2230626/150‑1 (MQ=255)
                                         gCATTAACACGAAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTACCGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCTTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTaa                                                                                                             <  2:921995/150‑1 (MQ=255)
                                           aTTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTACCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCGGGAATAAAATTGCCCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATg                                                                                                           <  2:3078747/150‑1 (MQ=255)
                                                cacGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATAAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGAc                                                                                                      <  1:2992663/150‑1 (MQ=255)
                                                  cGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTACCGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGCCCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCa                                                                                                    <  2:411304/150‑1 (MQ=255)
                                                   gCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCGGGAATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCat                                                                                                   <  1:676907/150‑1 (MQ=255)
                                                    gaaaGTTTCGCGCCCCAGCCCACGGTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCata                                                                                                  <  1:1596292/149‑1 (MQ=255)
                                                     aaaGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCatat                                                                                                 >  1:4544505/1‑150 (MQ=255)
                                                      aattttCCCGCCCCAGCCCGGGCTGTATCCCTTTGACCGCCTGGGGATGTTTTTCCAGCCAGCCGAGCCCCGTTAATACGGCAGGGAACAATTGTTCATCCCGCCTGTCAAAATCGTTTTAATTAGCATACCTGTAATGGTGACCATAtg                                                                                                <  2:187937/147‑1 (MQ=255)
                                                       aGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATAtgt                                                                                               >  2:3876275/1‑150 (MQ=255)
                                                         tttCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGc                                                                                             >  2:2525400/1‑150 (MQ=255)
                                                         tttCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGc                                                                                             >  1:4461279/1‑150 (MQ=255)
                                                           tCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAg                                                                                           >  1:2348043/1‑150 (MQ=255)
                                                               gcTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCa                                                                                            <  1:4569661/145‑1 (MQ=255)
                                                               gcTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCt                                                                                       >  1:2249551/1‑150 (MQ=255)
                                                                cccgacccccCGCTGCATCCCTTTAAACGCCTGAGATTGTTTTTCCAGCCAGGCAGGCGCATGTTATACGTCGGGGATCAATTTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTa                                                                                      <  2:2736898/144‑1 (MQ=255)
                                                                cTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTa                                                                                      >  1:1817673/1‑150 (MQ=255)
                                                                cTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTa                                                                                      >  2:2265708/1‑150 (MQ=255)
                                                                 tCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTAc                                                                                     >  2:4592410/1‑150 (MQ=255)
                                                                 cccAGCCCCCGTTGTACCCCTTTAAACCCCGGAGGATTTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCCCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTAc                                                                                     <  2:4454698/149‑1 (MQ=255)
                                                                    aGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAAt                                                                                  >  1:2922322/1‑150 (MQ=255)
                                                                      cccACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGCGCAGGCCTACAATTa                                                                                <  1:3675663/150‑1 (MQ=255)
                                                                       ccACGCTGTACCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCGGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAg                                                                               <  2:3265339/150‑1 (MQ=255)
                                                                        cgcGCTGTCCCCCTTTTAACGCCTGAGGATTTTTTTCCAGCCAGGCCAGCGCTTGTGATACGTCCGGGATCAAATTGCCCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGt                                                                              <  2:834946/148‑1 (MQ=255)
                                                                          cGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGCCCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGc                                                                            <  2:772043/150‑1 (MQ=255)
                                                                          cGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGc                                                                            >  2:128874/1‑150 (MQ=255)
                                                                          cGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGc                                                                            >  2:1530334/1‑150 (MQ=255)
                                                                                aTCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCcaccac                                                                      >  1:2894249/1‑150 (MQ=255)
                                                                                  ccccTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCAcat                                                                    >  2:1441171/1‑150 (MQ=255)
                                                                                     cTTTAACGCCTGAGGATGTTTTCCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGCCCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCAcatcat                                                                 <  2:3326896/150‑1 (MQ=255)
                                                                                         aaCGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGccc                                                             <  1:3259955/150‑1 (MQ=255)
                                                                                          ccGCCAGAGGATTTTTTTCCAGCGAGCCCGGGGCTTGTGATACGTCCGGTAACAATTTGCCCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCt                                                            <  2:1771226/149‑1 (MQ=255)
                                                                                          aCGCTTGAGGTTTTTTTTCCGGCCTGCCGAGCGCCTGTGATACCGGGGGGATAAAGTTGCCCCCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGTCCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCt                                                            <  2:4407736/150‑1 (MQ=255)
                                                                                          aCGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCt                                                            >  1:4866425/1‑150 (MQ=255)
                                                                                             ccTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGaa                                                         >  1:3372276/1‑150 (MQ=255)
                                                                                              cTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGCCCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAAc                                                        <  1:2851468/150‑1 (MQ=255)
                                                                                              cTGAGGATGTTTTCCCACCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGCCCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAAc                                                        <  2:545619/150‑1 (MQ=255)
                                                                                               ggggattGTTTTCCCAGCACGGCCAGCGGTTGTATTATGTCCGGGTCAAAATTGCCCTCCCGCCTGTCAAAATCGTGTTAATTAGCATAACTGTAATGGTGACCATAGGTGCAGGCCTACAATTAGTGCCCCCACATCATGCTCTGAACg                                                       <  2:2018791/144‑1 (MQ=255)
                                                                                                  ggATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTc                                                    >  1:1138712/1‑150 (MQ=255)
                                                                                                     tGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCgct                                                 <  1:1003347/150‑1 (MQ=255)
                                                                                                        ttttCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCgctgct                                              <  1:3839631/150‑1 (MQ=255)
                                                                                                         tttCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGctgctg                                             >  1:1417270/1‑150 (MQ=255)
                                                                                                               gccaggccagCGCCTGTGATACGTCCGGGATCAAATTGCCCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGctgcaactg                                       <  2:2123401/150‑1 (MQ=255)
                                                                                                               accaggccggCGCTTGGGTTGCGTCCGGTATAAAATTCCCCCCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGctgcaactg                                       <  2:958399/149‑1 (MQ=255)
                                                                                                                   ggcaagCGGCTGTGATACGCCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCtgcaactgcaac                                   <  2:3400125/150‑1 (MQ=255)
                                                                                                                      cagCGCCTGTGATACGTCCGGGATCAAATTGCCCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATa                                <  2:4193162/150‑1 (MQ=255)
                                                                                                                       agCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAg                               >  1:2540235/1‑150 (MQ=255)
                                                                                                                          gcCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGt                            <  1:3248178/150‑1 (MQ=255)
                                                                                                                            cTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTaa                          >  1:3738097/1‑150 (MQ=255)
                                                                                                                             tgtgATACGTCCGGGATCAAATTGCCCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAAc                         <  2:2661978/150‑1 (MQ=255)
                                                                                                                              gtgtTACTTCGGGAATAAATTTGCCCCCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACg                        <  2:2420983/150‑1 (MQ=255)
                                                                                                                              gtgATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACg                        <  1:3901336/150‑1 (MQ=255)
                                                                                                                              gtgATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACg                        >  2:2826512/1‑150 (MQ=255)
                                                                                                                               tgATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGc                       >  1:606499/1‑150 (MQ=255)
                                                                                                                                aaTACGTCCGGGATAAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCt                      <  2:4935380/149‑1 (MQ=255)
                                                                                                                                   ccgcccGGAATAAAATTGACTCCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTa                   <  1:476255/146‑1 (MQ=255)
                                                                                                                                   aCGTCGGGGATCAAATTGCCCCCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTa                   <  2:3151908/150‑1 (MQ=255)
                                                                                                                                     gcccGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTAcc                 <  2:2558704/148‑1 (MQ=255)
                                                                                                                                     gTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTAcc                 >  1:4082322/1‑150 (MQ=255)
                                                                                                                                     gTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTAcc                 >  1:1379680/1‑150 (MQ=255)
                                                                                                                                        cGGGATCAAATTGACCCCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAg              <  1:380480/150‑1 (MQ=255)
                                                                                                                                         gggATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAgg             <  1:805288/150‑1 (MQ=255)
                                                                                                                                         gcattCAAATGGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAgg             <  2:3376591/146‑1 (MQ=255)
                                                                                                                                           gATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGca           >  2:3831167/1‑150 (MQ=255)
                                                                                                                                            aTCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGcac          >  2:3554644/1‑150 (MQ=255)
                                                                                                                                             tCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGcaca         <  2:4774754/150‑1 (MQ=255)
                                                                                                                                             tCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGcaca         <  1:2049252/150‑1 (MQ=255)
                                                                                                                                               aaaTTGCCCCCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACaaa       <  2:4645101/150‑1 (MQ=255)
                                                                                                                                                aaTTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACaaaa      >  1:3510200/1‑150 (MQ=255)
                                                                                                                                                 aTTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACaaaaa     >  2:1791505/1‑150 (MQ=255)
                                                                                                                                                  ttGCCCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACaaaaaa    <  1:3579819/150‑1 (MQ=255)
                                                                                                                                                   tGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACaaaaaaa   >  1:1438444/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTAATGCTTCAGGATGACCTGTTGTTGCCAGTGTGCCATCAGCATTAACACGCAAAGTTTCGCGCTCCAGCCCACGCTGTATCCCCTTTAACGCCTGAGGATGTTTTTCCAGCCAGGCCAGCGCCTGTGATACGTCCGGGATCAAATTGACCTCCCGCCTGTCAAAATCGTTTTAATTAGCATAACTGTAATGGTGACCATATGTGCAGGCCTACAATTAGTGCCACCACATCATGCCCTGAACGGTCGCTGCTGCAACTGCAACATAGCGTAACGCTTTACCAAGGCACAAAAAAAA  >  AP012306/2326830‑2327127

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