Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,233	0	G	A	48.3%	46.1 / 33.0	59	pseudogene (60/75 nt)	ycgH	pseudogene; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (22/8);  new base A (0/28);  total (22/37)
Fisher's exact test for biased strand distribution p-value = 1.11e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222083‑1222382                                                                                                                                                        |                                                                                                                                                       gtgtGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTg                                                                                                                                                        >  2:532163/1‑151 (MQ=255)      aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAata                                                                                                                                                   <  1:1209068/151‑1 (MQ=255)                       ggCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCccgcc                                                                                                                                 >  1:400525/1‑151 (MQ=255)                         cGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGt                                                                                                                               >  1:1213855/1‑151 (MQ=255)                           aGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGtt                                                                                                                              >  2:261303/1‑150 (MQ=255)                               aTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCt                                                                                                                          >  2:764597/1‑150 (MQ=255)                                tttGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTa                                                                                                                         >  1:284478/1‑150 (MQ=255)                                ttgggggaAAATTAAACTTGGCGAAAAGGGCACCCGA‑AACCCAGCGGCCAGCGTCGGGGGAATGTGGGGGTGCGGCTGGGTGATAATGGCTACAATGCCCCCGCAGTGATGGGGGGCCTAAAATATAA‑TTTCTGATCCCGCCGTTAGCTa                                                                                                                         <  2:400518/147‑1 (MQ=255)                                         aaTTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACcgcg                                                                                                               >  2:425995/1‑151 (MQ=255)                                         aaTTAAACTTGGCGTAAAGGGCACCCTA‑AACCCAGCGGCCAGCGTCTGGGGATTGGGGGGGGGCAGGTGGGTGATAATGGCTAAAATGACCCCGCAGTGATGGGGGGCATAAAATAAAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACcgcg                                                                                                               <  1:975644/151‑1 (MQ=255)                                            tAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGcgtcg                                                                                                            >  1:314654/1‑151 (MQ=255)                                              aaCTTGGCGTAAATGGCACCCTA‑AACCCGGCGGCCAGCGTCGGGGGAATGGGGGGGTGCAGGGGGGTAATAAGGGCTAAAATGACCCCGCAGTGATGGGGGGCTTAAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTa                                                                                                          <  1:580973/151‑1 (MQ=255)                                                cTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAtt                                                                                                        >  2:906021/1‑151 (MQ=255)                                                 ttGGCGAAAATGCAACCCTA‑AACCCAGCGGCCAGCGTCGGGGGAAGGGGGGCGGGCAGCTGGTTAAAAATGGCTAAAACGACCCGGAGGTGATGGTGGCCTTAAAATAAAT‑TTTTTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTAtt                                                                                                        <  2:1220718/150‑1 (MQ=255)                                                  gggggaaaaggccacccTG‑AACCCAGGGTCCGGCGTCGGGGGAATGTGGGCGGGCAGGTGGGTGATAATGCCTAAAATGACCCCGCAGTGTTGGTGGCCCTGAAATATAA‑TTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc                                                                                                       <  1:147485/135‑1 (MQ=255)                                                   ggctaaaaGGGCAACCGA‑AATCCGGCGCCCAGCGTCGGGGGAATGGGGGCGTGCAGTGGGGTAATAATGGCTACAATGACCCCGCAGTGATGGTGGGCCTGAAATAAAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAtt                                                                                                        <  1:1280288/143‑1 (MQ=255)                                                   ggCGAAAAGGACACCCTG‑AACCCAGCGGCCAGCGTCGGGGAAATGGGGGCGGGCAGCGGGGTGATAATGTCTAAAAGGACCCCGCAGTGGGGGGGGCCCAAAAAAAAAG‑GTTCTGTCCCCGCCGTTAGCTAAAAAACCGCGTCGTAtt                                                                                                        <  1:1143481/148‑1 (MQ=255)                                                    gCGTAAATGGCACCCTG‑AACCCAGCGGCCAGCGTCGGGGAAATGGGGGGGTGCAGCTGGGTAATAATGGCTACAATGCCCCCGCAGTGAGGGGGGGCCTAAAAAATAA‑TTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                    <  1:142822/151‑1 (MQ=255)                                                    gCGTAAATGGCACCCTA‑AACCCAGCGGCCAGCGTCTGGGGAATGGGGGCGTGCAGGTGGGTGATAATGGCTACAATGCCCCCGCAGTGATGGTGGGCCTAAAATATAA‑TTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                    <  1:754351/151‑1 (MQ=255)                                                        aaagggaacccgg‑AACCCAGCGGCCAGGGTCTGGGAATGGGGGGCGTGCAGTTGGGTAATATTGGTTACAAGGCCCCGGCGGTGATGGTGGGCTTAAAATATAG‑GTTCTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                   <  2:985173/136‑1 (MQ=255)                                                        aaaTGGGCACCTGAAACCCAGCGCCCAGGGTGGGGGAAGGGTGGGCGGGCGGTTGGGTGATAATGCCTAAAATGACCCCGCGGTAGGGGTGGGCCTAAAATATAA‑TTTCTGAGCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                    <  2:688321/148‑1 (MQ=255)                                                        aaaTGGCACCCTG‑AACCCGGCGGCCAGCGTCTGGGGAATGGGGGCGTGCAGCGGGGTGATAATGGCTACAATGCCCCCGCGGTGATGGTGGCCCTGAAATATAA‑TTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg                                                                                                <  1:760013/151‑1 (MQ=255)                                                        aaaTGGCACCCTA‑AACCCGGCGCCCGGCGTTGGGGAAATGGGGGGGGGCGGTGGGGTAAAAATGGCTAAAATACCCCCGAGGTGATGGT‑GCCCAAAAAAATAT‑TTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc                                                                                                 <  1:44401/149‑1 (MQ=255)                                                        aaaTGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg                                                                                                >  1:1018466/1‑151 (MQ=255)                                                                 ccTG‑AACCCAGCGGCCGGCGTGTGGGAAATGGGGGCGGGCAGCGGGGTTATAATGGCTAAAAGGCCCCCGCGGGGAGGGTGGGCCTAAAATATAA‑GTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa                                                                                           <  1:484274/147‑1 (MQ=255)                                                                  cgagattCCAGGGGCCAGCTTCTGGGGAATGGGGGCGTGCGGGTGGTTAAAAAGGGCTGAAAGCACACCGTAAGGATGGTGGCCTTAAAATAAAT‑TTTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACa                                                                                         <  1:1077984/143‑1 (MQ=255)                                                                      aacccAGGGCCCAGCTTCGGGGGAAGGGGGGGGTCAAGCGGGTTAATAATGGCAACAATCACCCCGTTGTGATGGTGGGCTTAAAAAAAAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa                                                                                    <  2:321538/147‑1 (MQ=255)                                                                      aacccAGCGGCCAGCGTGTGGGGAATGGGGGCGGGCAGCTGGGTGATAAGGGCTAAAATGACCCCGCAGTGATGGTGGGCCTAAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGaa                                                                                   <  1:1134688/148‑1 (MQ=255)                                                                       acccggCGGCCAGCGTCGGGGAAATGGGGGCGTGCAGCTGGGTGATAATGCCTACAATGACACCGGAGTGTTGG‑GGGCCTTAAAAAAAATGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                                        <  1:532171/140‑1 (MQ=255)                                                                            gCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt                                                                              >  1:872636/1‑150 (MQ=255)                                                                             cGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                            >  1:327074/1‑151 (MQ=255)                                                                             cGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                            >  2:142747/1‑151 (MQ=255)                                                                                ccAGCGTCGGGGGAATGTGGGCGTGCAGCTGGGTGATAAGGGCAACAATGACCCCGCGGTGAGGGGGGCCTTAAAATATAA‑TTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                          <  2:1037379/150‑1 (MQ=255)                                                                                     gTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGcc                                                                    >  1:1221079/1‑151 (MQ=255)                                                                                      tggggggATGGTGGGGGTGCAGCTGGTTAATAATGTCTAAAATGCCCCCGCAGTGTGGGGGGGCTTGAAATATAT‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                                     <  2:572877/146‑1 (MQ=255)                                                                                      tcgggggATTGTGGGGGTGCGGCGGGGTAAAAATGCCTAAAATGACCCCGAAGTGAGGGGGGGAATGAAATAAAA‑GTTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt                                                                       <  1:350978/144‑1 (MQ=255)                                                                                         ggggAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt                                                                 >  2:795908/1‑150 (MQ=255)                                                                                         ggggAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt                                                                 >  2:458727/1‑150 (MQ=255)                                                                                           ggAATGGGGGCGGGGAGCTGGGTAAAATTGACTAAATTGACCCCGAAGTAAGGGTGGGCATAAAAAATAT‑TTTCTGATTCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttttt                                                                       <  2:481502/142‑1 (MQ=255)                                                                                                  ggggGTGCAGGTGGGGAAAAAGGGCTACAATGCCCCCGCAGTGATGGGGGGCCTGAAATATAA‑GTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt                                                        <  1:31876/150‑1 (MQ=255)                                                                                                      gTGCAGGGGGGTGATAATGGCTACAATGACCCGGCAGTGATGGTGGGCTTAAAATATAA‑TTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgt                                                     <  2:282891/149‑1 (MQ=255)                                                                                                               gggaaaaaggcataaaaTGCCCCCGGAGGGTTGGGGGGCTTAAAATAAAT‑TTTCGGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTg                                          <  2:969716/136‑1 (MQ=255)                                                                                                                gtaaTAATGGCTACAATCACACCGGGGAGATGGTGGGCCAGAAAAATAT‑TTGTTGATTCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCaa                                             <  2:86163/144‑1 (MQ=255)                                                                                                                 tGATAATGGCTACAATGACCCCGCAGTGTGGGTGGGCCTAAAATATAA‑TTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGc                                         <  1:799575/150‑1 (MQ=255)                                                                                                                   aaaaTGGCTAAAATGCCCCCGCAGTGATGGGGGGCTTAAAATATAT‑TTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGc                                         <  2:947393/146‑1 (MQ=255)                                                                                                                   aTAATGGCTACAATGACCCCGCAGTGATGGTGGGCCTAAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCg                                       <  1:10051/150‑1 (MQ=255)                                                                                                                    tAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg                                     >  2:479937/1‑151 (MQ=255)                                                                                                                     aagggCTAAAATCACCCCGAGGTGATGGGGGGCTTAAAATATAA‑TTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg                                     <  2:486911/147‑1 (MQ=255)                                                                                                                            aCAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAc                               >  2:1274268/1‑149 (MQ=255)                                                                                                                              aaTGACCCCGCAGTGAGGGTGGGCCTAAAATATAA‑TTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAg                             <  1:1129662/149‑1 (MQ=255)                                                                                                                               aTGACCCCGAAGTGAGGGGGGGCAAAAAATATAG‑GTTCTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAt                          <  2:566264/151‑1 (MQ=255)                                                                                                                                tGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAtt                         >  1:436688/1‑151 (MQ=255)                                                                                                                                 gACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAttt                        >  2:213490/1‑151 (MQ=255)                                                                                                                                   cacCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGc                      >  2:1096703/1‑151 (MQ=255)                                                                                                                                      cGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCggg                   >  2:205833/1‑151 (MQ=255)                                                                                                                                       gcggtggtggtgGGTTTAAAATATAA‑TTTCTGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAttt                        <  2:200866/142‑1 (MQ=255)                                                                                                                                          ttGATGGGGGGCTTAAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt                  <  2:314649/147‑1 (MQ=255)                                                                                                                                              tggtggGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgc           >  1:1104963/1‑151 (MQ=255)                                                                                                                                               ggtggCCCTAAAATAAAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCtt               <  1:1219494/146‑1 (MQ=255)                                                                                                                                                     ccTAAAATATAA‑GTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAcc    <  1:818065/151‑1 (MQ=255)                                                                                                                                                        gAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTg  >  1:1179423/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                       GTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTG  >  NC_000913/1222083‑1222382

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