Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_bme_pgi639,7240TG27.0% 168.5 / 79.3 100intergenic (+216/‑29)ahpC/ahpFalkyl hydroperoxide reductase, C22 subunit/alkyl hydroperoxide reductase, F52a subunit, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base T (39/34);  new base G (15/12);  total (54/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
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 105 total aligned reads displayed.

GTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAGGAACTGT  >  NC_000913_3_bme_pgi/639574‑639870
                                                                                                                                                      |                                                                                                                                                  
gtAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCat                                                                                                                                                    >  2:370329/1‑151 (MQ=255)
   aGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCatgat                                                                                                                                                 >  1:29880/1‑151 (MQ=255)
       cccGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCatgatgatg                                                                                                                                             <  1:889152/151‑1 (MQ=255)
        ccGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGt                                                                                                                                            <  2:866860/151‑1 (MQ=255)
          ggggATTCACAGGCTATGCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTCCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGtt                                                                                                                                           <  4:558535/150‑1 (MQ=255)
          ggggATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGttt                                                                                                                                          <  2:852704/151‑1 (MQ=255)
           gggATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTa                                                                                                                                         >  2:376676/1‑151 (MQ=255)
           gggATCCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTACCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGATTa                                                                                                                                         <  2:505889/151‑1 (MQ=255)
             gATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAACTTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTaaa                                                                                                                                       >  1:660526/1‑151 (MQ=255)
              aTTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTaaa                                                                                                                                       <  3:550292/150‑1 (MQ=255)
              aTTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAg                                                                                                                                      <  1:214547/151‑1 (MQ=255)
               ttCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGc                                                                                                                                     <  1:298436/151‑1 (MQ=255)
               gagcaaGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGc                                                                                                                                     >  1:726974/6‑151 (MQ=255)
                 cacaGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGcc                                                                                                                                    <  3:84571/150‑1 (MQ=255)
                   caGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAg                                                                                                                                 <  2:190178/151‑1 (MQ=255)
                     ggCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGga                                                                                                                               <  2:548314/151‑1 (MQ=255)
                       cTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGgag                                                                                                                              >  4:495759/1‑150 (MQ=255)
                           ccgccTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATaa                                                                                                                          >  4:132271/1‑150 (MQ=255)
                           ccgccTTGATCTGACGCGAAATCCTTCAGAAATTCCACTAATTCTTAGGGAGCTGCGGCACCCGATTTTCTTCCCCGCACCATGGTGCAAGCGGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATaaa                                                                                                                         >  2:806943/1‑151 (MQ=255)
                            cgccTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCaggat                                                                                                                                                 >  3:485373/1‑126 (MQ=255)
                            cgccTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCaggat                                                                                                                                                 <  4:485373/126‑1 (MQ=255)
                            cgccTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAAc                                                                                                                        <  1:182455/151‑1 (MQ=255)
                            cgccTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAAc                                                                                                                        >  2:658653/1‑151 (MQ=255)
                             gccTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAAc                                                                                                                        >  3:164657/1‑150 (MQ=255)
                                 tGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCCCCCGCTTCTATTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGc                                                                                                                    <  3:117711/150‑1 (MQ=255)
                                   ctctGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCg                                                                                                                 >  2:397728/1‑151 (MQ=255)
                                     ctGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGac                                                                                                               >  2:27689/1‑151 (MQ=255)
                                            gAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAat                                                                                                         <  3:35530/150‑1 (MQ=255)
                                             aaaTACTTCGGAAATTCCCCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAatat                                                                                                       <  2:760667/151‑1 (MQ=255)
                                                 aCTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGaaa                                                                                                   <  1:192409/151‑1 (MQ=255)
                                                   ttCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGaaaa                                                                                                  >  4:572853/1‑150 (MQ=255)
                                                   ttCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaac                                                                                                 >  1:460254/1‑151 (MQ=255)
                                                        aaaTTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaactcaa                                                                                             <  4:498222/150‑1 (MQ=255)
                                                              aCCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAgg                                                                                      <  2:751349/151‑1 (MQ=255)
                                                                cTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGc                                                                                     >  3:530196/1‑150 (MQ=255)
                                                                 aagaaaatcGGGTGCCGCAGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCt                                                                                    >  3:396765/8‑150 (MQ=255)
                                                                    aaaatcGGGTGCCGCAGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTAcc                                                                                >  1:481194/5‑151 (MQ=255)
                                                                     aaatcGGGTGCCGCAGCACCCGATTTTCTTCCCCGAACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCt                                                                               >  2:818078/4‑151 (MQ=255)
                                                                       ttcGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCtt                                                                              <  3:399084/150‑1 (MQ=255)
                                                                         cGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTga                                                                            <  3:57749/150‑1 (MQ=255)
                                                                         cGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTga                                                                            <  4:211970/150‑1 (MQ=255)
                                                                         cGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTga                                                                            >  4:126168/1‑150 (MQ=255)
                                                                          gggTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTgaga                                                                          >  1:175317/1‑151 (MQ=255)
                                                                          gggTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTg                                                                             <  4:134174/148‑1 (MQ=255)
                                                                          gggTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTg                                                                             >  3:134174/1‑148 (MQ=255)
                                                                            gTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGaa                                                                         >  3:120336/1‑150 (MQ=255)
                                                                                 gcGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAATGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGAc                                                                   <  1:370329/151‑1 (MQ=255)
                                                                                  cGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGAc                                                                   >  3:461739/1‑150 (MQ=255)
                                                                                  cGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGAc                                                                   >  3:261428/1‑150 (MQ=255)
                                                                                    gCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCaa                                                                >  2:133821/1‑151 (MQ=255)
                                                                                      aCCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGc                                                              >  1:152589/1‑151 (MQ=255)
                                                                                        ccGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaac                                                                                                 <  4:62159/114‑1 (MQ=255)
                                                                                        ccGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAaac                                                                                                 >  3:62159/1‑114 (MQ=255)
                                                                                         cGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCt                                                            >  4:240665/1‑150 (MQ=255)
                                                                                         cGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTg                                                           >  2:397771/1‑151 (MQ=255)
                                                                                         cGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGgaga                                                                                                                             >  3:470674/1‑85 (MQ=255)
                                                                                         cGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGgaga                                                                                                                             <  4:470674/85‑1 (MQ=255)
                                                                                            ttttCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGtt                                                         <  3:132271/150‑1 (MQ=255)
                                                                                             tttCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTCTTg                                                        <  4:461739/150‑1 (MQ=255)
                                                                                                cttcCCCGCACCATGATGCAAGCTGCAACCAGGTAGCCGCAGAGGCCTCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAACAATTGACCACGCCTGTTGCGtt                                                    >  2:266063/1‑151 (MQ=255)
                                                                                                 ttcCCCGCACCATTATGCAAGCTGCAGCCAGTTAGCCGCATAGGCCTTTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTa                                                   <  2:82328/151‑1 (MQ=255)
                                                                                                 ttcCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATaaa                                                                                                                         >  4:435688/1‑81 (MQ=255)
                                                                                                 ttcCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATaaa                                                                                                                         <  3:435688/81‑1 (MQ=255)
                                                                                                 ttcCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTa                                                   <  1:376676/151‑1 (MQ=255)
                                                                                                    cccGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAAtt                                                <  1:492810/151‑1 (MQ=255)
                                                                                                    cccGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAAtt                                                >  1:611956/1‑151 (MQ=255)
                                                                                                     ccGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAAtt                                                <  4:12921/150‑1 (MQ=255)
                                                                                                       gCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAg                                                                                       <  4:438621/109‑1 (MQ=255)
                                                                                                       gCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAg                                                                                       >  3:438621/1‑109 (MQ=255)
                                                                                                           cATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGc                                         >  2:284118/1‑151 (MQ=255)
                                                                                                               atgCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTgg                                      >  4:419083/1‑150 (MQ=255)
                                                                                                                tgCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGAt                                    <  1:191712/151‑1 (MQ=255)
                                                                                                                 gCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTg                                                                             <  2:947622/109‑1 (MQ=255)
                                                                                                                 gCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTg                                                                             >  1:947622/1‑109 (MQ=255)
                                                                                                                 gCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATg                                   <  2:660526/151‑1 (MQ=255)
                                                                                                                 gCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATg                                   >  1:608205/1‑151 (MQ=255)
                                                                                                                   aaGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGa                                  >  3:567625/1‑150 (MQ=255)
                                                                                                                      cTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGac                                                                                                               <  1:269901/70‑1 (MQ=255)
                                                                                                                      cTGCATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGac                                                                                                               >  2:269901/1‑70 (MQ=255)
                                                                                                                         cATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCt                           >  1:572461/1‑151 (MQ=255)
                                                                                                                         cATCCAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAGACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCt                           >  2:730422/1‑151 (MQ=255)
                                                                                                                            ccAGGTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTaa                         <  4:396765/150‑1 (MQ=255)
                                                                                                                               ggTAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCg                     >  2:566068/1‑151 (MQ=255)
                                                                                                                                gTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCg                     <  3:353318/150‑1 (MQ=255)
                                                                                                                                 tAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGc                   <  2:481194/151‑1 (MQ=255)
                                                                                                                                 tAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAGACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGc                   <  1:730422/151‑1 (MQ=255)
                                                                                                                                 tAGCCGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGc                   >  1:3417/1‑151 (MQ=255)
                                                                                                                                  aGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCa                  <  2:29880/151‑1 (MQ=255)
                                                                                                                                   gCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCa                  <  4:261428/150‑1 (MQ=255)
                                                                                                                                   gCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAg                 >  2:17258/1‑151 (MQ=255)
                                                                                                                                    ccGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGa                >  1:785764/1‑151 (MQ=255)
                                                                                                                                    ccGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACCTGCTCGAAACAAATATGAAAACTCCACGCGAGCCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGa                <  1:519332/151‑1 (MQ=255)
                                                                                                                                     cGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGa                >  3:142822/1‑150 (MQ=255)
                                                                                                                                     cGCAGAGGCCGCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGaa               >  2:366172/1‑151 (MQ=255)
                                                                                                                                         gagGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCa           <  2:726974/151‑1 (MQ=255)
                                                                                                                                            gCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAg         >  3:111215/1‑150 (MQ=255)
                                                                                                                                            gCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAg         <  4:120336/150‑1 (MQ=255)
                                                                                                                                              cGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAGGaa      >  2:401403/1‑151 (MQ=255)
                                                                                                                                               gCTTGCAGGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAGGAAc     <  1:366172/151‑1 (MQ=255)
                                                                                                                                                  tGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAGGAACTGt  >  2:540943/1‑151 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GTAAGCTCCCGGGGATTCACAGGCTAGCCGCCTTGCTCTGACGCGAAATACTTCGGAAATTCACCTAATTCTTCGGGTGCTGCGGCACCCGATTTTCTTCCCCGCACCATGATGCAAGCTGCATCCAGGTAGCCGCAGAGGCCGCTTGCATGATGATGTTTAAAGCCCAGGAGATAAACATGCTCGACACAAATATGAAAACTCAACTCAAGGCTTACCTTGAGAAATTGACCAAGCCTGTTGAGTTAATTGCCACGCTGGATGACAGCGCTAAATCGGCAGAAATCAAGGAACTGT  >  NC_000913_3_bme_pgi/639574‑639870

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