Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,9170AC46.2% 92.7 / 18.2 65V43G (GTG→GGG) ssuBaliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (31/4);  new base C (0/30);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 2.07e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTT  >  NC_000913/993768‑994058
                                                                                                                                                           |                                                                                                                                             
gTAGCAGACGCGCATCCTGGAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCccac                                                                                                                                              >  1:789752/1‑151 (MQ=255)
 tAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCccacc                                                                                                                                             >  2:1072015/1‑151 (MQ=255)
  aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCccacca                                                                                                                                            <  1:162730/151‑1 (MQ=255)
           gcATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCcaca                                                                                                                                    >  2:770971/1‑150 (MQ=255)
            cATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACaa                                                                                                                                   >  2:139623/1‑150 (MQ=255)
              tCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACt                                                                                                                                >  1:237866/1‑151 (MQ=255)
                      catcttCCGCGTTTTTTCTTGATTTCCAC‑CGGGCGGTGTGGTCCCCGCTAACCCCTCC‑CCGGGGTTGGGCTTTC‑CCACACCTGCC‑AGCGGGCGCAGCAGG‑TTCTTTTCG‑CCCCCCCCGCTGGGCCCCACCACCGCCACAAACTGACCTg                                                                                                                          <  2:26134/149‑1 (MQ=255)
                                 tttCTTCCTAATTTTCAG‑CCAGGGGTGTGGTCCCCGTAAACTCCGCC‑CCTGGTTTGGGCTTTC‑CCAGCCCCGCC‑AGCAGGCGCAGCAGG‑TTTTTTTGC‑CCCCCCCCCCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATg                                                                                                                <  2:435109/146‑1 (MQ=255)
                                 tgtCTTCC‑GAATTTCCGCCCAGGGGTGTGGCCCCCCATACCCCATCC‑CCGGGGGTTGGC‑TTT‑CCCGACCCGCCAAGCGGGCGCAGCAGG‑GTCTTTTTA‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                               <  2:468891/149‑1 (MQ=255)
                                      tcttaattttCAG‑CCAGGGGTGTGGTGCCCGCTAACACATCG‑CCTGGGGTTGGCGTTC‑CCAGCCCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTC‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                           <  2:82900/141‑1 (MQ=255)
                                        cTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                      >  2:106517/1‑151 (MQ=255)
                                         tGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                     >  2:484641/1‑151 (MQ=255)
                                          gAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                    >  2:649572/1‑151 (MQ=255)
                                           gaatcccccccagcGGTGGGGCCCCCCCAAACACCGCC‑CCTGGGGTGGGCTTTC‑CCAGCCCCGCC‑AGCGGCCGCAGAGGG‑TTATTTTCC‑CCACCCCCGCTGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                      <  1:384533/140‑1 (MQ=255)
                                           attttCAG‑CCAGGGGTGTGGTCCCCGCAAACACACCG‑CCTGGGTTTGGCTTTC‑CCAGCCCGCCC‑AGCGGGCGCAGCGGG‑GTATTTTTC‑CCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                    <  2:1027950/148‑1 (MQ=255)
                                           aaTTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTgg                                                                                                   >  1:137410/1‑151 (MQ=255)
                                               tCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                >  1:639552/1‑150 (MQ=255)
                                                    ccgggggTGGGGTCCCCCCAACCACCCCG‑CCGGCGGTGGCCTTTC‑CCAGCCCTGCA‑AGCGGGCGCGGAGGG‑GTCTTTTGG‑CCCCCCCCGCGCCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGa                                                                                            <  1:370412/145‑1 (MQ=255)
                                                     cgggggTGTGGTCCCCACTACCCCACCG‑CCGGCGTTGGGCTTTT‑CCAGACCTCCC‑GGCAGGCGCGGCGGG‑GTATTTTTC‑CCCCCCCCGCCGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                    <  2:401629/137‑1 (MQ=255)
                                                     cagcGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                          >  2:29282/1‑151 (MQ=255)
                                                         ggggtggggCCCCCCCACACACCGCCCTGGGGGTGGCTTTTCCCACCCCCGCC‑AGCAGGCCCA‑CAGGTGTACTTTCC‑CCGCCCCCGCCCCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                             <  1:198435/143‑1 (MQ=255)
                                                         ggTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                      >  2:752746/1‑151 (MQ=255)
                                                         ggTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                      >  1:154613/1‑151 (MQ=255)
                                                            gtggtgCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAtttt                                                                                   >  2:268535/1‑151 (MQ=255)
                                                              ggtgCCCGCTAACCCCCCG‑CCTGGGTTTGGCTTTC‑CCAGCCCCCCC‑AGCGGCCGCAGCGGG‑TTACTTTCC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                    <  2:359472/148‑1 (MQ=255)
                                                                  cccGCTAACACATCG‑CGTGGGTTGGGCTTTC‑CCAGCCCTGCC‑AGCAGGCGCAGCAGG‑GTTCTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgcg                                                                              <  2:249986/150‑1 (MQ=255)
                                                                  cccGATAACCCACCG‑CCGGGGGTGGGCTTTC‑CCAGCCCTGCC‑AGCGGGCGCAGCGGG‑TTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCcgc                                                                               <  2:740326/149‑1 (MQ=255)
                                                                     gCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAAt                                                                          >  2:550134/1‑151 (MQ=255)
                                                                        aaCACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                       >  1:24072/1‑151 (MQ=255)
                                                                        aaCACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                       >  1:892392/1‑151 (MQ=255)
                                                                         accccccccccggCGTTTGGCTTTC‑CCAGCCCTCCC‑AGCGGGCGCAGCGGG‑TTATTTTCC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGt                                                                        <  1:47481/138‑1 (MQ=255)
                                                                         acacATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                      >  1:394783/1‑151 (MQ=255)
                                                                             aTCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                  >  1:1056043/1‑151 (MQ=255)
                                                                                   cTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGc                                                             >  2:448926/1‑151 (MQ=255)
                                                                                     gCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                           >  2:671167/1‑151 (MQ=255)
                                                                                       gttgggCTTTC‑CCAGCCCGGCC‑AGCAGGCGCAGCAGG‑GTACTTTTC‑CCACCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                  <  2:623119/138‑1 (MQ=255)
                                                                                         tgggCTTTC‑CCAGCCCTGCC‑AGCGGGCGCGGAGGG‑GTTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAtt                                                          <  2:404587/146‑1 (MQ=255)
                                                                                           ggCTTTC‑CCAGCCCTGCC‑AGCGGGCGCAGCAGG‑GTTCTTTTC‑CCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGca                                                     <  2:389433/151‑1 (MQ=255)
                                                                                                  cccggccccgCC‑AGGGGGCGCAGCGGG‑TTTCTTTCC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                <  1:576271/141‑1 (MQ=255)
                                                                                                   ccagccccccc‑AGCGGGCGCAGAAGG‑GTTTTTTTC‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGCTTGCTTACTGCATTGAGCAACAAt                                                <  1:619433/140‑1 (MQ=255)
                                                                                                   ccAGCCCGGCC‑AGCGGGGGCAGCGGG‑TTATTTGCG‑CCCCCCCGGCGGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGAAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAAt                                                <  1:632896/149‑1 (MQ=255)
                                                                                                     aGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCg                                            >  1:34891/1‑151 (MQ=255)
                                                                                                        ccgccc‑AGCGGGCGCGGCAGG‑TTCCTTTGC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                          <  2:140681/146‑1 (MQ=255)
                                                                                                         ccgccccgcggccgcAGCGGG‑GTCCTTTTC‑CCCCCCCCCCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGc                                         <  2:175596/139‑1 (MQ=255)
                                                                                                            cc‑AGCGGGCGAAGCAGG‑TTCCTTTCC‑CCCCCCCCGCGGCCCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                      <  1:612804/150‑1 (MQ=255)
                                                                                                             ccggcgggCGCAGCGGG‑GTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCt                                      <  1:253042/144‑1 (MQ=255)
                                                                                                               agcagGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                   >  1:804279/1‑151 (MQ=255)
                                                                                                                gcagGCGCAGCGGG‑GTTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                   <  2:249009/150‑1 (MQ=255)
                                                                                                                   ggcgcggaggg‑gTATTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGt                                   <  1:845333/138‑1 (MQ=255)
                                                                                                                   ggCGCAGCAGG‑TTTCTTTTC‑CCCCCCCCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                <  2:1057044/150‑1 (MQ=255)
                                                                                                                       cagcagG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            >  2:852072/1‑150 (MQ=255)
                                                                                                                       cagcagG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            >  1:14652/1‑150 (MQ=255)
                                                                                                                        agcagG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            >  2:753749/1‑149 (MQ=255)
                                                                                                                               gttattttGCCCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTa                     <  1:1091049/147‑1 (MQ=255)
                                                                                                                                     ttccccccccccgccgCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                <  1:293568/140‑1 (MQ=255)
                                                                                                                                     ttcccccccaccGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                   <  1:1083723/140‑1 (MQ=255)
                                                                                                                                      tgcccacccCCGCTCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                  <  1:810143/144‑1 (MQ=255)
                                                                                                                                        cccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg            <  2:645248/144‑1 (MQ=255)
                                                                                                                                         ccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc             <  2:1046397/143‑1 (MQ=255)
                                                                                                                                           accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg         >  1:73289/1‑151 (MQ=255)
                                                                                                                                           accaccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTgg         >  2:507115/1‑151 (MQ=255)
                                                                                                                                             caccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc       >  1:985356/1‑151 (MQ=255)
                                                                                                                                               ccgccgCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcct      <  1:211029/150‑1 (MQ=255)
                                                                                                                                               ccGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctc     <  1:385340/151‑1 (MQ=255)
                                                                                                                                                cGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcc    >  1:97424/1‑151 (MQ=255)
                                                                                                                                                 gCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct   >  2:1048284/1‑151 (MQ=255)
                                                                                                                                                  cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt  >  1:468583/1‑151 (MQ=255)
                                                                                                                                                  cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt  >  1:762031/1‑151 (MQ=255)
                                                                                                                                                           |                                                                                                                                             
GTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAG‑CCAGCGGTGTGGTGCCCGCTAACACATCG‑CCTGCGGTTGGCGTTT‑CCAGACCTGCC‑AGCAGGCGCAGCAGG‑GTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTT  >  NC_000913/993768‑994058

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