Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	861,659	C→A	6.8%	G36C (GGC→TGC)	PP_0742 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	861,659	0	C	A	6.8%	236.0 / 4.1	73	G36C (GGC→TGC)	PP_0742	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/37);  new base A (2/3);  total (33/40)
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.82e-01

TTGGGGTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAGCAGACTACGCTGTCAGG  >  NC_002947/861513‑861808                                                                                                                                                     |                                                                                                                                                      ttGGGGTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTg                                                                                                                                                    >  2:495303/1‑150 (MQ=255) ttGGGGTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTg                                                                                                                                                    >  1:479236/1‑150 (MQ=255) ttGGGGTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTg                                                                                                                                                    >  1:375323/1‑150 (MQ=255)   ggggTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGc                                                                                                                                                  <  2:1116068/150‑1 (MQ=255)   ggggTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGCGCCGTACTGTGCGCAACCAGCAACTG‑GCCCTGGC‑TCCCTACGCTTGCGGGCCATGCGCCCGCCACTGCCGTgga                                                                                                                                                  <  1:328571/150‑2 (MQ=255)      gTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGCAGACGCGTCGGGCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCtt                                                                                                                                               >  2:595241/1‑150 (MQ=255)      gTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCtt                                                                                                                                               >  1:1032222/1‑150 (MQ=255)      gTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACGGCAGTGGcggg                                                                                                                                               <  1:184571/150‑4 (MQ=255)          ccACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTc                                                                                                                                           <  2:479236/150‑1 (MQ=255)            aCTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGa                                                                                                                                         >  1:900256/1‑150 (MQ=255)            aCTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGa                                                                                                                                         >  2:150508/1‑150 (MQ=255)             cTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGaa                                                                                                                                        <  1:1067327/150‑1 (MQ=255)               gCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGTCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACa                                                                                                                                      >  1:46420/1‑150 (MQ=255)                cTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGATGAAACGTCGGGCACGCGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACaa                                                                                                                                     <  2:264346/150‑1 (MQ=255)                cTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCCCGGGTAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACaa                                                                                                                                     <  1:617054/150‑1 (MQ=255)                cTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACaa                                                                                                                                     >  2:403554/1‑150 (MQ=255)                     cccAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCa                                                                                                                                <  1:188583/150‑1 (MQ=255)                                 ccACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGCATCAGTGCACGGATTGCCCgg                                                                                                                    >  2:795651/1‑150 (MQ=255)                                   aCCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGcccaggc                                                                                                                  >  1:982573/1‑146 (MQ=255)                                     ccATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCtt                                                                                                                <  1:495303/150‑1 (MQ=255)                                      cATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTc                                                                                                               >  1:256847/1‑150 (MQ=255)                                       aTTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCAGGcggact                                                                                                              <  2:901892/150‑6 (MQ=255)                                        ttGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGAGCCGGCGTCTCGGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTg                                                                                                             <  2:1002562/150‑1 (MQ=255)                                            tgcACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAg                                                                                                         >  2:84237/1‑150 (MQ=255)                                             gcACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAgg                                                                                                        >  1:100455/1‑150 (MQ=255)                                              cACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGc                                                                                                       <  1:84237/150‑1 (MQ=255)                                                    aTCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTg                                                                                                                                             >  2:391872/1‑106 (MQ=255)                                                    aTCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTg                                                                                                                                             <  1:391872/106‑1 (MQ=255)                                                       gacgaAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGggcggc                                                                                              <  2:785369/150‑1 (MQ=255)                                                        acgaAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGc                                                                                                                                                  <  1:322806/97‑1 (MQ=255)                                                        acgaAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGc                                                                                                                                                  >  2:322806/1‑97 (MQ=255)                                                          gaAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGa                                                                                           <  1:1171955/150‑1 (MQ=255)                                                               gTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGCGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAAc                                                                                      >  2:696950/1‑150 (MQ=255)                                                               gTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAAc                                                                                      >  2:749634/1‑150 (MQ=255)                                                                   gTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGc                                                                                  <  2:256847/150‑1 (MQ=255)                                                                         ggCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGcc                                                                            <  1:998694/150‑1 (MQ=255)                                                                                aaGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCg                                                                     >  1:821106/1‑150 (MQ=255)                                                                                      tgtgCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAAc                                                                                      <  1:637710/127‑1 (MQ=255)                                                                                      tgtgCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAAc                                                                                      >  2:637710/1‑127 (MQ=255)                                                                                        tgCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGc                                                             <  1:150508/150‑1 (MQ=255)                                                                                                cAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAgg                                                       <  2:731704/150‑1 (MQ=255)                                                                                                aaGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAgg                                                       >  1:1191147/2‑150 (MQ=255)                                                                                                  ggTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTa                                                     <  1:403554/150‑1 (MQ=255)                                                                                                         cgcgcCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCat                                              <  2:1032222/150‑1 (MQ=255)                                                                                                            gcCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCatgat                                           >  2:628822/1‑150 (MQ=255)                                                                                                            gcCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCatgat                                           >  2:88289/1‑150 (MQ=255)                                                                                                             ccGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCTGATGATGGCGATAGGTACTGTCatgatg                                          <  1:961268/150‑1 (MQ=255)                                                                                                              cGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATgg                                         >  1:874449/1‑150 (MQ=255)                                                                                                               ggCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGc                                        <  1:1020107/150‑1 (MQ=255)                                                                                                               ggCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCAATTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGc                                        <  1:434625/150‑1 (MQ=255)                                                                                                                gCTTCGCTGCGCTTGCTGGCCATGCGCCCGCAACTACCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGcc                                       >  1:157650/1‑150 (MQ=255)                                                                                                                  ttCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCTTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCctc                                     >  1:173583/1‑150 (MQ=255)                                                                                                                     gCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCga                                                                    <  1:767087/116‑1 (MQ=255)                                                                                                                     gCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCga                                                                    >  2:767087/1‑116 (MQ=255)                                                                                                                        gcgcTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCaaa                               <  1:628822/150‑1 (MQ=255)                                                                                                                         cgcTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGTATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAAc                              <  2:900256/150‑1 (MQ=255)                                                                                                                         cgcTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAAc                              <  1:595241/150‑1 (MQ=255)                                                                                                                                 ggCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCtgt                      <  2:1006402/150‑1 (MQ=255)                                                                                                                                      tGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAGCa                 <  2:1191147/150‑1 (MQ=255)                                                                                                                                      tGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAGCa                 <  1:130195/150‑1 (MQ=255)                                                                                                                                       acaaaggccACGGCAGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAg                   >  1:88359/7‑147 (MQ=255)                                                                                                                                       acaaaggccACGGCAGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAg                   <  2:88359/141‑1 (MQ=255)                                                                                                                                        cgcccgccACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCTTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAGCAGa               <  2:173583/150‑1 (MQ=255)                                                                                                                                         aaaggccACGGCAGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCtt                                   <  2:466749/125‑1 (MQ=255)                                                                                                                                         aaaggccACGGCAGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCtt                                   >  1:466749/5‑129 (MQ=255)                                                                                                                                           ccgccACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGc                                                                         <  2:571982/89‑1 (MQ=255)                                                                                                                                           ccgccACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGc                                                                         >  1:571982/1‑89 (MQ=255)                                                                                                                                           ccgccACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATgg                                                              >  2:314209/1‑100 (MQ=255)                                                                                                                                           ccgccACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATgg                                                              <  1:314209/100‑1 (MQ=255)                                                                                                                                            cgccACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAGCAGACTAc           >  2:232090/1‑150 (MQ=255)                                                                                                                                             gccACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAGCAGACTACg          >  1:675453/1‑150 (MQ=255)                                                                                                                                             gccACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAGCAGACTACg          <  2:100455/150‑1 (MQ=255)                                                                                                                                               cACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAGCAGACTACGCt        <  2:675453/150‑1 (MQ=255)                                                                                                                                                     cGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAGCAGACTACGCTGTCggg  <  2:908969/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TTGGGGTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAGCAGACTACGCTGTCAGG  >  NC_002947/861513‑861808

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