Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	918,981	A→G	5.9%	L40P (CTG→CCG)	opdC ←	histidine‑specific outer membrane porin D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	918,981	0	A	G	5.9%	202.6 / 2.3	67	L40P (CTG→CCG)	opdC	histidine‑specific outer membrane porin D
Reads supporting (aligned to +/- strand):  ref base A (29/34);  new base G (2/2);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

AGCATGGCATGAGCATCCAGGCCAAAACCGACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGCTCCTTGGATCTGACGGTG  >  NC_002947/918834‑919127                                                                                                                                                    |                                                                                                                                                    aGCATGGCATGAGCATCCAGGCCAAAACCGACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAg                                                                                                                                                    <  1:222669/149‑1 (MQ=255) aGCATGGCATGAGCATCCAGGCCAAAACCGACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAg                                                                                                                                                    <  2:154255/149‑1 (MQ=255)   catggcatgAGCATCCAGGCCAAAACCGACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGt                                                                                                                                                  >  8:266610/1‑149 (MQ=255)           gAGCATCCAGGCCAAAACCGACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCGGGTTG‑AAGGAc                                                                                                                                         >  4:118902/1‑149 (MQ=255)              cATCCAGGCCAAAACCGACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGaa                                                                                                                                      <  3:145443/149‑1 (MQ=255)               aTCCAGGCC‑AAACCGCCCACACCCGGGGTGTAGACCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTAGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAAcc                                                                                                                                    >  4:326300/1‑149 (MQ=255)                  cAGGCCAAAACCGACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAg                                                                                                                                                    <  7:321175/132‑1 (MQ=255)                  cAGGCCAAAACCCACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAg                                                                                                                                                    >  8:321175/1‑132 (MQ=255)                   aGGCCAAAACCGACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTc                                                                                                                                 >  7:201632/1‑149 (MQ=255)                   aGGCCAAAACCGACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTc                                                                                                                                 <  3:319845/149‑1 (MQ=255)                     gCCAAAACCGACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGa                                                                                                                               <  4:71562/149‑1 (MQ=255)                     gCCAAAACAGACCACCCCCGGGGTGGAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGa                                                                                                                               <  4:66426/149‑1 (MQ=255)                        aaaaCCGACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑aa                                                                                                                                             <  4:237829/132‑1 (MQ=255)                        aaaaCCGACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑aa                                                                                                                                             >  3:237829/1‑132 (MQ=255)                         aaaCCGACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGaa                                                                                                                           <  1:288566/149‑1 (MQ=255)                            ccGACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGc                                                                                                                         >  4:51147/1‑148 (MQ=255)                                    acCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCt                                                                                                                >  2:319279/1‑149 (MQ=255)                                       cGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGt                                                                                                             >  4:103972/1‑149 (MQ=255)                                        ggggTGTAGCTCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTc                                                                                                            <  1:319279/149‑1 (MQ=255)                                             gtAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCa                                                                                                        <  2:43407/148‑1 (MQ=255)                                             gtAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCa                                                                                                        >  1:43407/1‑148 (MQ=255)                                                 cccGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGc                                                                                                   <  3:258412/149‑1 (MQ=255)                                                  ccGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCGGGTTG‑AAGGACGAAc                                                                                                                                     >  1:130177/1‑114 (MQ=255)                                                  ccGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCGGGTTG‑AAGGACGAAc                                                                                                                                     <  2:130177/114‑1 (MQ=255)                                                  ccGAATCGAAACGCGCTACCGCGCCCTGGCCCCAGTCGCGTGTGTCGCGCACGCCGGCGCTGTGTTTGTCGCGTTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGcc                                                                                                  <  8:317794/149‑1 (MQ=255)                                                    gACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTAGATGAAGCCATCAGCCTTGTCTTCATCTGCCTg                                                                                                >  3:205608/1‑149 (MQ=255)                                                     aCTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGAAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTgg                                                                                               >  3:70643/1‑149 (MQ=255)                                                         gAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGggcgg                                                                                           >  7:114049/1‑149 (MQ=255)                                                          aaaCGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCGGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGggcggc                                                                                          <  3:118902/149‑1 (MQ=255)                                                          aaaCGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCtt                                                                                                               <  6:306812/128‑1 (MQ=255)                                                          aaaCGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCtt                                                                                                               >  5:306812/1‑128 (MQ=255)                                                               cgcTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAg                                                                                     <  3:51147/149‑1 (MQ=255)                                                               cgcTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAg                                                                                     >  3:274452/1‑149 (MQ=255)                                                                     cgcgCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCg                                                                               >  6:136548/1‑149 (MQ=255)                                                                      gggcCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGt                                                                              <  1:238286/147‑1 (MQ=255)                                                                        gcCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTg                                                                            <  3:103972/149‑1 (MQ=255)                                                                          ccTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCa                                                                          <  8:114049/149‑1 (MQ=255)                                                                          ccTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCa                                                                          <  2:187752/149‑1 (MQ=255)                                                                            tGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGggcgg                                                                                           <  2:345051/130‑1 (MQ=255)                                                                            tGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGggcgg                                                                                           >  1:345051/1‑130 (MQ=255)                                                                                 ccACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACg                                                                   <  7:266610/149‑1 (MQ=255)                                                                                 ccACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACg                                                                   >  5:95226/1‑149 (MQ=255)                                                                                    cTCGCGGGTGTCGCGGACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCaa                                                                <  3:326300/149‑1 (MQ=255)                                                                                    cTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTAGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCaa                                                                <  4:205608/149‑1 (MQ=255)                                                                                       gcgGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAg                                                                                                                                            >  3:306280/1‑70 (MQ=255)                                                                                        cgGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCTATGAATCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCATTGaa                                                            <  5:265619/149‑1 (MQ=255)                                                                                        cgGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGaa                                                            <  8:201632/149‑1 (MQ=255)                                                                                           gtgtCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGATGTAGGCATTACGAAAATGCAGGTTGAAAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGaa                                                            >  1:112850/1‑147 (MQ=255)                                                                                             gtCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCg                                                       >  6:176308/1‑149 (MQ=255)                                                                                                                  ttggTCGCCGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACCAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAg                                  <  4:70643/146‑1 (MQ=255)                                                                                                                  gTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAg                                  <  6:95226/149‑1 (MQ=255)                                                                                                                      tCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCAc                              <  2:112850/149‑1 (MQ=255)                                                                                                                        gcgGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAGAGCGGAAGGCACCg                            >  8:134084/1‑149 (MQ=255)                                                                                                                                aGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGc                    >  4:349723/1‑149 (MQ=255)                                                                                                                                  tAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGCTc                  >  3:238206/1‑149 (MQ=255)                                                                                                                                  tAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGCTc                  >  3:243099/1‑149 (MQ=255)                                                                                                                                     gCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGCTCCtt               >  5:375454/1‑149 (MQ=255)                                                                                                                                       aTTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGCTCCTTg              >  8:94656/1‑148 (MQ=255)                                                                                                                                          aCGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGCTCCTTGGAt           <  5:176308/148‑1 (MQ=255)                                                                                                                                           cGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGCTCCTTGGATCt         <  3:349723/149‑1 (MQ=255)                                                                                                                                           cGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGCTCCTTGGATCt         <  5:136548/149‑1 (MQ=255)                                                                                                                                            gAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGCTCCTTGGAt           >  7:196771/1‑146 (MQ=255)                                                                                                                                            gAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGCTCCTTGGAt           <  8:196771/146‑1 (MQ=255)                                                                                                                                            gAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGACGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGCTCCTTGGATCt         >  2:309847/1‑148 (MQ=255)                                                                                                                                              aaaTGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGCTCCTTGGATCTGAc      <  4:274452/149‑1 (MQ=255)                                                                                                                                                  gCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGCTCCTTGGATCTGACGGTg  >  4:219281/1‑149 (MQ=255)                                                                                                                                                   cAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGCTCCTTGGATCTGACGGTg  <  7:94656/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    AGCATGGCATGAGCATCCAGGCCAAAACCGACCACACCCGGGGTGTAGCCCGACTCGAAACGCGCTACCGCGCCCTGGCCCCACTCGCGGGTGTCGCGCACGCCGGCGTTGTGGTTGTCGCGGTTGAAGTAGGCATTACGAAAATGCAGGTTG‑AAGGACGAACCTTCGATGAAGCCATCAGCCTTGTCTTCATCTGCCTGGGCGGCGAAAGGGATCGTTGCAGTCAACGCAATGAAAAGCGGGGTAAAACGAAACGCGGAAGGCACCGGTACTAGCTCCTTGGATCTGACGGTG  >  NC_002947/918834‑919127

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