Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,001,876	C→G	10.3%	A186G (GCC→GGC)	PP_2625 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,001,876	0	C	G	10.3%	197.1 / 13.9	69	A186G (GCC→GGC)	PP_2625	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/30);  new base G (5/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-01

CCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCGGTGGTAGCGCAGGTA  >  NC_002947/3001733‑3002023                                                                                                                                                |                                                                                                                                                    ccGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGa                                                                                                                                               >  1:1132521/1‑150 (MQ=255)   ggCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGc                                                                                                                                             <  1:9183/150‑1 (MQ=255)         aGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGa                                                                                                                                          <  1:475065/147‑1 (MQ=255)         aGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGa                                                                                                                                          >  2:475065/1‑147 (MQ=255)         aGCTAATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCa                                                                                                                                       >  2:766357/1‑150 (MQ=255)             gATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGaa                                                                                                                                   <  1:986122/150‑1 (MQ=255)               ttGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGc                                                                                                                                 <  1:538477/150‑1 (MQ=255)                  ccAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGCGGGCAGGCCCCGCCCGGCATCGAg                                                                                                                                              >  2:1018746/1‑134 (MQ=38)                  ccAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGGCGGGCAGGCCCCGCCCGGCATCGAg                                                                                                                                              <  1:1018746/134‑1 (MQ=38)                  ccAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGa                                                                                                                              >  2:4232/1‑150 (MQ=255)                    acgCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATc                                                                                                                            <  1:211464/148‑1 (MQ=255)                            gATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCa                                                                                                                    >  1:1281277/1‑150 (MQ=255)                              tCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGc                                                                                                                  <  2:1281277/150‑1 (MQ=255)                                 aGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGAAGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTa                                                                                                               >  1:1381357/1‑150 (MQ=255)                                  ggCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTAt                                                                                                              >  1:1377401/1‑150 (MQ=255)                                   gCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATc                                                                                                             >  2:501940/1‑150 (MQ=255)                                      aGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAg                                                                                                          <  2:885522/150‑1 (MQ=255)                                           ggTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTtgt                                                                                                     <  1:501940/150‑1 (MQ=255)                                           ggTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGAAGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTtgt                                                                                                     <  2:1381357/150‑1 (MQ=255)                                            gTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTtgtg                                                                                                    >  1:389585/1‑150 (MQ=255)                                              cAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTgtgag                                                                                                  <  2:603662/150‑1 (MQ=255)                                                  gTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTTGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAg                                                                                                          <  2:389344/138‑1 (MQ=255)                                                  gTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTTGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAg                                                                                                          >  1:389344/1‑138 (MQ=255)                                                   tACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATc                                                                                             <  2:1027333/150‑1 (MQ=255)                                                      gCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACg                                                                                          >  1:202243/1‑150 (MQ=255)                                                         tGAGCGGGGAGGGCTGGTCTACCCGTGAGACGACCTGGCAAGGGATGGAGAACGGCCAAGTGCTCGATGCAGGGCGGGGCGTGCCAGCCATCAAGCAGATCCGCAAGCTGATCGGTGACCAACGTATCGAGGTTGTGAGTGAACACAGcc                                                                                       >  2:29660/1‑149 (MQ=255)                                                           aGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAg                                                                                     <  1:506719/150‑1 (MQ=255)                                                              ggggAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTc                                                                                  <  2:202243/150‑1 (MQ=255)                                                                 gAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGAc                                                                               <  2:893065/150‑1 (MQ=255)                                                                     gCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGc                                                                           >  1:873326/1‑150 (MQ=255)                                                                            tACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCtt                                                                    <  1:766357/150‑1 (MQ=255)                                                                              cccGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCg                                                                  >  2:574723/1‑150 (MQ=255)                                                                                     gaCGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGAt                                                                <  2:1076616/145‑1 (MQ=255)                                                                                     gaCGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGAt                                                                >  1:1076616/1‑145 (MQ=255)                                                                                       cGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAg                                                         <  2:1309946/150‑1 (MQ=255)                                                                                        gCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGCTCCCGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAgg                                                        >  2:337093/1‑150 (MQ=255)                                                                                         cccTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGc                                                       >  1:808370/1‑150 (MQ=255)                                                                                          ccTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTtgtg                                                                                                    >  2:1417043/1‑104 (MQ=255)                                                                                          ccTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTtgtg                                                                                                    <  1:1417043/104‑1 (MQ=255)                                                                                          ccTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTAGCCGGGCTCACCTTCGATGCGCAAGGcc                                                      >  2:330132/1‑150 (MQ=255)                                                                                            tGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTg                                                    >  1:685841/1‑150 (MQ=255)                                                                                             ggCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTgg                                                   <  2:1132521/150‑1 (MQ=255)                                                                                              gCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGCTCTACGGGCTGACCTTCGATGCGCAAGGCCTGGc                                                  <  1:337093/150‑1 (MQ=255)                                                                                                  gTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGccccc                                              <  2:808370/150‑1 (MQ=255)                                                                                                     aTGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTg                                           >  2:407489/1‑150 (MQ=255)                                                                                                         agaACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCg                                                                                                                           >  2:17764/1‑66 (MQ=255)                                                                                                         agaACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCg                                                                                                                           <  1:17764/66‑1 (MQ=255)                                                                                                           aaCGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCa                                                                                            <  2:723220/95‑1 (MQ=255)                                                                                                           aaCGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCa                                                                                            >  1:723220/1‑95 (MQ=255)                                                                                                              cctgatctGCTCGATGGCGGGCAGGCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCgg                                  >  2:471911/8‑150 (MQ=37)                                                                                                               gCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCggg                                 >  1:105781/1‑150 (MQ=255)                                                                                                                ccAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGgc                                <  2:389585/150‑1 (MQ=255)                                                                                                                       cTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGgcg                               >  2:57804/1‑144 (MQ=255)                                                                                                                       cTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGgcg                               <  1:57804/144‑1 (MQ=255)                                                                                                                        tCGATGGCGGGCAGGCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCtt                                            <  1:801778/130‑1 (MQ=38)                                                                                                                        tCGATGGCGGGCAGGCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCtt                                            >  2:801778/1‑130 (MQ=38)                                                                                                                         cGATGGCGGGCAGGCCCCGCCCGACATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGACCAGCCGGGCGCGAAGcac                       >  2:171026/1‑150 (MQ=255)                                                                                                                          gATGGCGGGCAGGCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGcaca                      >  2:716718/1‑150 (MQ=255)                                                                                                                            tggcgggcaggCCCCGCCCGGCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCACGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATg                    >  2:336616/4‑150 (MQ=255)                                                                                                                             gCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGt                   <  1:574723/150‑1 (MQ=255)                                                                                                                                gggcgggGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCgg                >  2:861094/1‑150 (MQ=255)                                                                                                                                     gggCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCGGTGGTa           >  1:241631/1‑150 (MQ=255)                                                                                                                                      ggCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTc                                                                   <  2:516221/93‑1 (MQ=255)                                                                                                                                      ggCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTc                                                                   >  1:516221/1‑93 (MQ=255)                                                                                                                                         cTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCGGTGGTAgcgc       >  2:1227914/1‑150 (MQ=255)                                                                                                                                          tGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCGGTGGTAGCGCa      <  1:4232/150‑1 (MQ=255)                                                                                                                                             ccGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCgg                                  >  2:1369792/1‑119 (MQ=255)                                                                                                                                             ccGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCgg                                  <  1:1369792/119‑1 (MQ=255)                                                                                                                                              cGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCGGTGGTAGCGCAGGTa  >  2:1335/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                    CCGGCTACAGCTGATTGCCAGGCGTCCGATCGAGGCAAGCTGGGTCAAGGTACGCCTGAGCGGGGAGGGCTGGTCTACCCGTGAGACGCCCTGGCAAGTGATGGAGAACGGCCAAGTGCTCGATGCCGGGCGGGGCCTGCCCGCCATCGAGCAGATCAGGAAGCTGATCGGTGCCCAGCGTATCGAGGTTGTGAGTGATCACGGCGAGGTCGACGGGCTGACCTTCGATGCGCAAGGCCTGGCCCCCTTGATCAGCCGGGCGCGAAGCACATGTCGGTGGTAGCGCAGGTA  >  NC_002947/3001733‑3002023

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