Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	666,102	A→C	6.8%	M990L (ATG→CTG) ‡	yfhM →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	666,102	0	A	C	6.8%	185.3 / 4.8	59	M990L (ATG→CTG) ‡	yfhM	lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (30/25);  new base C (2/2);  total (32/27)
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.40e-01

ACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGTCGAAATCACCAC  >  NC_002947/665954‑666250                                                                                                                                                     |                                                                                                                                                     aCCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCa                                                                                                                                                      >  2:130526/1‑149 (MQ=255)      acacTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGcc                                                                                                                                                 <  6:336463/149‑1 (MQ=255)         cTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGaa                                                                                                                                              <  2:312014/149‑1 (MQ=255)          tACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAg                                                                                                                                             >  3:36650/1‑149 (MQ=255)            cGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCCATGCCGAAgtg                                                                                                                                           <  2:275357/149‑3 (MQ=255)             ggTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGc                                                                                                                                           >  3:201489/1‑148 (MQ=255)             ggTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCa                                                                                                                                          <  1:130526/149‑1 (MQ=255)                    aTAGTCGCCCAGCAGCGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGt                                                                                                                                   >  1:219252/1‑149 (MQ=255)                    aTAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGt                                                                                                                                   >  2:108426/1‑149 (MQ=255)                       gTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTgg                                                                                                                                 <  8:224859/148‑1 (MQ=255)                          gCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTcgc                                                                                                                             >  4:194013/1‑149 (MQ=255)                           cccAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTcgcc                                                                                                                            >  7:141942/1‑149 (MQ=255)                            ccAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTcgccg                                                                                                                           >  7:276730/1‑149 (MQ=255)                                cagAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTcgccgcgcc                                                                                                                       >  5:293756/1‑149 (MQ=255)                                cagAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTcgccgcgcc                                                                                                                       <  3:194013/149‑1 (MQ=255)                                   agCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGtt                                                                                                                    <  2:206152/149‑1 (MQ=255)                                         ccGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGAtt                                                                                                                <  8:141942/147‑1 (MQ=255)                                                               gTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTc                                                                                        <  3:325919/149‑1 (MQ=255)                                                                tGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTcc                                                                                       >  2:25210/1‑149 (MQ=255)                                                                    ggCCAGGCTACGGTCGACATCCCCGACTTCAAAGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACATCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCAGGc                                                                                   >  7:211617/1‑149 (MQ=255)                                                                     gCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGcc                                                                                  <  4:201489/149‑1 (MQ=255)                                                                     gCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGACGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGcc                                                                                  >  7:9755/1‑149 (MQ=255)                                                                      ccAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGc                                                                                                         <  7:249254/125‑1 (MQ=255)                                                                      ccAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGc                                                                                                         >  8:249254/1‑125 (MQ=255)                                                                      ccAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCg                                                                                 >  1:83442/1‑149 (MQ=255)                                                                          gCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTgg                                                                             <  8:64722/149‑1 (MQ=255)                                                                                   gACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTAc                                                                    <  3:19464/149‑1 (MQ=255)                                                                                   gACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTAc                                                                    <  2:278484/149‑1 (MQ=255)                                                                                   gACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTAc                                                                    >  7:90804/1‑149 (MQ=255)                                                                                         cccGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCt                                                              >  3:116692/1‑149 (MQ=255)                                                                                          ccGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTg                                                             >  3:20695/1‑149 (MQ=255)                                                                                          ccGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTg                                                             <  2:83442/149‑1 (MQ=255)                                                                                          ccGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTg                                                             >  1:88014/1‑149 (MQ=255)                                                                                          ccGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTg                                                             >  1:52555/1‑149 (MQ=255)                                                                                           cGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCTGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTgg                                                            <  2:219252/149‑1 (MQ=255)                                                                                               ttCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTg                                                                              <  6:86945/127‑1 (MQ=255)                                                                                                 cAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTg                                                       <  4:36650/148‑1 (MQ=255)                                                                                                     ggCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCt                                                  >  7:77667/1‑149 (MQ=255)                                                                                                     ggCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCt                                                  >  2:94593/1‑149 (MQ=255)                                                                                                       cGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTgg                                                <  5:192289/149‑1 (MQ=255)                                                                                                             gCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAAcc                                          >  7:213370/1‑149 (MQ=255)                                                                                                              cGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCt                                         <  8:276730/149‑1 (MQ=255)                                                                                                                   gATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCg                                    >  4:186890/1‑149 (MQ=255)                                                                                                                       gcgcAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCt                                         <  3:306545/140‑1 (MQ=255)                                                                                                                       gcgcAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCt                                         >  4:306545/1‑140 (MQ=255)                                                                                                                        cgcAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCcgc                               <  4:116692/149‑1 (MQ=255)                                                                                                                               cTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAg                         >  8:270626/1‑148 (MQ=255)                                                                                                                                catgccGAAGTGCTCTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGc                        >  2:140165/5‑148 (MQ=255)                                                                                                                                    ccGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACt                    >  1:83230/1‑148 (MQ=255)                                                                                                                                    ccGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACt                    >  3:173185/1‑148 (MQ=255)                                                                                                                                    ccGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACt                    >  2:67042/1‑148 (MQ=255)                                                                                                                                     cGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAAc                                           <  6:338994/124‑1 (MQ=255)                                                                                                                                     cGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAAc                                           >  5:338994/1‑124 (MQ=255)                                                                                                                                       aaGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGc                <  4:173185/149‑1 (MQ=255)                                                                                                                                          ttgccTTCGGCCATGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGTc             >  2:203702/5‑149 (MQ=255)                                                                                                                                          agCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGTc             <  8:77667/149‑1 (MQ=255)                                                                                                                                           tgccTTCGGCCATGCCGAAGTCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGt              <  4:231829/144‑1 (MQ=255)                                                                                                                                           tgccTTCGGCCATGCCGAAGTCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGt              >  3:231829/4‑147 (MQ=255)                                                                                                                                                 cGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGGAACTGAGCGTCTAAAt        <  1:67042/147‑1 (MQ=255)                                                                                                                                                    cATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGTCGAAATcacca   >  3:204996/1‑149 (MQ=255)                                                                                                                                                     aTGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGTCGAAATcaccac  >  6:68109/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ACCCAACACTACGGTAACCATAGTCGCCCAGCAGAGCCTGCCGGTGACGCTGGACGACAAGGGTGAGGGCCAGGCTACGGTCGACATCCCCGACTTCAACGGCGAATTGCGGCTGATGGCGCAGGCCTGGACCGAAGAGCACTTCGGCATGGCCGAAGGCAAGACCGTGGTCGCCGCGCCGTTGATTGCCGAGCTTTCGGCCCCGCGCTTCCTGGCCGGTGGCGACCGTACCAGCCTGGCACTGGACCTGGCCAACCTGTCCGGCCGCGCCCAGCAACTGAGCGTCGAAATCACCAC  >  NC_002947/665954‑666250

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