Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,913	0	G	C	39.7%	137.0 / 20.8	74	G44G (GGC→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (41/3);  new base C (1/28);  total (43/31)
Fisher's exact test for biased strand distribution p-value = 9.95e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTT  >  NC_000913/993770‑994058                                                                                                                                                  |                                                                                                                                                  aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCccacca                                                                                                                                            <  2:1074802/151‑1 (MQ=255) aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCccacca                                                                                                                                            <  2:1139942/151‑1 (MQ=255)   cAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCaccacc                                                                                                                                          >  2:973478/1‑151 (MQ=255)    aGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCg                                                                                                                                         >  1:759794/1‑151 (MQ=255)       cgcgCATCCTGAAACATCATTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGcc                                                                                                                                       >  1:1488769/1‑150 (MQ=255)           cATCCTGAAACATCATTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACaa                                                                                                                                   >  1:328750/1‑150 (MQ=255)                 gacaCATCATTCGCGTCTTTTCCTG‑AATTTCCGCCAGCGGTGTGGTGCCCGCTCACACATTGCCTGCGCTTGGCGCTTCCCAACCTGCCCGCCCGCCCAGCAGCGTACTTTTG‑CCACCACCGCTGCGCTTCACCACCGCCACAAACTGAcc                                                                                                                            >  2:1422874/4‑151 (MQ=255)                 aaaaCATCATTCGGGTTTC‑TCCTGAAATTCCGGCCGGGGGTGTGGGCCCCGCAACCACACCGCGGGGGGTTGGTTTTCCCACCCCCGCCAGCAGGCGCGGGGGGTTACTTTTC‑CCCCCCCCGCGGGGCCCCACCACCGCCACAAACTGAcc                                                                                                                            <  2:543608/150‑1 (MQ=255)                  aaaCATCATTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCt                                                                                                                           >  2:324668/1‑151 (MQ=255)                  aaaCATCATTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCt                                                                                                                           >  2:921544/1‑151 (MQ=255)                  aaaCATCATTCGCACGCCTTCCTG‑AATTACAGCCACCAGTGTGGGGCCCGCTACCACATCGCCTGCGTTTAGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGCACTTTTG‑CGACCACCCATGCGTCCCACCACACCcacgcaccgac                                                                                                                             >  2:318919/1‑141 (MQ=255)                   aacaccttccgcgTGTCTCCCTA‑ATTTCCAGCCAGGGGTGTGGGCCCCGCAAACACACCCCCGGGGGTTGGCTTTCCCAGCCCCGCCAGCGGGCGCAGCGGGGTTCTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTg                                                                                                                          <  1:704553/142‑1 (MQ=255)                    aCATCATTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGc                                                                                                                         >  2:1337118/1‑151 (MQ=255)                      atcatTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                       >  2:691892/1‑151 (MQ=255)                                tgtCTTCCTA‑AATTCCAGCCAGCGGTGTGGTGCCCGCTAACACACCGCCTGGGTTGGGCTTTCCCAGCCCGGCCAGCGGGCGCAGCGGGTTTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATg                                                                                                                <  1:872725/148‑1 (MQ=255)                                  tCTTCCTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAAt                                                                                                           >  1:335261/1‑151 (MQ=255)                                   cTTCCTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATc                                                                                                          >  2:1350854/1‑151 (MQ=255)                                    ttCCTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAAGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGGATATGTAAATcc                                                                                                         >  2:1424771/1‑151 (MQ=255)                                      ccTG‑ATTTTCAGCCGGGGGTGTGGGCCCCGCAAACCCATCGCCGGCGGTGGGCTTTCCCAGCCCTGCCAGCGGGCGCAGCGGGTTTCTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCa                                                                                                        <  1:1321814/150‑1 (MQ=255)                                      ccTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                         >  1:470470/1‑149 (MQ=255)                                         g‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGCACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                    >  2:604537/1‑151 (MQ=255)                                           attttCACCCAGGGGTGGGGTCCCCGCAAACCCACCGCCGGCGGTGGGCGTTCCCAGCCCTGCCAGAAGGCGCGGCGGGTTCTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATg                                                                                                                <  1:843902/136‑1 (MQ=255)                                           atttccacccgggggTGGGGTCCCCGCAAACACCCCCCCTGGGTTTGGCTTTCCCAGCCCTCCCGGCGGGCGCAGC‑GGGTTCTTTTGCCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTg                                                                                                    <  1:1565428/137‑1 (MQ=255)                                              ttCAGCCCGCCGTGTGGTGCCCGCTAACACATCTCCTGCGCTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTC‑CCAGCACCGCTCCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATTTAAATCCAGTTGGtt                                                                                                 >  2:1092021/1‑150 (MQ=255)                                               tCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                >  1:668486/1‑150 (MQ=255)                                                cagccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                >  1:559599/1‑149 (MQ=255)                                                cagccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                >  1:520790/1‑149 (MQ=255)                                                   ccgggggggtggtcCCCCCAAACCCTCCGCCTGGGGTGGGCGTTCCCACCCCTGCCAGCAGGCGCAGCGGGGTATTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                <  1:93316/138‑1 (MQ=255)                                                   cccggggTGTGGCCCCCCATAACACCCCGCCTGGGTTGGGCTTTCCCACCCCCGCCAGCGGGCCCAGCGGGGTCCTTTCC‑CCACCCCCGCTCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTc                                                                                                <  2:1313777/141‑1 (MQ=255)                                                   ccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGAc                                                                                           >  2:775897/1‑151 (MQ=255)                                                    cggggggggggtcCCCCCAAACCCACCGCCGGGGGTTGGCTTTCCCACCCCGCCCGGCGGGCGCAGAGGGGTTTTTTCG‑CCCCCACCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGa                                                                                            <  2:696532/140‑1 (MQ=255)                                                    cagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                          >  1:575039/1‑151 (MQ=255)                                                    cagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                          >  1:1513982/1‑151 (MQ=255)                                                    cagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                          >  1:1135512/1‑151 (MQ=255)                                                     agcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                          >  2:141038/1‑150 (MQ=255)                                                        gggggggtcCCCGCAACAACACCGCCGGGGTTTGGCTTTCCCACCCCTGCCAGCAGCCGCAGCGGGTTTCTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                      <  1:1369953/146‑1 (MQ=255)                                                                    gCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAAt                                                                          >  2:1325996/1‑151 (MQ=255)                                                                       aaCACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAAT‑CAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                      >  1:182590/1‑151 (MQ=255)                                                                           caTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCt                                                                   >  2:1211086/1‑151 (MQ=255)                                                                            attgCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                  >  2:110224/4‑151 (MQ=255)                                                                            aTCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                  >  2:68833/1‑151 (MQ=255)                                                                              cGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTAc                                                                >  1:1094058/1‑151 (MQ=255)                                                                               gtcTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCAAAAACTGACCTGACGGAATATGTAAATCCAGTTGGCTCAGGACGACATTTTCCGCGTAATGTTTGCTTACt                                                               >  2:1397838/3‑151 (MQ=255)                                                                               gCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACt                                                               >  2:68749/1‑151 (MQ=255)                                                                                   gcgttgggCTTTCCCAGCCCGGCCAGCAGGCGCAGCAGGTTTTTTTTG‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                              <  1:115427/142‑1 (MQ=255)                                                                                     gtttGGCTTTCCCAGCCCTCCCAGCGGGCGCAGCGGGTTCTTTTTC‑CCCCCCCCGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTg                                                              <  1:896156/144‑1 (MQ=255)                                                                                     gtttGGCTTTCCCACCCCGGCCAGCGGGCGCGGCGGGGTTTTTTTC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTg                                                         <  1:1376251/149‑1 (MQ=255)                                                                                       tgggCTTTCCCAGCCCTGCCGGCGGGCGCAGCGGGGTTTTTTTC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                    <  1:416275/136‑1 (MQ=255)                                                                                        tGGAGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGc                                                      >  2:242831/1‑151 (MQ=255)                                                                                             ttcccAGCCCTGCCGGCGGGCGGAGCGGGGTACTTTGC‑CCCCCCCCGGTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                   <  1:657755/146‑1 (MQ=255)                                                                                                               ggCGCGGCGGGGTTTTTTGC‑CCCCCCCCGCCGCGCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                <  1:1211850/150‑1 (MQ=255)                                                                                                                gcgcAGCGGGGTACTTTTC‑CCCCCCCCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAg                               <  1:640884/150‑1 (MQ=255)                                                                                                                 cgcAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                             >  1:189277/1‑151 (MQ=255)                                                                                                                   cagcagGTTTCTTCCG‑CCCCCCCCGCGGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                              <  2:987686/148‑1 (MQ=255)                                                                                                                    agcagGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            >  1:343896/1‑149 (MQ=255)                                                                                                                    agcagGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            >  2:388676/1‑149 (MQ=255)                                                                                                                      gggggTATTTTTGCCCCCCCCCGCGCCGCCCCGCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                            <  2:759781/146‑1 (MQ=255)                                                                                                                      cagGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                         >  1:1107473/1‑150 (MQ=255)                                                                                                                       agGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                         >  1:777452/1‑149 (MQ=255)                                                                                                                        ggttttttttC‑CCCCCCCCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGa                              <  1:820682/137‑1 (MQ=255)                                                                                                                        gggtctttttC‑CCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGc                         <  1:793733/142‑1 (MQ=255)                                                                                                                        gggtctttgcg‑CCCCCCCCGCGGCCCCCCACCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                             <  2:1483730/134‑1 (MQ=255)                                                                                                                           tatttttC‑CCCCCCCGGCTGCGCCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTa                     <  2:17225/146‑1 (MQ=255)                                                                                                                           tACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                   >  1:502254/1‑151 (MQ=255)                                                                                                                               tttgcccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                   <  1:489728/137‑1 (MQ=255)                                                                                                                               tttG‑CCACCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                  <  1:111528/148‑1 (MQ=255)                                                                                                                               tgtacccccccccgccgCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCa                 <  1:763864/139‑1 (MQ=255)                                                                                                                                ttgcccccccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                  <  1:1145170/138‑1 (MQ=255)                                                                                                                                  g‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg            >  2:1210199/1‑151 (MQ=255)                                                                                                                                    ccccccccGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCg            <  1:1096202/144‑1 (MQ=255)                                                                                                                                        ccccGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc       <  1:163290/149‑1 (MQ=255)                                                                                                                                        caccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc       >  1:7930/1‑151 (MQ=255)                                                                                                                                            gCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcct   >  2:916311/1‑151 (MQ=255)                                                                                                                                             cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt  >  1:1521411/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                  AGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTG‑AATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCACCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTT  >  NC_000913/993770‑994058

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