Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,213	0	A	C	45.3%	106.8 / 10.4	64	pseudogene (40/75 nt)	ycgH	pseudogene; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base A (27/8);  new base C (0/29);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 4.46e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGC  >  NC_000913/1222069‑1222362                                                                                                                                                  |                                                                                                                                                       tGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGc                                                                                                                                                    >  1:1419226/1‑149 (MQ=255)  gTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAg                                                                                                                                                  <  2:1210065/150‑1 (MQ=255)   tctcAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTg                                                                                                                                                >  2:1355094/1‑151 (MQ=255)    ctcAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTg                                                                                                                                                >  1:282894/1‑150 (MQ=255)      cAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTg                                                                                                                                                <  1:1440705/148‑1 (MQ=255)      cAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGAt                                                                                                                                              <  1:503501/150‑1 (MQ=255)        aTGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGAtggt                                                                                                                                           <  2:1252958/151‑1 (MQ=255)        aTGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGAtggt                                                                                                                                           >  1:1191833/1‑151 (MQ=255)        aTGAACGGTGTGAAAGACGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGAtgg                                                                                                                                            >  2:11204/1‑150 (MQ=255)         tGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGAtggtg                                                                                                                                          >  2:1143450/1‑151 (MQ=255)              ggTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCt                                                                                                                                     >  1:962380/1‑151 (MQ=255)               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                                                        >  2:142747/1‑151 (MQ=255)                                                                                           cGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                        >  1:327074/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)                                                                                                                                                  |                                                                                                                                                       TGTCTCAATGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTG‑AATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGGTGGGCCTGAAATATAA‑GTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGC  >  NC_000913/1222069‑1222362

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