Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	993,917	0	A	C	41.3%	91.9 / 16.5	63	V43G (GTG→GGG)	ssuB	aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base A (28/9);  new base C (0/26);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 2.32e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993769‑994063                                                                                                                                                        |                                                                                                                                                   tAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCccacc                                                                                                                                                  >  1:412640/1‑151 (MQ=255)  aGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCccacc                                                                                                                                                  >  2:576288/1‑150 (MQ=255)   gCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCcaccac                                                                                                                                                >  2:9849/1‑151 (MQ=255)                 tGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGa                                                                                                                                   >  1:465619/1‑150 (MQ=255)                  gAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGAcc                                                                                                                                 >  1:797071/1‑151 (MQ=255)                       atcatTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCg                                                                                                                            >  1:339028/1‑151 (MQ=255)                              gcgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATg                                                                                                                     >  2:1033411/1‑151 (MQ=255)                               cgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                    >  1:975869/1‑151 (MQ=255)                               cgTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGt                                                                                                                    >  2:125659/1‑151 (MQ=255)                                  gtCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTaa                                                                                                                  >  1:503465/1‑150 (MQ=255)                                      tCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATcc                                                                                                              >  2:535042/1‑150 (MQ=255)                                       ccTGAATTTCAGCCAGGGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCCCCC‑CCGCTGCGGCCCCCCACCGCCACAGACTGACCTGCCGGAATATGTAAATCCAg                                                                                                            <  2:799318/151‑1 (MQ=255)                                         tGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                          >  1:637385/1‑151 (MQ=255)                                           aaTTCCACCCGGGGGTGTGGTGCCCGCTAACACAGCGCTTGGGGTTGGCT‑TTTCCACCCCAGCAAGGGGGCCCACCGGGGTATTTTCG‑CCACCC‑CCGCGGCGCCCCCCCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGtt                                                                                                          <  2:302353/149‑1 (MQ=255)                                            aTTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGt                                                                                                       >  1:1130113/1‑151 (MQ=255)                                             tttCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGtt                                                                                                      >  2:579556/1‑151 (MQ=255)                                             ttccagccggcGGTGTGGTCCCCGCTAACACACCGCCGGGGGTTGGCT‑TTCCCAGCCCTGCCAGCAGGCGCAGC‑GGGTTCTTTTGCCCCCCC‑CCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGt                                                                                                           <  2:428115/143‑1 (MQ=255)                                                 agccagcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAgg                                                                                                  >  1:1150554/1‑151 (MQ=255)                                                     agcGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACg                                                                                               >  2:593422/1‑150 (MQ=255)                                                        ggTGTGGTCCCCGTTAACACATCCCCTGCGTTGGGCT‑TTCCCAGCCCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCCCCC‑CCGCGGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGatat                                                                                           <  2:92907/151‑1 (MQ=255)                                                          ggtggcccccccaaaCACCCCGCGTGCGGTTGGCT‑TTCCCAGCCCTGCCGGCGGCCGCAGAAGGGTATTTTGG‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                         <  1:393914/138‑1 (MQ=255)                                                           gtggtgCCCGCAAACACACCCCCTGTGTTTGGCT‑TTCCCACCCCCCCCAGCGGGCGCAGCAGGTTTCTTTCC‑CCCCCC‑CCGCTGCCCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGat                                                                                             <  1:271890/146‑1 (MQ=255)                                                             ggtcCCCGCTAACCCACCGCCTGGGGTTGGCT‑TTTCCAGCCCTGCCAGCGGGCGCAGCGGGGTATTTTCC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATAttt                                                                                         <  2:869372/148‑1 (MQ=255)                                                             gggccccGCAAACCCACCGCTGGCGGTGGGCT‑TTCCCAGCCCGCCCCGCGGGCGGAGGGGGTTTTTTTTC‑CCCCCC‑CCGCGGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTcc                                                                                      <  1:939940/147‑1 (MQ=255)                                                                    gCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAAt                                                                               >  2:394748/1‑151 (MQ=255)                                                                       accccaccGCCTGGGTTTGGCT‑TTCCCAGCCCTGCCAGCAGGCGCAGCGGGGTCCTTTTC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                            <  2:282288/144‑1 (MQ=255)                                                                       aaCACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGtt                                                                            >  2:1148027/1‑151 (MQ=255)                                                                         caaatccccggggtttGGCT‑TTCCCAGCCCTCCCACCGGGCGCAGCGGGTTTTTTTTC‑CCCCCC‑CCGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGttt                                                                           <  2:872484/136‑1 (MQ=255)                                                                          cctccGCCTGGGTTGGGCT‑TTCCCACCCCTGCCAGCGGGCGCAGCGGGGTCTTTTTC‑CCCCCA‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGc                                                                         <  2:549917/147‑1 (MQ=255)                                                                          cacaccgcTGGGGTTTGCGTTTTCCACCCCGCCCAGCGGGCGCGGCGGGGTTTTTTTC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTg                                                                          <  1:1007683/144‑1 (MQ=255)                                                                           caTCGCCTGCGGGTGGCG‑TTTGCAGACCTGCCAGCAGGCGCAGCAGGGCACTTTTG‑CCACCA‑CCGCTTCGGCCCACCACCGTCACAAACTGACCTCCCGGACTCTGTAAATCCAGTTGGCTCAGGACGATATTTTCTGCGTAATGttt                                                                           >  1:1143523/1‑148 (MQ=255)                                                                             tCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                       >  1:700979/1‑150 (MQ=255)                                                                                   gCGGTGGGCT‑TTCCCAGCCCTGCCAGCGGGCGCAGCAGGTTTCTTTGC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAt                                                                <  2:1061813/151‑1 (MQ=255)                                                                                    cGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCAtt                                                               >  1:347505/1‑151 (MQ=255)                                                                                       tgggCT‑TTCCCAGCCCTCCCAGCAGGCGCAGCAGGGTTCTTTGC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCa                                                                 <  1:1032959/144‑1 (MQ=255)                                                                                        tgggctttcccaccccTCCCAGCGGGCGCAGGGGGTTACTTTTC‑CCCCCC‑CCGCTGCGCCCCCCCTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGa                                                             <  1:973802/136‑1 (MQ=255)                                                                                        tGGCT‑TTTCCAGCCCTGCCAGCAGCCGCAGCGGGGTCTTTTTG‑CCACCC‑CCGCTGCGGCCCCCCCCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg                                                            <  2:685441/150‑1 (MQ=255)                                                                                          ggtttttCCAGCCCGGCCAGCGGCCGAAGCGGGTTACTTTCC‑CCCCCC‑CCCGTGGGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAg                                                            <  2:966618/145‑1 (MQ=255)                                                                                           cG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGcaac                                                        >  1:528780/1‑151 (MQ=255)                                                                                              ttcccgccccTGCCGGCGGGCGCGGAAGGTTCTTTTTC‑CCACCC‑CGGCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCtt                                                                       <  1:390247/126‑1 (MQ=255)                                                                                                tCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATg                                                    >  1:429929/1‑151 (MQ=255)                                                                                                 ccagcccgcccAGCGGCCGCAGCAGGGTCTTTTCG‑CCCCCC‑CCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATgg                                                   <  2:492491/142‑1 (MQ=255)                                                                                                     aCCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                               >  1:871985/1‑151 (MQ=255)                                                                                                      ccgcccAGCGGGCGAAGCGGGGTATTTTCC‑CCCCCC‑CCTCTGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTg                                               <  2:742374/146‑1 (MQ=255)                                                                                                       cTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGcc                                             >  2:268051/1‑151 (MQ=255)                                                                                                                ggCGCAGCGGGTTACTTTTG‑CCACCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                      <  2:241935/149‑1 (MQ=255)                                                                                                                    cggagtgttttttttC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTc                                      <  1:464725/134‑1 (MQ=255)                                                                                                                    cagcagtttttttttG‑CCCCCC‑CCGCCGCGCCCCACCGCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGAc                                  <  1:636496/138‑1 (MQ=255)                                                                                                                         gggtttttttC‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAgaa                             <  2:105207/143‑3 (MQ=255)                                                                                                                         gggTTCTTTTG‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCa                             <  2:118681/149‑1 (MQ=255)                                                                                                                         gggTTCTTTTC‑CCCCCCCCCGCGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACg                                 <  2:772851/146‑1 (MQ=255)                                                                                                                          ggTACTTTTG‑CCCCCC‑CCGCTGCGGCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTa                          <  2:548707/151‑1 (MQ=255)                                                                                                                          ggTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGt                           >  2:679765/1‑150 (MQ=255)                                                                                                                           gTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAt                         >  2:992414/1‑151 (MQ=255)                                                                                                                               ttttG‑CCCCCC‑CCGCTGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTAtt                        <  2:953408/148‑1 (MQ=255)                                                                                                                               ttttC‑CCCCCC‑CCGTCGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                       <  1:647667/149‑1 (MQ=255)                                                                                                                                tttgcccaccc‑CCGGGGCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                       <  1:1115882/144‑1 (MQ=255)                                                                                                                                tttC‑CCCCCC‑CCGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCAt                     <  2:359744/150‑1 (MQ=255)                                                                                                                                 tttccccccccccGCTCCGCCCCCCCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCa                                     <  1:364205/124‑1 (MQ=255)                                                                                                                                      cccccgcgggggcgcCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTc                       <  1:1042043/133‑1 (MQ=255)                                                                                                                                      ccccccccGCGGCGCCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATAcc                  <  2:919417/143‑1 (MQ=255)                                                                                                                                      cccccacctcggcccccccccTCCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGt                <  1:983471/136‑1 (MQ=255)                                                                                                                                         ca‑ccGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcc            >  2:353604/1‑151 (MQ=255)                                                                                                                                             cGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGcctcc         >  1:640528/1‑151 (MQ=255)                                                                                                                                               cTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCtt       >  1:1152738/1‑151 (MQ=255)                                                                                                                                                    gCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAt  >  1:520348/1‑151 (MQ=255)                                                                                                                                                     cccACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCa     <  1:47473/147‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                   TAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCG‑TTTCCAGACCTGCCAGCAGGCGCAGCAGGGTACTTTTG‑CCACCA‑CCGCTGCGGCCCACCACCGCCACAAACTGACCTGCCGGAATATGTAAATCCAGTTGGTTCAGGACGATATTTTCCGCGTAATGTTTGCTTACTGCATTGAGCAACAATGGCGTGCCCTGGTTCAGACGAGCAGTATTCATACCGTGGCCTCCTTCAAAT  >  NC_000913/993769‑994063

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