Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	462,067	C→T	S51S (AGC→AGT)	ppiD →	periplasmic folding chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	462,067	0	C	T	97.3%	104.2 / NA	37	S51S (AGC→AGT)	ppiD	periplasmic folding chaperone
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/22);  minor base G (0/1);  total (14/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CCGTTTCTTGTCACAATAGGCCTTTGCGCGCATCGATACGTTGCGTGAGGTACACAGTCATCTACAGCGGAGTGTTGTTACACCATGATGGACAGCTTACGCACGGCTGCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGCCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTCGCGACCCCAGCCTTCCAGGTTGACGGCAAATT  >  NC_000913/461831‑462312                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      ccGTTTCTTGTCACAATAGGCCTTTGCGCGCATCGATACGTTGCGTGAGGTACACAGTCATCTACAGCGGAGTGTTGTTACACCATGATGGACAGCTTACGCACGGCCGCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCga                                                                                                                                                                                                                                         <  1:214171/251‑1 (MQ=255)              aaatgggCCTTTG‑GCGGATCGGTTCGTTTCGTGAGGGTCACAGTTATCTTCAACGGAGGGTTTGTACACCATGATGGTCATCCTACGCACGGCTGCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACa                                                                                                                                                                                                                           <  2:224383/246‑1 (MQ=255)                           cgcgcATCGATTCGTTGCGTGAGGTACACAGTCATCTACAGCGGAGTGTTGTTACACCATGATGGACAGCTTACGCACGGCTGCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATGGACCGGCGTGAGTGGTTCCCGTATTTGCGGATGCAATAACTACTCCGCAAAAGTGAATGCCCAGTAAACCAGTCGTGGGCATTGCTATCACGCCTTCAACAGCGATCGTAATc                                                                                                                                                                                                               <  2:255219/251‑1 (MQ=255)                                aTCGATACGTTGCGTGAGGTACACAGTCATCTACAGCGGAGTGTTGTTACACCATGATGGACAGCTTACGCACGGCTGCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATg                                                                                                                                                                                                          >  1:112513/1‑251 (MQ=255)                                      aCGTTGCGTGAGGTACACAGTCATCTACAGCGGAGTGTTGTTACACCATGATGGACAGCTTACGCACGGCTGCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTAACTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAACGAGCGTAATCGCATgcagca                                                                                                                                                                                                     >  1:43551/1‑250 (MQ=255)                                          tGCGTGAGGTACACAGTCATCTACAGCGGAGTGTTGTTACACCATGATGGACAGCTTACGCACGGCTGCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGGGGGCAATTCGAGAACCCCTTCAACAGCGAGCGTAATCGCCTGCCGAAACAGc                                                                                                                                                                                                >  2:193029/1‑251 (MQ=255)                                             gTGAGGTACACAGGCATCTTCAGCGGAGTGTTGTTACACCATGATGGACAGCTTACGCACGGCTGCAAACAGTCGCGTGCTCAAGATTATTTTCGGTATCATTATCGTGGCGTGCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGGCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTgg                                                                                                                                                                                             <  2:235903/251‑1 (MQ=255)                                                    acacAGTCATCTACAGCGGAGTGTTGTTACACCATGATGGACAGCTTACGCACGTCTTCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTTCCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGTGCGATCa                                                                                                                                                                                      <  1:187527/251‑1 (MQ=255)                                                    acacAGTCATCTACAGCGGAGTGTTGTTACACCATGAGGGACAGCTTACGCACGGCTGCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCa                                                                                                                                                                                      <  1:143328/251‑1 (MQ=255)                                                                                           gACAGCTTACGCACGGCTGCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAAc                                                                                                                                               <  1:126112/251‑1 (MQ=255)                                                                                                   aCGCACGGCTGCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGt                                                                                                                                        <  2:217780/250‑1 (MQ=255)                                                                                                                        cGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGa                                                                                                                  <  2:189837/251‑1 (MQ=255)                                                                                                                         gTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGAt                                                                                                                 >  1:123494/1‑251 (MQ=255)                                                                                                                              cAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTg                                                                                                                            <  2:3010/235‑1 (MQ=255)                                                                                                                              cAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTg                                                                                                                            >  1:3010/1‑235 (MQ=255)                                                                                                                              cAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATcgacg                                                                                                            <  1:106489/251‑1 (MQ=255)                                                                                                                                 gATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGAAGAgg                                                                                                         >  1:87520/1‑251 (MQ=255)                                                                                                                                    tattTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGgcg                                                                                                       <  1:85785/250‑1 (MQ=255)                                                                                                                                          cGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTgg                                                                                                <  1:106923/251‑1 (MQ=255)                                                                                                                                                       gtctCGTTCATATTGACCGGCGGTAGTGGTTTCCTGACTGGCGGAGGCAATATCTACGCCGCAGAAGTGAATGACCAGGAACTCATTCGTTGGCAATTCGAGAACGCTTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTTGCCGCGAACGAAGGCTATATGAAAACCCCGCGTCAACAGGTGCTGATTCTTCTGATCGAAGATGCGCTGCTGGATCAGTagccacg                                                                                   <  2:254016/248‑7 (MQ=255)                                                                                                                                                          tCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATg                                                                                                                                                    >  2:69622/1‑183 (MQ=255)                                                                                                                                                          tCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATg                                                                                                                                                    <  1:69622/183‑1 (MQ=255)                                                                                                                                                                    tGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGCACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTg                                                                      >  1:249472/1‑251 (MQ=255)                                                                                                                                                                      aCCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTggg                                                                    <  1:227131/251‑1 (MQ=255)                                                                                                                                                                           cgttaGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAATTGAATGACCAGTAAATCAGTCGTGGGCAATTCAAGAACGCCTTCAACAGCGAGCGTAATCGCTTGCGGCAACAGCTGGGCGATCAATGCTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGGCAACAGGTTCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCa                                                               <  1:254438/247‑1 (MQ=255)                                                                                                                                                                                     tACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGc                                                      >  1:20293/1‑250 (MQ=255)                                                                                                                                                                                                cGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGGTCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCGGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAgg                                          >  1:74360/1‑251 (MQ=255)                                                                                                                                                                                                       aaTAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGAtttt                                   <  1:177199/251‑1 (MQ=255)                                                                                                                                                                                                             tACGCCGCAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTCGCGAc                             >  1:33702/1‑251 (MQ=255)                                                                                                                                                                                                                   gcAAAAGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTCGCGACCCCAGc                       >  2:164021/1‑251 (MQ=255)                                                                                                                                                                                                                     aaaaGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTCGCGACCCCAGCCt                     <  2:123494/251‑1 (MQ=255)                                                                                                                                                                                                                       aaGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTCGCGACCCCAGCCTTc                   <  2:83091/251‑1 (MQ=255)                                                                                                                                                                                                                       aaGTGAATGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGAGTTTCGCGACCCCAGCCTTc                   <  2:112513/251‑1 (MQ=255)                                                                                                                                                                                                                            aaTGACCAGGAAATCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTCGCGACCCCAGCCTTCCAGGt              <  1:204641/251‑1 (MQ=255)                                                                                                                                                                                                                                      aaaTCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGGCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTCGCGACCCCAGCCTTCCAGGTTGACGGCaaa    <  1:193029/251‑1 (MQ=255)                                                                                                                                                                                                                                       aaTCAGTCGTGGGCAATTAGAGAACGCATTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTCGCGCCCCCAGCCTTCCAGGTTGACGGAAAAt   >  2:252314/1‑251 (MQ=255)                                                                                                                                                                                                                                        aTCAGTCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTCATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGGACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTCGCGACCCCAGCCTTCCAGGTTGAGGGCAAAtt  >  2:185996/1‑251 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                      CCGTTTCTTGTCACAATAGGCCTTTGCGCGCATCGATACGTTGCGTGAGGTACACAGTCATCTACAGCGGAGTGTTGTTACACCATGATGGACAGCTTACGCACGGCTGCAAACAGTCTCGTGCTCAAGATTATTTTCGGTATCATTATCGTGTCGTTCATATTGACCGGCGTGAGTGGTTACCTGATTGGCGGAGGCAATAACTACGCCGCAAAAGTGAATGACCAGGAAATCAGCCGTGGGCAATTCGAGAACGCCTTCAACAGCGAGCGTAATCGCATGCAGCAACAGCTGGGCGATCAATACTCCGAGCTGGCAGCGAACGAAGGCTATATGAAAACCCTGCGTCAACAGGTGCTGAATCGTCTGATCGACGAGGCGCTGCTGGATCAGTACGCACGTGAGCTGAAACTGGGTATCAGCGATGAGCAGGTTAAACAGGCGATTTTCGCGACCCCAGCCTTCCAGGTTGACGGCAAATT  >  NC_000913/461831‑462312

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