Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,332,268	T→C	6.5%	intergenic (+3/‑46)	dksA → / → gluQ	RNA polymerase‑binding transcription factor DksA/glutamyl‑Q tRNA(Asp) synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,332,268	0	T	C	6.5%	177.7 / 2.6	61	intergenic (+3/‑46)	dksA/gluQ	RNA polymerase‑binding transcription factor DksA/glutamyl‑Q tRNA(Asp) synthetase
Reads supporting (aligned to +/- strand):  ref base T (34/23);  new base C (2/2);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.96e-01

GATTCAGGACGAAGAGTACGGCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTC  >  NC_002947/5332121‑5332414                                                                                                                                                     |                                                                                                                                                   gaTTCAGGACGAAGAGTACGGCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGAtcttt                                                                                                                                                  <  1:1963441/150‑1 (MQ=255)   ttCAGGATGAAGAGTACGGCTGGTGCGAGTCGGGCGGCATCGAGATCGGCCTGCGCCGCCTGGAGGCCCGCCCTACCGCCGACCTGTGCTGTGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGAtctttct                                                                                                                                                <  1:1148000/150‑1 (MQ=255)    tCAGGACGAAGAGTACGGCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGAtctttctt                                                                                                                                               <  1:107058/150‑1 (MQ=255)       ggACGAAGAGTACGGCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCa                                                                                                                                            >  1:1650317/1‑150 (MQ=255)       ggACGAAGAGTACGGCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCCCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCAAAAAAAAAA‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCa                                                                                                                                            <  2:2121722/150‑1 (MQ=255)         aCGAAGAGTACGGCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCAcc                                                                                                                                          >  2:394599/1‑150 (MQ=255)           gAAGAGTACGGCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGAtctttc                                                                                                                                                 >  1:1821343/1‑141 (MQ=255)           gAAGAGTACGGCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGAtctttc                                                                                                                                                 <  2:1821343/141‑1 (MQ=255)           gAAGAGTACGGCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGAtcttt                                                                                                                                                  <  2:1935363/140‑1 (MQ=255)           gAAGAGTACGGCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGAtcttt                                                                                                                                                  >  1:1935363/1‑140 (MQ=255)           gAAGAGTACGGCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTg                                                                                                                                        >  2:1792509/1‑150 (MQ=255)           gAAGAGTACGGCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTg                                                                                                                                        >  1:1712142/1‑150 (MQ=255)              gagTACGGCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGAtcttt                                                                                                                                                  <  1:1312556/137‑1 (MQ=255)              gagTACAGCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGAtcttt                                                                                                                                                  >  2:1312556/1‑137 (MQ=255)                    ggCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGgcg                                                                                                                               >  2:1199635/1‑150 (MQ=255)                      cTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGca                                                                                                                             <  2:1215339/150‑1 (MQ=255)                            cGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATGATgc                                                                                                                       <  1:554215/150‑1 (MQ=255)                                  gTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCAACTGAAAGGAGAATGCCGAGAAAAAGG‑AAAATTCGGGCGACAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCg                                                                                                                 >  2:273738/1‑150 (MQ=255)                                       gCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCa                                                                                                            <  1:303172/150‑1 (MQ=255)                                         aTCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAtt                                                                                                          >  1:861972/1‑150 (MQ=255)                                          tCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAttt                                                                                                         >  1:792683/1‑150 (MQ=255)                                          tCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAttt                                                                                                         <  1:142897/150‑1 (MQ=255)                                             agaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTtata                                                                                                      <  1:244332/150‑1 (MQ=255)                                              gaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGTAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGcccc                                                                                                                  >  2:1418811/1‑137 (MQ=255)                                              gaTCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGcccc                                                                                                                  <  1:1418811/137‑1 (MQ=255)                                                   gCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGTGCGTGGCGCATCATGCGCCCCGGTTCATTTATATGCCCa                                                                                                >  2:1330033/1‑150 (MQ=255)                                                   gCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTta                                                                                                        >  2:176272/1‑142 (MQ=255)                                                   gCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTta                                                                                                        <  1:176272/142‑1 (MQ=255)                                                   gCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCa                                                                                                <  1:1532026/150‑1 (MQ=255)                                                    ccTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGAtcttt                                                                                                                                                  <  1:2066961/99‑1 (MQ=255)                                                    ccTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGAtcttt                                                                                                                                                  >  2:2066961/1‑99 (MQ=255)                                                        cgccgcCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGca                                                                                                                             >  2:1788637/1‑116 (MQ=255)                                                        cgccgcCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGca                                                                                                                             <  1:1788637/116‑1 (MQ=255)                                                           cgccTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGAcc                                                                                        <  1:643661/150‑1 (MQ=255)                                                             ccTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCAt                                                                                      <  2:792683/150‑1 (MQ=255)                                                               tGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGAtctttct                                                                                                                                                >  1:2320121/1‑90 (MQ=255)                                                               tGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGAtctttct                                                                                                                                                <  2:2320121/90‑1 (MQ=255)                                                                 gAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGAcc                                                                                  <  2:964404/150‑1 (MQ=255)                                                                    gcccgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATgc                                                                                                                       >  1:271988/1‑110 (MQ=255)                                                                    ccccgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATgc                                                                                                                       <  2:271988/109‑1 (MQ=255)                                                                     cccgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGACGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACt                                                                              >  1:1610255/1‑150 (MQ=255)                                                                      ccgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTc                                                                             >  2:916696/1‑150 (MQ=255)                                                                       cgcccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCCTTTATATGCCCAGCCTGACCATGACCGACGcc                                                                            >  1:509798/1‑150 (MQ=255)                                                                         cccTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTCATGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAg                                                                          <  1:1330033/150‑1 (MQ=255)                                                                          ccgaCCTCCGACCTGTCCTTCGACTCCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGc                                                                         <  2:509798/147‑1 (MQ=255)                                                                           cTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAt                                                                                                           <  1:1293061/115‑1 (MQ=255)                                                                           cTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCAt                                                                                                           >  2:1293061/1‑115 (MQ=255)                                                                            tACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGAtctttc                                                                                                                                                 <  1:1572317/76‑1 (MQ=255)                                                                            tACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGAtctttc                                                                                                                                                 >  2:1572317/1‑76 (MQ=255)                                                                                      cTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCCTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGgcgc                                                             >  1:1245952/1‑150 (MQ=255)                                                                                      cTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGAAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCAGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGgcgc                                                             >  1:1043623/1‑150 (MQ=255)                                                                                        gtgCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCt                                                                        <  1:1369080/137‑1 (MQ=255)                                                                                        gtgCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCt                                                                        >  2:1369080/1‑137 (MQ=255)                                                                                           cTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGc                                                        >  1:833747/1‑150 (MQ=255)                                                                                               gACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGtt                                                                                                               >  2:157722/1‑91 (MQ=255)                                                                                               gACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGtt                                                                                                               <  1:157722/91‑1 (MQ=255)                                                                                                    cAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccacccc                                               >  2:834750/1‑150 (MQ=255)                                                                                                      aGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCCTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGcccccaccccca                                             <  2:1245952/150‑1 (MQ=255)                                                                                                       ggAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAg                                            >  2:1691855/1‑150 (MQ=255)                                                                                                        gacaTTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGc                                           >  1:2400849/4‑150 (MQ=255)                                                                                                          aaTTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCgg                                         >  2:418643/1‑150 (MQ=255)                                                                                                                cGAGAAAAAGGCAAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGCGGCGCATCATGCGCCCCGTTTCCTTTATATGCCCAGCCTGACCCTGACCGACTCCCGCTACATCGCGCGCTTCGCCCCCACCCCAAGCAGCTTcc                                    >  2:280348/1‑150 (MQ=255)                                                                                                                cGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                   >  2:1496531/1‑150 (MQ=255)                                                                                                                cGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCt                                   <  1:1199635/150‑1 (MQ=255)                                                                                                                cGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCCTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGc                                                                                              <  2:2151935/91‑1 (MQ=255)                                                                                                                cGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCCTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGc                                                                                              >  1:2151935/1‑91 (MQ=255)                                                                                                                            aaaaaaCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCt                        >  2:297028/1‑150 (MQ=255)                                                                                                                            aaaaaaCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACCGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCt                        <  2:1650317/150‑1 (MQ=255)                                                                                                                               aaaCGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGc                     <  1:1792509/150‑1 (MQ=255)                                                                                                                                aaCGGTCGGCAAAGGCGGATCTTTCTTCCACGTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCt                    >  1:195894/1‑150 (MQ=255)                                                                                                                                    gTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTc                <  1:916696/150‑1 (MQ=255)                                                                                                                                         cAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTcgccgc           >  1:264455/1‑150 (MQ=255)                                                                                                                                         cAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTcgccgc           >  1:2046528/1‑150 (MQ=255)                                                                                                                                              gCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCATGCACTTCGGCTCGCTGGTCGCCGCCCTCg      >  2:1671964/1‑150 (MQ=255)                                                                                                                                                tGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCTCAGCTTAAAAAATAAGCCCAGCCGGGCCATGACCTACTACAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcc    >  1:1855434/1‑150 (MQ=255)                                                                                                                                                  aTCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGcctc  >  1:274991/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GATTCAGGACGAAGAGTACGGCTGGTGCGAGTCGTGCGGCATCGAGATCGGCCTGCGCCGCCTGGAAGCCCGCCCTACCGCCGACCTGTGCTTCGACTGCAAGGAAATTGCCGAGAAAAAGG‑AAAAAACGGTCGGCAAAGGCTGATCTTTCTTCCACCTGAACGGGGCGCATCATGCGCCCCGTTTCATTTATATGCCCAGCCTGACCATGACCGACTCCAGCTACATCGGGCGCTTCGCCCCCACCCCCAGCGGCTTCCTGCACTTCGGCTCGCTGGTCGCCGCCCTCGCCTC  >  NC_002947/5332121‑5332414

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