Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,233	0	G	A	49.2%	55.8 / 41.9	65	pseudogene (60/75 nt)	ycgH	pseudogene; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (29/4);  new base A (0/32);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.83e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGC  >  NC_000913/1222084‑1222375                                                                                                                                                        |                                                                                                                                               tgtgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTg                                                                                                                                                >  2:1183893/1‑150 (MQ=255)    gAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGaaa                                                                                                                                             <  1:1356833/150‑1 (MQ=255)     aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAata                                                                                                                                           <  2:1316540/151‑1 (MQ=255)         tCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAAtata                                                                                                                                         >  1:318259/1‑149 (MQ=255)                aGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGAt                                                                                                                                >  1:388166/1‑151 (MQ=255)                          aGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGtt                                                                                                                       >  1:1484755/1‑150 (MQ=255)                            tAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGc                                                                                                                    >  2:39387/1‑151 (MQ=255)                             aaTTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCt                                                                                                                   >  1:591070/1‑151 (MQ=255)                                            aaaCTTGGCGTAAAGGCCACCCTGAATCCGGCGGCCGGCGTCTGGGGAA‑TGTGGGTGGGCAGGTGGGTAAAATTGTTTAAAATGCCACGGCATTGTTGG‑TGGGCCGAAAATATAATTTTTAACCCCGCCGTTAGCTAAAAAACCGcgtcg                                                                                                     <  2:823922/150‑1 (MQ=255)                                              aCTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAt                                                                                                  >  2:1509660/1‑151 (MQ=255)                                                 tgacGTAAAGGCCACCCTAACCCCGGCGGCCCGGTTTTGGGAAT‑TGTGGCGGTGCATTTGGGTAATAATGCTTACAACGACCCCGGAGTGGTGG‑GGGGCTTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGcgtcgt                                                                                                    <  2:332494/143‑1 (MQ=255)                                                     gTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACCATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                             >  1:522397/1‑149 (MQ=255)                                                      tAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc                                                                                          >  1:1257389/1‑151 (MQ=255)                                                      ctctccctatttggAACCCGGCGGCCGGCGTCGGGGGAA‑TGGGGGGGTGCAGTGGGGTAAAAATGGCTACAATGACCCCGGAGTGATGG‑TGGGCGTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt                                                                                              <  1:516767/134‑1 (MQ=255)                                                       aaaTGGCACCCTGAACCCGGCGCCCAGCGTTTGGGGAA‑TGTGGGCGTGCAGCGGGGTAATAATGTTTACAATAACCCCGCAGTGTTGG‑TGGGCCTAAAATAAAATTTCTGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg                                                                                         <  2:1497735/151‑1 (MQ=255)                                                       aaaTGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg                                                                                         >  2:1138471/1‑151 (MQ=255)                                                       aaaGGGCACCCGAAACCCGGCGGCCGGCTTCGGGGGAA‑TGTGGGGGTGCAGCTGGGTAATAATGGCTACAATGCCCCCGCAGTGATGG‑GGGGCCTAAAATATAAGTTCTGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc                                                                                          <  2:499075/150‑1 (MQ=255)                                                        aaTGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc                                                                                        >  1:557117/1‑151 (MQ=255)                                                          tGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCgg                                                                                      >  1:480340/1‑151 (MQ=255)                                                           ggCACCCTGAATCCGGCGGCCGGCGTCTGGGGAA‑TGTGGGCGGGCAGCTGGGTGATAATGGTTACAATGACCCCGCAGTGATGG‑TGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGt                                                                                     <  2:1349129/151‑1 (MQ=255)                                                              acccTGAACCCAGCGGCCGGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTAATAATGGCTACAATCCCCCCGCAGTGAGGG‑GGGCCTTAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCgg                                                                                      <  1:945224/145‑1 (MQ=255)                                                                ccTGAAGCGAGCGCCCCGCGTTGGGGAAA‑GGTGGGCGTGCAGGTGGGGGAAAATGGCTAAAAGGACCCGGCAGTGTTGG‑GGGGCGAAAAATAAAAGTTTTCATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTa                                                                                    <  2:1117282/147‑1 (MQ=255)                                                                ccTAAATCCGGCGGCCGGCGTCGGGGGAT‑TGGGGGCGTGCTGTTGGGTGATAATGGCTACAATCACCCCGGAGTGTGGG‑GGGGCTTGAAATAAAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata                                                                                <  1:440398/151‑1 (MQ=255)                                                                 cTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat                                                                               >  2:880422/1‑151 (MQ=255)                                                                 cTGAACCCAGCGGCCGCCTTCGGGGGAA‑TGGGGGCGGGCGGTTGGGTAATATTGCCTAAAATGCCACCGGAGTGAGGG‑GGGCCTTAAAATATAATTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                                 <  1:783867/149‑1 (MQ=255)                                                                  tGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCAt                                                                                              >  2:516840/1‑135 (MQ=255)                                                                    aacccgggggCCGGCTTCTGGGGAA‑TGGGGGCGGGCAGGGGGGTGATAATGGCTACAATGACCCCGCGGTGATGG‑TGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                              <  2:1090756/141‑1 (MQ=255)                                                                       ccggCGGCCGGGGTCGGGGGAA‑TGTGGGGGGGCAGCTGGTTAATAATGGCAACAATGCCACCGGAGTGATGG‑TGGGCAAAAAATAAAAGTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAata                                                                         <  1:958070/148‑1 (MQ=255)                                                                       ccAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat                                                                          >  1:827284/1‑150 (MQ=255)                                                                        cAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAatat                                                                        >  1:824990/1‑151 (MQ=255)                                                                         aGCGGCCTGCGTCTGGGGAT‑GGTGGGGGTGCAGGTGGTTGATAATGGCTAAAGTGCCCCCGAAGTGTTGG‑TGGTCATAAAATATATTTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat                                                                               <  1:218473/143‑1 (MQ=255)                                                                           cGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtattat                                                                     >  2:594331/1‑151 (MQ=255)                                                                            ggCCGGTGTCGGGGGAA‑TGGGGGCGTACGGCTGGGTAAAAGTGACTAAAAGGACCCCACAG‑GTTGGTTGGGCCTAAAATATAGTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACat                                                                                 <  2:244279/138‑1 (MQ=255)                                                                             gCCGGCGTCAGGGAAG‑TGGGGGCGGGCAGCGGGGTAATAATGCATACAATACCACTGGGGGGGTGG‑GGGGTGAAAAAAATAATTTCGCGTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                   <  1:1258829/151‑1 (MQ=255)                                                                             gCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                   >  2:1080866/1‑151 (MQ=255)                                                                               cAGCGTCGGGGGAA‑TGTGGGCGTGCGGCTGGGTGATAAGGGCTAAAATGACCCCGAAGTGGTGG‑TGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAtttt                                                                  <  2:1384796/150‑1 (MQ=255)                                                                                 gCGTCTGGGGAG‑TGTGGGCGTGCAGCTGGGTGATAATGGCTATAATTACACCGCAGTGATGG‑GGGGCCTGAAATATAAGTTCTGATGCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                                 >  1:1172713/1‑149 (MQ=255)                                                                                 gCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTAAAATGACCCCGCAGTGAGGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                                 <  2:713666/149‑1 (MQ=255)                                                                                  cGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                              >  1:201678/1‑151 (MQ=255)                                                                                   gTCGGGGGAG‑TGGGGGCGGGCGGCTGGGTAATAATGGCTAAAATGACCCCGCAGTGAGGG‑GGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                              <  2:1408441/150‑1 (MQ=255)                                                                                     tgggggAAGTGTGGGCGTGCAGCTGGTTGATAAGGCCTAAAATCCCCCCGAGGGGGTGG‑TGCGCTTAAAATATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGc                                                              <  1:899215/147‑1 (MQ=255)                                                                                     cTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg                                                            >  2:142275/1‑150 (MQ=255)                                                                                       ggggAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt                                                          >  1:54305/1‑150 (MQ=255)                                                                                       ggggAA‑TGGGGGGGTGCAGGTGGGTGATAATGGCTAAAATGACCCCGCAGTGTGGG‑TGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATa                                                         <  2:825065/151‑1 (MQ=255)                                                                                       ggggAA‑TGGGGGCGTGCAGCTGGGTGATAATGGCTACAATGACCCCGCAGTGATGG‑GGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg                                                            <  2:1485489/148‑1 (MQ=255)                                                                                           aa‑TGTGGGCGTGCAGCGGGGTGATAATGGCTAAAATGACCCCGCAGTGATGG‑GGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCa                                                      <  2:346057/150‑1 (MQ=255)                                                                                              tgggggcgggcgtgtGGGTGAAAATGGCTAAAATGACCCCGCAGTGATGG‑GGGGCTTAAAATATAATTTTTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACg                                                    <  2:647145/136‑1 (MQ=255)                                                                                                  ggCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATgg                                               >  1:506213/1‑151 (MQ=255)                                                                                                    cGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg                                             >  2:1238077/1‑151 (MQ=255)                                                                                                      tGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtc                                           >  1:34278/1‑151 (MQ=255)                                                                                                                  aTAATGGCTAAAATGACACCGCGGTGATGG‑TGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGt                               <  2:1406062/151‑1 (MQ=255)                                                                                                                    aaTGGCTAAAATGACACCGCAGTGGTGG‑GGGGCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGa                             <  2:272072/151‑1 (MQ=255)                                                                                                                       ggCTAAAATGACCCCGTAGTGGTGG‑TGGTCTTAAAATAAAAGTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCa                          <  2:1121824/151‑1 (MQ=255)                                                                                                                           aCAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAc                        >  1:383658/1‑149 (MQ=255)                                                                                                                            caagccccccGCAGTGATGG‑GGGCCCTAAAATATATTTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACa                       <  2:806669/141‑1 (MQ=255)                                                                                                                             aaTGACCCCGCAGGGATGG‑GGGGCTTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACa                       <  2:577103/148‑1 (MQ=255)                                                                                                                             aaTGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                     >  1:128848/1‑150 (MQ=255)                                                                                                                                  ccccGCGGTGGGGG‑TGGGCCTAAAATATAAGTTCTGTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTg                <  2:309165/148‑1 (MQ=255)                                                                                                                                    ccGCGGTGATGG‑GGGGCATAAAATATATTTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCg              <  1:924338/150‑1 (MQ=255)                                                                                                                                     cGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg             >  1:1301720/1‑150 (MQ=255)                                                                                                                                      gaaGTGATGG‑GGGGCCTAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCggg            <  2:534527/148‑1 (MQ=255)                                                                                                                                       cAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTc          >  1:1450853/1‑151 (MQ=255)                                                                                                                                              gg‑tggGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTa      <  2:13886/148‑1 (MQ=255)                                                                                                                                              gg‑tggGCCTAAAAAATAATTTTGGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTa      <  2:1251534/148‑1 (MQ=255)                                                                                                                                                 tggGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcgc  >  2:477602/1‑150 (MQ=255)                                                                                                                                                  gggCCTAAAAAAAATGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcgc  <  1:263151/149‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                               TGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGGCAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTGATGG‑TGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGC  >  NC_000913/1222084‑1222375

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