Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,752,333	A→G	6.0%	L517P (CTC→CCC)	PP_3316 ←	chaperone‑associated ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,752,333	0	A	G	6.0%	201.6 / 2.2	68	L517P (CTC→CCC)	PP_3316	chaperone‑associated ATPase
Reads supporting (aligned to +/- strand):  ref base A (33/30);  new base G (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

CTTGCTGATCTGCGCCGCATTGTCATATTGCTTGCGCGAAGCCGCGTAGTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTCGCGCAGGCCGCGCAGGATCATGATG  >  NC_002947/3752185‑3752476                                                                                                                                                     |                                                                                                                                                  cTTGCTGATCTGCGCCGCATTGTCATATTGCTTGCGCGAAGCCGCGTAGTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGa                                                                                                                                                   <  7:4566/149‑1 (MQ=255)    gCTGATCTGCGCCGCATTGTCATATTGCTTGCGCGAAGCCGCGTAGTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGc                                                                                                                                                 <  7:107721/148‑1 (MQ=255)    gCTGATCTGCGCCGCATTGTCATATTGCTTGCGCGAAGCCGCGTAGTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGc                                                                                                                                                 <  3:106513/148‑1 (MQ=255)      tGATCTGCGCCGCATTGTCATATTGCTTGCGCGAAGCCGCGTAGTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAg                                                                                                                                                  >  5:220882/1‑145 (MQ=255)      tGATCTGCGCCGCATTGTCATATTGCTTGCGCGAAGCCGCGTAGTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAg                                                                                                                                                  >  4:109304/1‑145 (MQ=255)           tGCGCCGCATTGTCATATTGCTTGCGCGAAGCCGCGTAGTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGAt                                                                                                                                         >  7:116098/1‑149 (MQ=255)                            ttgcttgcGCGAAGCCGCGTAGTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCa                                                                                                                                                >  5:12295/1‑125 (MQ=255)                            ttgcttgcGCGAAGCCGCGTAGTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCa                                                                                                                                                <  6:12295/125‑1 (MQ=255)                                ttgcGCGAAGCCGCGTAGTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAgcg                                                                                                                   >  7:270318/1‑149 (MQ=255)                                  gcgcGAAGCCGCGTAGTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAgcgc                                                                                                                  <  1:256162/148‑1 (MQ=255)                                   cgcgAAGCCGCGTAGTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAgcgcg                                                                                                                 <  6:220882/148‑1 (MQ=255)                                      gAAGCCGCGTAGTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCt                                                                                                             >  2:96253/1‑149 (MQ=255)                                      gAAGCCGCGTAGTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTCATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCt                                                                                                             >  7:282802/1‑149 (MQ=255)                                                 gTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGt                                                                                                 <  3:109304/149‑1 (MQ=255)                                                    cTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCg                                                                                               >  6:138034/1‑148 (MQ=255)                                                    cTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATAGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCg                                                                                               >  6:13745/1‑148 (MQ=255)                                                       cTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGa                                                                                             <  5:138034/147‑1 (MQ=255)                                                              cgcAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTcggc                                                                                    >  1:21511/1‑149 (MQ=255)                                                                      gTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCga                                                                            >  2:244974/1‑149 (MQ=255)                                                                      gTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑AAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCga                                                                            >  7:25825/1‑149 (MQ=255)                                                                         cAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCgatg                                                                          <  8:270318/148‑1 (MQ=255)                                                                         cAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCgatg                                                                          <  4:124286/148‑1 (MQ=255)                                                                           aCTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATc                                                                       <  8:116098/149‑1 (MQ=255)                                                                               cGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGc                                                                     <  1:244974/147‑1 (MQ=255)                                                                               cGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGAAACGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGcct                                                                   >  1:79742/1‑149 (MQ=255)                                                                                  cTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGcctcct                                                                >  1:190432/1‑149 (MQ=255)                                                                                      aTTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGt                                                            >  3:250493/1‑149 (MQ=255)                                                                                          ccTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGa                                                                                             >  2:173364/1‑112 (MQ=255)                                                                                          ccTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGa                                                                                             <  1:173364/112‑1 (MQ=255)                                                                                           cTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGGACAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGAt                                                         <  1:96253/148‑1 (MQ=255)                                                                                           cTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATc                                                                       <  2:234026/133‑1 (MQ=255)                                                                                           cTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATc                                                                       >  1:234026/1‑133 (MQ=255)                                                                                            tGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCt                                                      >  5:313723/1‑149 (MQ=255)                                                                                             gCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCgat                                                                           <  1:308467/127‑1 (MQ=255)                                                                                             gCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCgat                                                                           >  2:308467/1‑127 (MQ=255)                                                                                              cACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGGGCAGATCGAt                                                                                                                                         <  2:291918/66‑1 (MQ=255)                                                                                              cACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGGGCAGATCGAt                                                                                                                                         >  1:291918/1‑66 (MQ=255)                                                                                              cACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGGGCAGATCGAt                                                                                                                                         >  1:296708/1‑66 (MQ=255)                                                                                              cACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGGGCAGATCGAt                                                                                                                                         <  2:296708/66‑1 (MQ=255)                                                                                              cACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCt                                                      <  6:214100/147‑1 (MQ=255)                                                                                                    aaCCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAg                                                                               >  2:23779/1‑116 (MQ=255)                                                                                                    aaCCGGACGCGCCGGCGGCGATAACTTGACGCGGGCGGCATCCTGATCGAGCACATCGATTGCCTTGTCCGG‑CAGAAAGCGCGCGCTCGACGTAGCGGTCGGACAGTTCGGCGGCAg                                                                               <  1:23779/117‑1 (MQ=255)                                                                                                     aCCGGACGCGCCGTCGCCGATAACTTGAAGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTtgtg                                             <  5:13745/149‑1 (MQ=255)                                                                                                         gACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTgc                                         >  7:263649/1‑149 (MQ=255)                                                                                                          aCGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAg                                                                                                     <  5:321984/88‑1 (MQ=255)                                                                                                          aCGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAg                                                                                                     >  6:321984/1‑88 (MQ=255)                                                                                                            gcgcCGTCGCCGATAACTTGACGCTGGCGGCAGCCTGATGGTGCAGATCGATGGCCTTGTCCGG‑CAGATAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCt                                      <  8:282802/149‑1 (MQ=255)                                                                                                            gagcCGTCGCCGATAACTTGACGAGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGAT‑GACGTAGCGGTCGGATAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCt                                      >  2:21550/3‑149 (MQ=255)                                                                                                             cgcCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGcct                                                                   <  4:47528/119‑1 (MQ=255)                                                                                                             cgcCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGcct                                                                   >  3:47528/1‑119 (MQ=255)                                                                                                                  tCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGATAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAg                                <  1:21550/149‑1 (MQ=255)                                                                                                                       gttAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCg                                                              >  2:315047/3‑114 (MQ=255)                                                                                                                        aTAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTcg                          >  7:103161/1‑149 (MQ=255)                                                                                                                         tAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTcgc                         <  2:68673/149‑1 (MQ=255)                                                                                                                            cTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTCGCGCa                      <  8:25825/149‑1 (MQ=255)                                                                                                                             ttGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTCGCGCAg                     >  2:101730/1‑149 (MQ=255)                                                                                                                             ttGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTCGCGCAg                     >  1:239654/1‑149 (MQ=255)                                                                                                                               gACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGaa                                 <  1:223348/135‑1 (MQ=255)                                                                                                                               gACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGaa                                 >  2:223348/1‑135 (MQ=255)                                                                                                                                    gggcggCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTCGCGCAGGCcgcg               >  4:304298/1‑148 (MQ=255)                                                                                                                                    gggcggCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTCGCGCAGGCcgcg               >  3:180125/1‑148 (MQ=255)                                                                                                                                    gggcggCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTCGCGCAGGCcgcg               >  4:181103/1‑148 (MQ=255)                                                                                                                                        ggcAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTCGCGCAGGCCGCGCAGGa          <  4:180125/149‑1 (MQ=255)                                                                                                                                          cAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTCGCGCAGGCCGCGCAGGATc        <  2:21511/149‑1 (MQ=255)                                                                                                                                            gCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTCGCGCAGGCCGCGCAGGATCat      <  2:190432/149‑1 (MQ=255)                                                                                                                                             ccTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTCGCGCAGGCCGCGCAGGATCat      >  4:299583/1‑148 (MQ=255)                                                                                                                                             ccTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTCGCGCAGGCCGCGCAGGATCat      <  3:299583/148‑1 (MQ=255)                                                                                                                                                 aTCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTCGCGCAGGCCGCGCAGGATCatgatg  <  3:304298/148‑1 (MQ=255)                                                                                                                                                 aTCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTCGCGCAGGCCGCGCAGGATCatgat   <  8:263649/147‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CTTGCTGATCTGCGCCGCATTGTCATATTGCTTGCGCGAAGCCGCGTAGTCCTGCTCGCGGCGCAACTGGTGCAACTCCGACTCCATTTCCTGCACCGCAACCGGACGCGCCGTCGCCGATAACTTGACGCGGGCGGCAGCCTGATCGAGCAGATCGATGGCCTTGTCCGG‑CAGAAAGCGCGCGCT‑GACGTAGCGGTCGGACAGTTCGGCGGCAGCGATGATCGCCTCCTCGGTGATGCTGACCTTGTGGTGCGCTTCGAAGGTGTCGCGCAGGCCGCGCAGGATCATGATG  >  NC_002947/3752185‑3752476

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