Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	861,656	T→G	6.8%	S37R (AGT→CGT)	PP_0742 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	861,656	0	T	G	6.8%	210.9 / 7.3	74	S37R (AGT→CGT)	PP_0742	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/37);  new base G (2/3);  total (34/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 = 4.57e-01

GTACAGTTGGGGTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAGCAGACTACGCT  >  NC_002947/861507‑861802                                                                                                                                                        |                                                                                                                                                   gTACAGTTGGGGTTTCCACTGCTCGGCCCAGCAGGCGGCCGCCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACt                                                                                                                                                    <  2:1151810/150‑1 (MQ=255)  tACAGTTGGGGTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTg                                                                                                                                                   >  2:188583/1‑150 (MQ=255)       ttGGGGTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTg                                                                                                                                              >  1:375323/1‑150 (MQ=255)       ttGGGGTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTg                                                                                                                                              >  2:495303/1‑150 (MQ=255)       ttGGGGTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTg                                                                                                                                              >  1:479236/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                                                                                                                                            >  2:322806/1‑97 (MQ=255)                                                              acgaAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGc                                                                                                                                            <  1:322806/97‑1 (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                                                                                >  2:637710/1‑127 (MQ=255)                                                                                            tgtgCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAAc                                                                                <  1:637710/127‑1 (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                                                              >  2:767087/1‑116 (MQ=255)                                                                                                                           gCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCga                                                              <  1:767087/116‑1 (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           <  1:130195/150‑1 (MQ=255)                                                                                                                                            tGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAGCa           <  2:1191147/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    <  2:100455/150‑1 (MQ=255)                                                                                                                                                   gccACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAGCAGACTACg    >  1:675453/1‑150 (MQ=255)                                                                                                                                                     cACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAGCAGACTACGCt  <  2:675453/150‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                   GTACAGTTGGGGTTTCCACTGCTCGGCCCAGCCGGCGGCCACCCATTGCTGCACCTGATCGACGAAGCGTCGGTCACGGGCAGTGAAGTACTGTGCGC‑‑CCAGGTCGAGCGCGCCGGCTTCGCTGCGCTTGCTGGCCATGCGCCCGCCACTGCCGTGGCCTTTGTCGAACAAGTGCACGGATTGCCCGGCCTTCTGCAAGGCCTGGGCGGCTGACAAACCGGCGATGCCGGCCCCGATGATGGCGATAGGTACTGTCATGATGGCCTCTTCAAACCTTGCTGTGAGCAGACTACGCT  >  NC_002947/861507‑861802

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: ‑