Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,333,896	G→T	6.0%	R556S (CGC→AGC)	PP_2932 ←	amidase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,333,896	0	G	T	6.0%	213.4 / 3.7	67	R556S (CGC→AGC)	PP_2932	amidase family protein
Reads supporting (aligned to +/- strand):  ref base G (37/26);  new base T (2/2);  total (39/28)
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.81e-01

GGCGTTCTTCGCTGGGGGGGATAGGCATGGCAGACCGCAAAATCAACTGTTTGTGCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGCACGGTCACGGTCGGCACGCCCAGGTGACG  >  NC_002947/3333748‑3334040                                                                                                                                                     |                                                                                                                                                  ggCGTTCTTCGCTGGGGGGGATAGGCATGGCAGACCGCAAAATCAACTGTTTGTGCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGgcg                                                                                                                                                   >  8:131127/1‑149 (MQ=255)  gCGTTCTTCGCTGGGGGGGATAGGCATGGCAGACCGCAAAATCAACTGTTTGTGCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGgcgc                                                                                                                                                  <  5:334695/149‑1 (MQ=255)  gCGTTCTTCGCTGGGGGGGATAGGCATGGCAGACCGCAAAATCAACTGTTTGTGCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGgcgc                                                                                                                                                  >  4:174086/1‑149 (MQ=255)        ttcGCTGGGGGGGATAGGCATGGCAGACCGCAAAATCAACTGTTTGTGCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCg                                                                                                                                            >  7:24432/1‑149 (MQ=255)        ttcGCTGGGGGGGATAGGCATGGCAGACCGCAAAATCAACTGTTTGTGCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCg                                                                                                                                            >  3:74391/1‑149 (MQ=255)           gCTGGGGGGGATAGGCATGGCAGACCGCAAAATCAACTGTTTGTGCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTg                                                                                                                                         <  7:131127/149‑1 (MQ=255)                    gATAGGCATGGCAGACCGCAAAATCAACTGTTTGTGCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCTCGCGCCggatgg                                                                                                                                       <  6:105164/142‑5 (MQ=255)                    gATAGGCATGGCAGACCGCAAAATCAACTGTTTGTGCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCTCGCGCCggatgg                                                                                                                                       >  5:105164/1‑138 (MQ=255)                        ggCATGGCAGACCGCAAAATCAACTGTTTGTGCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCgg                                                                                                                            <  5:198276/149‑1 (MQ=255)                              gCAGACCGCAAAATCAACTGTTTGTGCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGgcagca                                                                                                                       >  2:192449/1‑148 (MQ=255)                                      cAAAATCAACTGTTTGTGCATTTTTACACATATAGTTCAAATTCCCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGAGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGACTCGAACGCGGCAGCAAAACTcag                                                                                                               >  3:190256/1‑148 (MQ=255)                                         aaTCAACTGTTTGTGCATTTTTACACATATAGTTGAGATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGCGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAGCGCGGCAGCAGAGCTCagcagc                                                                                                           >  6:275444/1‑149 (MQ=255)                                          aTCAACTGTTTGTGCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAgcg                                                                                                          >  1:266245/1‑149 (MQ=255)                                             aaCTGTTTGTGCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCTCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGt                                                                                                       >  2:36232/1‑149 (MQ=255)                                              aCTGTTTGTGCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGCCTCAGCCGAGCGCCGGCGTGCGCGGCGGAGCCCTCCGGCGCGAGCCGGTGGCCGCGAACGGGGCGGCAATCCTCAGGAGCGGGtt                                                                                                      >  4:265657/1‑149 (MQ=255)                                                  tttGTGCATTTTTACACAGATAGTTGAAAGTCGCGATCTGCTGTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGGGGCCTCGACCGCGGCGGCACCTCTCAGCAGCGGGGTGTCg                                                                                                  >  4:70302/1‑149 (MQ=255)                                                     gtgCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCta                                                                                               >  4:109758/1‑149 (MQ=255)                                                        cATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCTCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAg                                                                                            <  1:36232/149‑1 (MQ=255)                                                          tttttACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGc                                                                                          >  3:109652/1‑149 (MQ=255)                                                                cacaTATAGTTGAAATTCGCGATCTGCTATCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTccc                                                                                    >  5:116824/1‑149 (MQ=255)                                                                   atatAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTc                                                                                                                                    <  7:187286/98‑1 (MQ=255)                                                                   atatAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCgg                                                                                  >  7:99551/1‑148 (MQ=255)                                                                      taGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCaaa                                                                              <  1:192449/149‑1 (MQ=255)                                                                        gTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAgg                                                                            >  4:133493/1‑149 (MQ=255)                                                                        gTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAg                                                                             <  8:99551/148‑1 (MQ=255)                                                                         ttGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCTCGCGCCggat                                                                                                                                         <  6:75529/87‑3 (MQ=255)                                                                         ttGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCTCGCGCCggat                                                                                                                                         >  5:75529/1‑85 (MQ=255)                                                                         ttGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAgg                                                                            >  7:325964/1‑148 (MQ=255)                                                                         ttGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGt                                                                           >  5:104792/1‑149 (MQ=255)                                                                          tGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAAACATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGt                                                                           >  8:325464/1‑148 (MQ=255)                                                                             aaTTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGc                                                                       <  3:174086/149‑1 (MQ=255)                                                                               ttCGCGATCTGCTTTCGCGGGCCCGGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAGCCCTCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAGAGTCAGCCGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGccc                                                                     >  2:191535/1‑149 (MQ=255)                                                                                   cgATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCAcc                                                                 >  4:39322/1‑149 (MQ=255)                                                                                     aTCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCgg                                                               >  5:365153/1‑149 (MQ=255)                                                                                          cTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGc                                                              <  7:325464/145‑1 (MQ=255)                                                                                          cTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATAc                                                          <  4:74391/149‑1 (MQ=255)                                                                                          cTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATAc                                                          <  4:190256/149‑1 (MQ=255)                                                                                             tCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAATGTCAGCCCCACCGGCATAcc                                                         <  3:109758/147‑1 (MQ=255)                                                                                             tCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATAcc                                                         <  8:24432/147‑1 (MQ=255)                                                                                             tCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCg                                                        <  3:133493/148‑1 (MQ=255)                                                                                               gcgGGCCCTGTAAGACTCAGCCTAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGCGCCGGTGGCCTCTAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGTCAAAGGTAAGCCCCACCGGCATACCGata                                                     >  7:54103/1‑149 (MQ=255)                                                                                               gcgGGCCCTGTAAGACTCAGCAGAGCGCCGGCGTGCGCGGCGGAACCATACGGCGCGACACGGTGGCCTCGAACGCGGCAGCAAAAATCAGCAGCGCGTTGTAGATATAGGCCTTCCCGGCAAAGGTCAGCCCCACAGGCATACCGata                                                     >  7:52399/1‑149 (MQ=255)                                                                                                      cTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCg                                                                                                  <  8:173829/97‑1 (MQ=255)                                                                                                      cTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCg                                                                                                  >  7:173829/1‑97 (MQ=255)                                                                                                       tGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGc                                                                       >  2:255030/1‑123 (MQ=255)                                                                                                        gTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCg                                                                <  6:244627/129‑1 (MQ=255)                                                                                                        gTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCg                                                                >  5:244627/1‑129 (MQ=255)                                                                                                        gTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATc                                           >  6:344104/1‑149 (MQ=255)                                                                                                         tAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCa                                          >  7:249507/1‑149 (MQ=255)                                                                                                         tAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCa                                          >  8:102218/1‑149 (MQ=255)                                                                                                            gACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGTCAAAGGGCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGc                                       <  5:275444/149‑1 (MQ=255)                                                                                                               tCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCa                                    >  2:224479/1‑149 (MQ=255)                                                                                                               tCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCa                                    >  4:45287/1‑149 (MQ=255)                                                                                                                     gAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGca                              <  3:45287/149‑1 (MQ=255)                                                                                                                      aGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGcac                             <  4:109652/149‑1 (MQ=255)                                                                                                                           cGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCTGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGcacggtca                        <  2:266245/149‑1 (MQ=255)                                                                                                                               gTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGcacggtcacggt                    >  5:337397/1‑149 (MQ=255)                                                                                                                                tGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATACTCGGCCATCACGCCCATCGGcacggtcacggt                    <  1:224479/149‑1 (MQ=255)                                                                                                                                  cgcgGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGCACGGTCACGGTCgg                 <  8:325964/149‑1 (MQ=255)                                                                                                                                  cgcgGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGCACGGTCACGGTCgg                 >  8:175425/1‑149 (MQ=255)                                                                                                                                        ggAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGCACGGTCACGGTCGGCACGcc           <  5:344104/149‑1 (MQ=255)                                                                                                                                            ccATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGCACGGTCACTGTCGGCACGCCCatg       <  8:249507/149‑3 (MQ=255)                                                                                                                                            ccATCCGGCGCGAGCCGCTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTACCGGCAAAGGTCAGCCAACCCGGCATACCGATA‑TCGGCCAACCCGCCCCTCGGCACGGTCACGGTCGGCCCGCCCAgg       >  8:185003/1‑149 (MQ=255)                                                                                                                                             cATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGCACGGTCACGGTCGGCACGCCCAGGt      <  7:102218/149‑1 (MQ=255)                                                                                                                                              aTCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGCACGGTCACGGTCGGCACGCCCAGGt      >  2:16358/1‑148 (MQ=255)                                                                                                                                               tCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGCACGGTCACGGTCGGCACGc            <  2:188991/141‑1 (MQ=255)                                                                                                                                               tCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGCACGGTCACGGTCGGCACGc            >  1:188991/1‑141 (MQ=255)                                                                                                                                               tCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGCACGGTCACGGTCGGCACGCCCAGGTGa    >  7:11303/1‑149 (MQ=255)                                                                                                                                                 cGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCCGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGCACGGTCACGGTCGGCACGCCCAGGTGACg  <  7:185003/149‑1 (MQ=255)                                                                                                                                                  ggCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGCACGGTCACGGTCGGCACGCCCAGGTGACg  >  4:10319/1‑148 (MQ=255)                                                                                                                                                     gcgAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGat                                                      <  3:163910/94‑1 (MQ=255)                                                                                                                                                     gcgAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGat                                                      >  4:163910/1‑94 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GGCGTTCTTCGCTGGGGGGGATAGGCATGGCAGACCGCAAAATCAACTGTTTGTGCATTTTTACACATATAGTTGAAATTCGCGATCTGCTTTCGCGGGCCCTGTAAGACTCAGCCGAGCGCCGGCGTGCGCGGCGGAACCATCCGGCGCGAGCCGGTGGCCTCGAACGCGGCAGCAAAACTCAGCAGCGCGTTGTCGCTATAGGCCTTCCCGGCAAAGGTCAGCCCCACCGGCATACCGATA‑TCGGCCATCACGCCCATCGGCACGGTCACGGTCGGCACGCCCAGGTGACG  >  NC_002947/3333748‑3334040

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