Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	907,602	0	A	C	45.0%	84.4 / 13.1	60	G227G (GGT→GGG)	ybjT	putative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (23/10);  new base C (0/27);  total (23/37)
Fisher's exact test for biased strand distribution p-value = 4.07e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGC  >  NC_000913/907458‑907750                                                                                                                                                   |                                                                                                                                                     cgcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTc                                                                                                                                                <  2:935409/151‑1 (MQ=255)  gcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCa                                                                                                                                               <  1:896994/151‑1 (MQ=255)     gTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCg                                                                                                                                            <  1:223570/151‑1 (MQ=255)      tATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGc                                                                                                                                           >  1:568501/1‑151 (MQ=255)         catcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCa                                                                                                                                         <  1:237444/150‑1 (MQ=255)                    agATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCa                                                                                                                              >  2:911185/1‑150 (MQ=255)                     gttCGTGTTTCAGCCCCAAAACCAACGCCCTGGGGGTGGGGGGGGGTGGGGAAGTAATCACATTGAAACCCCACCCCGAACCCCAGCGGGGGGGGAGGGGG‑AGGGGAA‑TCACCCAGCGGCGTTTCCCCCTCACCGCCATAAAATGTTCa                                                                                                                              <  1:440599/147‑1 (MQ=255)                          tgtTCCAGCCCCGGAATCACCGCCCGGGGGGGGGTGGGCGGTAGGGAAGTAATCACATTGAAAACCCCCACCAAAACCCAGCGGGGGGGGAGGGGG‑GGGGGAA‑TCAACCAGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTg                                                                                                                         <  1:66832/149‑1 (MQ=255)                           gtTTCAGCCCCTGAATCACCGCCCGGGGGGTGGGGGGGGGTACGGAAGTAATCACATTGAGAACCCACACCAAAACCCAGGGGGTGGGGGGGGGG‑GTGGGGA‑TCAACCGGCGGCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGtt                                                                                                                       <  1:932027/150‑1 (MQ=255)                               cAGCCCCTGAACAAACCCCCGGGGGGTGGTGGGCGGTCGGAAAGTAATCACATGGAAAACCCACCCCAAACCCCGGGGGGGGGGGGGGGGG‑AGGGGGA‑AAAACCAGGGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTg                                                                                                                   <  2:939877/150‑1 (MQ=255)                                           cAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCAcc                                                                                                      >  2:383556/1‑151 (MQ=255)                                           cAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCAcc                                                                                                      >  2:489433/1‑151 (MQ=255)                                            accGCCCGGGCGGGGGGGGGGGGAGCGAAAGTAACCACTTTGAAAACCCCCCCCAAAACCCGGGGGGGGGGGGGGGGG‑AGGGGGA‑TAACCCGGGGGCTTTTCCCGCTCACCGCCAT‑AAATGTTCAAACTGTTGCTGATAACTGAGCAc                                                                                                       <  2:241391/146‑1 (MQ=255)                                            aaCGCCCTGGCGCTGGTGGGCGGTACGGAAGTAATCACATTGCGACACCCCACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCAcc                                                                                                      >  1:418685/1‑150 (MQ=255)                                              cGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctc                                                                                                    >  1:22642/1‑150 (MQ=255)                                              cGCCCTGGCGGTGGTGGGCGGTACGGAAATAAACCCTTTGAAAAACCCCACCGAAAGCCAGGGGGTGGGGGGGGGA‑AGGGGGA‑TCACCCAGGGGCGCTCCCCCTTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACct                                                                                                     <  2:532019/149‑1 (MQ=255)                                              cGCCCGGGCGGTGGGGGGGGGTAGGGAAGTAATCACATGGAAAACCCACCCCAAAATCCGGCGGGGGGGGGGGGGG‑AGGGGGA‑TCAACCAGCGGCGCTTCCCGCCCACCGCCATAAAATGATCAAACTGTTGCTGATAACTGAGCACctc                                                                                                    <  2:168711/150‑1 (MQ=255)                                                cccGGGGGGTGGGGGGGGGTGGGAAAGAAATCACAGTGAGAACCCACCCCAAAACCCAGGGGGGGGGGAGGGGG‑AGGGGAA‑TAACCCGGCGGCGCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctct                                                                                                   <  1:595868/149‑1 (MQ=255)                                                cccGGGGGGTGGGGGGGGGTAGGAAAGTAATCCCATTGAGAACCCACCCCAAAACCCGGGGGGGGGGGAGGGGG‑ATGGGAA‑TCAACCAGCGGCGCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTg                                                                                                  <  1:8450/150‑1 (MQ=255)                                                  cTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAAATGTTGCTGATAACTGAGCACCTCTGGTc                                                                                               >  1:19857/1‑151 (MQ=255)                                                       ggggggggggggTAGGGAATTAACCACATGGAAAACCCCCCCCAAAATCCGGGGGGGGGGGGGGGGG‑ATGGGGA‑CCACCCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggc                                                                                           <  2:823282/144‑1 (MQ=255)                                                             ggggggTACGGAAGAAATCACATAGAAAACCCACCCCAAAACCCGGGGGGTGGGGGGGGGG‑AGGGGGA‑TCACCCGGGGGCGCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTc                                                                                     <  2:493009/150‑1 (MQ=255)                                                                 ggcaCGAAAGTAACCAATTGGAAAAACCCCCCAAAAATCCAGCGGGGGGGGAGGGGG‑ATGGGAA‑TCACCCAGCGGCCTTAACCGTCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAg                                                                                 <  2:296937/147‑1 (MQ=255)                                                                    aCGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGc                                                                             >  1:16497/1‑151 (MQ=255)                                                                     cGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCg                                                                            >  1:81916/1‑151 (MQ=255)                                                                       gAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgt                                                                          >  2:983870/1‑151 (MQ=255)                                                                        aataaaccacaTGGAAAACCCCCCCCAAAACCCGGGGGGGGGGGGGGGGG‑AGGGGGA‑TCACCCGGCGGCGTTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCg                                                                            <  2:325559/141‑1 (MQ=255)                                                                         agaaaccacaTAGAAAACCCACCCCAAATCCCGGGGGGGGGGGGGGGGG‑AGGGGGA‑TCACCCAGCGGCGTTTACCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtgt                                                                        <  2:381064/145‑1 (MQ=255)                                                                          gTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGtt                                                                       >  2:441287/1‑151 (MQ=255)                                                                                aaagtGAGAACCCCCACAAAAACCCGGGGGGTGGGGAGGGGG‑AGGGGGA‑TCAACCAGCGCCGCTTACCTCCCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCt                                                                   <  1:303548/145‑1 (MQ=255)                                                                                   tagagaACCCCCCCCAAAACCCGGGGGGGGGGGGGGGGG‑AGGGGAA‑CCACCCGGCGGCGCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGcc                                                               <  1:963866/148‑1 (MQ=255)                                                                                    tGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCgg                                                             >  2:765468/1‑151 (MQ=255)                                                                                      agaAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgat                                                           >  2:212854/1‑151 (MQ=255)                                                                                       aaacccACACCAAAACCCAGGGGGGGGGGGGGGGG‑AGGGGGA‑TCAACCGGGGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatg                                                          <  2:326019/147‑1 (MQ=255)                                                                                        aacccccccagAAACCCGGGGGGGGGGGGGGGGG‑GGGGGGA‑CCACCCGCGGGCGCTACCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCTCCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCg                                                              <  1:955003/136‑1 (MQ=255)                                                                                         aacccccaccaaaaCCAGGGGGGGGGGGGGGGGAATGGGGA‑TCAACCAGCGCCCTTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatga                                                         <  2:430585/141‑1 (MQ=255)                                                                                           ccccacCAAAACCCGGGGGGGGGGGGGGGGG‑AGGGGGA‑TAAACCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatgat                                                        <  2:845310/146‑1 (MQ=255)                                                                                           ccACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCt                                                      >  1:934316/1‑151 (MQ=255)                                                                                            cccacCGAAACCCGGGGGGGGGGGGGGGGG‑AGGGGGA‑CCACCCAGCGGCGTTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTa                                                     <  1:213969/149‑1 (MQ=255)                                                                                             cacccaaaatgccaggggggggggggggg‑AAGGGGACACAACCGGCGCCGCTCCCCGCCCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTa                                                     <  2:166032/136‑1 (MQ=255)                                                                                              cacCGAAACCCGGGGGGGGGGGGGGGGG‑AGGGGAA‑TCAACCAGGGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTa                                                     <  2:383084/149‑1 (MQ=255)                                                                                                   aaccccggCGGGGGGGGAGGGGG‑GAGGGGA‑TCACCCAGGGCCATTCCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGc                                              <  2:568593/144‑1 (MQ=255)                                                                                                    aaccgagggggggggggggggg‑AGGGGAC‑CCACCCAGGGCCGTTTCCCGTCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACg                                               <  1:1003601/141‑1 (MQ=255)                                                                                                     aTCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCa                                             >  2:331161/1‑150 (MQ=255)                                                                                                       ccAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAcc                                          >  1:746301/1‑151 (MQ=255)                                                                                                          gggggggggggggggg‑AGGGGGA‑TCAACCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCa                                         <  2:365830/147‑1 (MQ=255)                                                                                                            gggggggggggggg‑AGGGGAA‑TCAACCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATa                                     <  2:710751/151‑1 (MQ=255)                                                                                                            gggTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGAt                                      >  1:842114/1‑150 (MQ=255)                                                                                                            gggTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGAt                                      >  2:683429/1‑150 (MQ=255)                                                                                                                        gg‑atgggga‑tCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtt                         >  1:969267/1‑151 (MQ=255)                                                                                                                             gggga‑tCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGttt                        <  2:111115/148‑1 (MQ=255)                                                                                                                              ggga‑tCAACCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCa                    <  1:489344/151‑1 (MQ=255)                                                                                                                              ggga‑cCACCCGGCGGCTCTCCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCa                    <  1:507718/151‑1 (MQ=255)                                                                                                                               gga‑tCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTcc                     >  2:232405/1‑149 (MQ=255)                                                                                                                                    cacccggcggcgCTTCCCCCCCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgcgc                <  2:218794/144‑1 (MQ=255)                                                                                                                                        cggcggctCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATgg           <  1:911051/149‑1 (MQ=255)                                                                                                                                         aGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATggg          >  2:640158/1‑151 (MQ=255)                                                                                                                                         aGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATggg          >  2:890330/1‑151 (MQ=255)                                                                                                                                                 ttACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCgt       <  2:838980/146‑1 (MQ=255)                                                                                                                                                 ttACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGc  >  2:74105/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAATCCAGCGGGTGGGGAGGGGG‑ATGGGGA‑TCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGC  >  NC_000913/907458‑907750

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