Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,475,136	C→G	8.8%	P148P (CCG→CCC)	PP_3964 ←	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,475,136	0	C	G	8.8%	200.6 / 8.2	68	P148P (CCG→CCC)	PP_3964	transposase
Reads supporting (aligned to +/- strand):  ref base C (28/34);  new base G (3/3);  total (31/37)
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.75e-01

AACTTGGCCACCATCACGTGGGCCAACAAACCTGTGGTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGACGAATCACCGTGCGTGGAAACTGCGGCGGCAAAGGCGCACGTTTCGGGGATTGTCGTG  >  NC_002947/4474990‑4475281                                                                                                                                                   |                                                                                                                                                  aaCTTGGCCACCATCACGTGGGCCAACAAACCTGTGGTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGc                                                                                                                                                <  2:32332/150‑1 (MQ=255)   cTTGGCCACCATCACGTGGGCCAACAAACCTGTGGTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCgt                                                                                                                                              <  2:756064/150‑1 (MQ=255)   cTTGGCCACCATCACGTGGGCCAACAAACCTGTGGTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCgt                                                                                                                                              <  2:89834/150‑1 (MQ=255)     tGGCCACCATCACGTGCGCCCACAACCCTGTTGTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCgtgt                                                                                                                                            <  1:593661/150‑1 (MQ=255)       gCCACCATCACGTGGGCCAACAAACCTGTGGTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTaa                                                                                                                                          <  1:492136/150‑1 (MQ=255)       gCCACCATCACGTGGGCCAACAAACCTGTGGTCGGGATACCTTTGTCGATAACCTGGACTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTaa                                                                                                                                          >  2:94544/1‑150 (MQ=255)         caccaTCACGTGGGCCAACAAACCTGTGGTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATc                                                                                                                                        <  1:76593/150‑1 (MQ=255)          accaTCACGTGGGCCAACAAACCTGTGGTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCGGGCgtgttta                                                                                                                                         <  2:948788/148‑4 (MQ=255)          accaTCACGTGGGCCAACAAACCTGTGGTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCGGGCgtgttta                                                                                                                                         >  1:948788/1‑145 (MQ=255)             aTCACGTGGGCCAACAAACCTGTGGTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCACGACGGTGAACACTCCCGGCGTGTAATCCAGc                                                                                                                                    >  1:1161426/1‑150 (MQ=255)                  gTGGGCCAACAAACCTGTGGTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCGGGCgtgttt                                                                                                                                          <  2:140637/139‑3 (MQ=255)                  gTGGGCCAACAAACCTGTGGTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCGGGCgtgttt                                                                                                                                          >  1:140637/1‑137 (MQ=255)                          acaaACCTGTGGTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGt                                                                                                                       >  1:742511/1‑150 (MQ=255)                           taatccTGTGTTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTc                                                                                                                      <  2:1222890/146‑1 (MQ=255)                                     gTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATg                                                                                                            >  1:317275/1‑150 (MQ=255)                                     gTCGGGATACCTTTGTCGATAACCTGGACTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATg                                                                                                            <  1:94544/150‑1 (MQ=255)                                          gATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGAc                                                                                                                         >  2:687788/1‑132 (MQ=255)                                          gATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGAc                                                                                                                         <  1:687788/132‑1 (MQ=255)                                          gATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTg                                                                                                       >  1:511891/1‑150 (MQ=255)                                            tACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGaa                                                                                                     <  2:968164/150‑1 (MQ=255)                                               cTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGtt                                                                                                  <  1:415587/150‑1 (MQ=255)                                                tttGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGATGTCTTCGCCGATGCGTTGAAGTTg                                                                                                 >  2:835604/1‑150 (MQ=255)                                                 ttGTCGATAACCTGGGCTGGCACCGGGGTCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTgg                                                                                                <  2:317275/150‑1 (MQ=255)                                                       aTAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCca                                                                                          >  2:1148259/1‑150 (MQ=255)                                                           ccTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGa                                                                                      <  2:742511/150‑1 (MQ=255)                                                            cTGGGCTGGCACCGGGGCCTGGATCAGGGTTTGGCACTGGGGACAGGCCCATGTGCCCGGCACATGTTGCTCGAAGGTAAACACGCCCGGCGTGTAATCAAGCTTGGAGCTGACGTCTTCGCCGATGAGTTGACGTTGGCAGCCACAGAc                                                                                     >  1:1161829/1‑150 (MQ=14)                                                            cTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTg                                                                                                 <  2:1158521/138‑1 (MQ=255)                                                            cTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTg                                                                                                 >  1:1158521/1‑138 (MQ=255)                                                            cTGGGCTGGAAGCGAGGCCGGGAGCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGAc                                                                                     >  2:1156640/1‑150 (MQ=14)                                                                    gCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTggg                                                                             <  2:1400850/150‑1 (MQ=255)                                                                     cACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGgt                                                                            <  1:76887/150‑1 (MQ=255)                                                                     cACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGgt                                                                            <  1:162000/150‑1 (MQ=255)                                                                      aCCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGtt                                                                                                  >  1:92634/1‑127 (MQ=255)                                                                      aCCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGtt                                                                                                  <  2:92634/127‑1 (MQ=255)                                                                      aCCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACATACGCACTGGgtg                                                                           >  1:693240/1‑150 (MQ=255)                                                                      aCCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGgtg                                                                           >  1:1177109/1‑150 (MQ=255)                                                                        cGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGtt                                                                         >  2:369571/1‑150 (MQ=255)                                                                          gggCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGtttt                                                                       >  2:106464/1‑150 (MQ=255)                                                                           ggCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTc                                                                      <  1:106819/150‑1 (MQ=255)                                                                             ccTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTgg                                                                              <  1:15534/140‑1 (MQ=255)                                                                             ccTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTgg                                                                              >  2:15534/1‑140 (MQ=255)                                                                             ccTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCgg                                                                    <  1:1148259/150‑1 (MQ=255)                                                                               tGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGTTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCgggg                                                                  <  2:1161829/150‑2 (MQ=255)                                                                                     aGGGTTTCGGACTGGCCACAGGCCCATTTGCCCCGCACATGTTGTTTGACGGTAAACACGCCCGGCGTGTAATCCAGCTTATCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGa                                                            <  1:82601/150‑1 (MQ=14)                                                                                     aGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGa                                                            <  1:369571/150‑1 (MQ=255)                                                                                            cGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGACGAATCa                                                     <  1:106464/150‑1 (MQ=255)                                                                                             gCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGACGAATCAc                                                    <  2:1177109/150‑1 (MQ=255)                                                                                             gCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGACGAATCAc                                                    >  1:129762/1‑150 (MQ=255)                                                                                                      aCAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGACGACTCACCGTGCGTgg                                           >  1:218776/1‑150 (MQ=255)                                                                                                           cccATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCca                                                                                          <  1:394307/98‑1 (MQ=255)                                                                                                           cccATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCca                                                                                          >  2:394307/1‑98 (MQ=255)                                                                                                           cccATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTggg                                                                             >  1:522946/1‑111 (MQ=255)                                                                                                           cccATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTggg                                                                             <  2:522946/111‑1 (MQ=255)                                                                                                                   ccccGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGa                                                                                                              <  2:1382331/70‑1 (MQ=40)                                                                                                                   ccccGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGa                                                                                                              >  1:1382331/1‑70 (MQ=40)                                                                                                                   ccccGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATTCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGACGAATCACCGTGCGTGGAAACTGCggcggc                              >  2:986745/1‑150 (MQ=255)                                                                                                                    cccGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGAGGATGGTGGCCGACCAAGACGGCCACCAGAACGATGGGTTCGTGACGAATCACCGTGCGTGGAAACTGCGGCGGCa                             >  1:674216/1‑150 (MQ=9)                                                                                                                     ccGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCg                                                                                                          <  2:553916/72‑1 (MQ=38)                                                                                                                     ccGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCg                                                                                                          >  1:553916/1‑72 (MQ=38)                                                                                                                                    gACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGACGAATCACCGTGCGTGGAAACTGCGGCGGCAAAGGCGCACGTTTCgg             >  2:1093389/1‑150 (MQ=255)                                                                                                                                    gACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGACGAATCACCGTGCGTGGAAACTGCGGCGGCAAAGGCGCACGTTTCgg             <  1:923571/150‑1 (MQ=255)                                                                                                                                         tAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACATACGCACTGGGTGTTTTCGGGTTCGTGACGAATCACCGTGCGTGGAAACTGCGGCGGCAAAGGCGCACGTTTCGGGGAtt        <  2:693240/150‑1 (MQ=255)                                                                                                                                         tAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGACGAATCACCGTGCGTGGAAACTGCGGCGGCAAAGGCGCACGTTTCGGGGAtt        >  2:3528/1‑150 (MQ=255)                                                                                                                                         attacacGCCGGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGACGAATCACCGTGCGTGGAAACTGCGGCGGCAAAgg                         >  2:1029480/4‑133 (MQ=255)                                                                                                                                         attacacGCCGGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGACGAATCACCGTGCGTGGAAACTGCGGCGGCAAAgg                         <  1:1029480/130‑1 (MQ=255)                                                                                                                                           aaCACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGACGAATCACCGTGCGTGGAAACTGCGGCGGCAAAGGCGCACGTTTCGGGGATTGt      <  2:129762/150‑1 (MQ=255)                                                                                                                                            acacGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGACGAATCACCGTGCGTGGAAACTGCGGCGGCAAAGGCGCACGTTTCGGGGATTGTc     >  2:31108/1‑150 (MQ=255)                                                                                                                                               cGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGACGAATCACCGTGCGTGGAAACTGCGGCGGCAAAGGCGCACGTTTCGGGGATTGTCGTg  <  2:218776/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  AACTTGGCCACCATCACGTGGGCCAACAAACCTGTGGTCGGGATACCTTTGTCGATAACCTGGGCTGGCACCGGGGCCTGGATCAGGGTTTCGCACTGGCGACAGGCCCATTTGCCCCGCACATGTTGCTCGACGGTAAACACGCCCGGCGTGTAATCCAGCTTCTCGCTGACGTCTTCGCCGATGCGTTGAAGTTGGCAGCCACAGACGCACTGGGTGTTTTCGGGTTCGTGACGAATCACCGTGCGTGGAAACTGCGGCGGCAAAGGCGCACGTTTCGGGGATTGTCGTG  >  NC_002947/4474990‑4475281

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