Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,964,429	A→G	28.8%	F310S (TTC→TCC)	fleQ ←	transcriptional regulator FleQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,964,429	0	A	G	28.8%	87.3 / 45.4	66	F310S (TTC→TCC)	fleQ	transcriptional regulator FleQ
Reads supporting (aligned to +/- strand):  ref base A (24/23);  new base G (8/11);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

AGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGC  >  NC_002947/4964286‑4964572                                                                                                                                                |                                                                                                                                                aGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACg                                                                                                                                            >  5:287727/1‑149 (MQ=255)    cGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGtt                                                                                                                                          >  1:67654/1‑148 (MQ=255)          gacACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGGACACGTTCAGGCgg                                                                                                                                   >  6:301046/1‑149 (MQ=255)               aGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGgtagta                                                                                                                              >  5:18962/1‑149 (MQ=255)               aGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGgtag                                                                                                                                <  5:85888/147‑1 (MQ=255)               aGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGgtag                                                                                                                                >  6:85888/1‑147 (MQ=255)                    cATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAgg                                                                                                                          >  4:77592/1‑148 (MQ=255)                      tGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCt                                                                                                                       <  8:284068/149‑1 (MQ=255)                       gATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCtt                                                                                                                      >  2:325444/1‑149 (MQ=255)                        aTCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGCCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTc                                                                                                                     <  7:190982/149‑1 (MQ=255)                           gAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTcgcg                                                                                                                  <  5:87927/149‑1 (MQ=255)                            aGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTcgcg                                                                                                                  >  6:45772/1‑148 (MQ=255)                               cGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTGCAGGTCTTCGCGGAAg                                                                                                              <  2:294894/149‑1 (MQ=255)                                 gAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCgg                                                                                                                                   >  4:149579/1‑126 (MQ=255)                                 gAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCgg                                                                                                                                   <  3:149579/126‑1 (MQ=255)                                  aaTTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGTATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTg                                                                                                           <  5:23507/149‑1 (MQ=255)                                  aaTTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTg                                                                                                           <  3:20696/149‑1 (MQ=255)                                   aTTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCAAGCAGCGGCGCCATTTCGATAGGGAACACGTTAAGGCGGTAGTACAGGTCTTCGCGGAAGGTGc                                                                                                          >  7:15458/1‑149 (MQ=255)                                    ttGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGcc                                                                                                         <  7:38381/149‑1 (MQ=255)                                    ttGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGCCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGcc                                                                                                         <  4:228780/149‑1 (MQ=255)                                       aaGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTc                                                                                                      >  4:51909/1‑149 (MQ=255)                                         gCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCt                                                                                                    >  2:195022/1‑149 (MQ=255)                                          cGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTc                                                                                                   >  2:283824/1‑149 (MQ=255)                                          cGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACAAGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTc                                                                                                   >  2:283733/1‑149 (MQ=255)                                                      cgcTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTc                                                                                       <  1:195022/149‑1 (MQ=255)                                                           ctcGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAgg                                                                                   <  5:301046/148‑1 (MQ=255)                                                              gTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTc                                                                                <  7:89710/148‑1 (MQ=255)                                                               tGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGCCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCtt                                                                              <  6:78864/149‑1 (MQ=255)                                                                 cTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTtgt                                                                            <  6:34941/149‑1 (MQ=255)                                                                   ccATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGc                                                                          >  4:87363/1‑149 (MQ=255)                                                                   ccATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGc                                                                          <  1:325444/149‑1 (MQ=255)                                                                    cATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGc                                                                          <  6:287727/148‑1 (MQ=255)                                                                    cATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCAGTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTtgtgag                                                                         >  6:138767/1‑147 (MQ=255)                                                                         gcgAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGt                                                                        <  4:101006/145‑1 (MQ=255)                                                                         gcgAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGt                                                                        >  3:101006/1‑145 (MQ=255)                                                                            aaaTCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGc                                                                  >  2:39055/1‑148 (MQ=255)                                                                                  gCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGatgatg                                                           <  3:77592/149‑1 (MQ=255)                                                                                    tCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATgcg                                                         >  5:250005/1‑149 (MQ=255)                                                                                      gTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCGTTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCa                                                       >  6:195162/1‑149 (MQ=255)                                                                                          atcaGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGt                                                    <  7:157936/148‑1 (MQ=255)                                                                                          atcaGCAACGGGATGCCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCgatgat                                                            >  6:320404/1‑140 (MQ=255)                                                                                          atcaGCAACGGGATGCCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCgatgat                                                            <  5:320404/140‑1 (MQ=255)                                                                                          atcaGCAACGGGATGCCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTc                                                   >  3:215186/1‑149 (MQ=255)                                                                                           tcaGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTc                                                   >  5:305804/1‑148 (MQ=255)                                                                                             aGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGAt                                                >  5:131233/1‑149 (MQ=255)                                                                                              gCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATg                                               >  3:181497/1‑149 (MQ=255)                                                                                                aaCGGGATGCCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGACGATGCt                                             >  2:305346/1‑149 (MQ=255)                                                                                                    ggATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTgg                                         <  3:51909/149‑1 (MQ=255)                                                                                                    ggATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTgg                                         >  4:305190/1‑149 (MQ=255)                                                                                                       tGTCTTCCTCCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTc                                      <  6:18962/149‑1 (MQ=255)                                                                                                         tCTTCCACCCGCTCACGCAGCGGCGCCGTTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCtt                                    <  5:195162/149‑1 (MQ=255)                                                                                                          cTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCttg                                   >  7:61283/1‑149 (MQ=255)                                                                                                           ttCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCttgt                                  >  1:339480/1‑149 (MQ=255)                                                                                                             ccACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCttgtt                                 <  3:311942/148‑1 (MQ=255)                                                                                                             caaCCCGCTCACGCAGCGGCGCCTTTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTGGCGGAAGGTGCCGTCCTCGATCATGTTTTCGAGGTTCTTGTGCGTCGCAGCGTTGATGCGGACGTCGATGCTTTGGGTCttgt                                  <  7:116862/145‑1 (MQ=255)                                                                                                                         gCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACgcgc                    >  2:119590/1‑149 (MQ=255)                                                                                                                               gcgcCATTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTTGGTCttgt                                  >  6:312828/1‑129 (MQ=255)                                                                                                                               gcgcCATTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTTGGTCttgt                                  <  5:312828/129‑1 (MQ=255)                                                                                                                               gcgcCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAg              <  8:61283/149‑1 (MQ=255)                                                                                                                                cgcCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAgg             <  4:181497/149‑1 (MQ=255)                                                                                                                                 gcCATTTCGATAGGGGACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGt            <  6:131233/149‑1 (MQ=255)                                                                                                                                 gcCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCAGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTgg                                         <  7:200566/120‑1 (MQ=255)                                                                                                                                 gcCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCAGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTgg                                         >  8:200566/1‑120 (MQ=255)                                                                                                                                        cGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTc      <  1:119590/148‑1 (MQ=255)                                                                                                                                         gATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCt    <  2:339480/149‑1 (MQ=255)                                                                                                                                           tAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTgc  <  3:305190/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                AGCCGTGACGACACAGCGACATGATCGAGGCGGAATTGAAGCGAATGGAACCGCGCTTCTCGTGCTCCATGCGCGAAATCAGCTCGTTCATCAGCAACGGGATGTCTTCCACCCGCTCACGCAGCGGCGCCATTTCGATAGGGAACACGTTCAGGCGGTAGTACAGGTCTTCGCGGAAGGTGCCGTCCTCGATCATGGTTTCGAGGTTCTTGTGCGTCGCAGCGATGATGCGCACGTCGATGCTTTGGGTCTTGTTGCTACCCACGCGCTCGAAGGTACGCTCCTGC  >  NC_002947/4964286‑4964572

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