Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,762,942	A→T	8.8%	N62I (AAC→ATC)	PP_3326 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,762,942	0	A	T	8.8%	182.4 / 8.1	68	N62I (AAC→ATC)	PP_3326	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/29);  new base T (2/4);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

GTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCAT  >  NC_002947/3762800‑3763087                                                                                                                                                |                                                                                                                                                   gTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGAT‑CATg                                                                                                                                               <  3:2173/148‑1 (MQ=255)  tCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATc                                                                                                                                            <  1:144877/149‑1 (MQ=255)      cGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGTCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCgg                                                                                                                                        >  1:114452/1‑149 (MQ=255)      cGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCgg                                                                                                                                        <  3:366755/149‑1 (MQ=255)       gTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGt                                                                                                                                       >  8:216081/1‑149 (MQ=255)       gTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGt                                                                                                                                       >  2:104083/1‑149 (MQ=255)        tCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGt                                                                                                                                       >  7:8823/1‑148 (MQ=255)          aaCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCAc                                                                                                                                    >  8:317348/1‑149 (MQ=255)              cAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTc                                                                                                                                <  8:194063/149‑1 (MQ=255)                 gcgcCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAT‑CATGTTCCCGGTCACgtcgaac                                                                                                                             <  8:194223/149‑7 (MQ=255)                  cgcCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAgc                                                                                                                            <  4:176200/149‑1 (MQ=255)                    cccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAgcgc                                                                                                                          <  5:38572/149‑1 (MQ=255)                     ccACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGcc                                                                                                                         >  6:69527/1‑149 (MQ=255)                          gCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCg                                                                                                                    <  3:114736/149‑1 (MQ=255)                            cTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGa                                                                                                                  <  8:8823/149‑1 (MQ=255)                               tgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGGTCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATc                                                                                                                <  5:34426/148‑1 (MQ=255)                                ggACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCAt                                                                                                              <  7:216081/149‑1 (MQ=255)                                 gACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATg                                                                                                             >  1:185049/1‑149 (MQ=255)                                             aCGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAAc                                                                                                 >  7:280277/1‑149 (MQ=255)                                             aCGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAAc                                                                                                 >  2:299392/1‑149 (MQ=255)                                                 ggCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGAT‑CATGTTCCCGGTCa                                                                                                                                     <  6:312517/110‑1 (MQ=38)                                                 ggCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGAT‑CATGTTCCCGGTCa                                                                                                                                     >  5:312517/1‑110 (MQ=38)                                                  gCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACt                                                                                             >  5:20558/1‑148 (MQ=255)                                                   cAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGcc                                                                                                           <  3:197743/133‑1 (MQ=255)                                                   cAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGcc                                                                                                           >  4:197743/1‑133 (MQ=255)                                                   cAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTa                                                                                           >  4:292052/1‑149 (MQ=255)                                                     atatCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACCTTGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAc                                                                                          <  1:104083/149‑1 (MQ=255)                                                      tatCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTAccc                                                                                        >  7:268578/1‑149 (MQ=255)                                                         cGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGc                                                                                     >  8:217266/1‑149 (MQ=255)                                                           aGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGATCGCCTCGCGGATCATGCCGTGGTGGCACCACTTACCCCGCCt                                                                                   >  5:51578/1‑149 (MQ=255)                                                                   tACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCAc                                                                            <  7:317348/148‑1 (MQ=255)                                                                     cGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTg                                                                         <  8:280277/149‑1 (MQ=255)                                                                        ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                                      <  2:114452/149‑1 (MQ=255)                                                                        ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                                      >  6:331465/1‑149 (MQ=255)                                                                         cACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGt                                                                     >  3:264809/1‑149 (MQ=255)                                                                                       tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACg                                                       >  7:18358/1‑149 (MQ=255)                                                                                         aCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACgg                                                      <  6:19554/148‑1 (MQ=255)                                                                                           gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCAc                                                                                              <  2:154198/106‑1 (MQ=255)                                                                                           gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCAc                                                                                              >  1:154198/1‑106 (MQ=255)                                                                                           gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgctgc                                                                <  7:320744/136‑1 (MQ=255)                                                                                           gCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtgctgc                                                                >  8:320744/1‑136 (MQ=255)                                                                                             cTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCtt                                                 <  7:324354/149‑1 (MQ=255)                                                                                               aCAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                                                <  3:350040/148‑1 (MQ=255)                                                                                               aCAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAAAGGGCCTTg                                                >  4:350040/1‑148 (MQ=255)                                                                                                 aaGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGc                                                                                 >  3:342109/1‑113 (MQ=255)                                                                                                 aaGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGc                                                                                 <  4:342109/113‑1 (MQ=255)                                                                                                 aaGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCCGGAA‑CATGATACCGGTAACGGGCCGGCGCCGGGCGGGTCATGCCGTGCTGGCACCACTTAACCCCCCTGCGGCACATGCCGTGCTGCGGGGCAACCGGCCTTGGt                                              >  6:50160/1‑148 (MQ=255)                                                                                                   ggCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTca                                            >  2:35134/1‑148 (MQ=255)                                                                                                      ggAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgc                                        >  3:269936/1‑149 (MQ=255)                                                                                                       gAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgcg                                       >  2:50827/1‑149 (MQ=255)                                                                                                       gAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAgcg                                       >  1:37519/1‑149 (MQ=255)                                                                                                            gATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCCCGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCg                                   >  7:38630/1‑148 (MQ=255)                                                                                                                aaGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGTCAACGGGCCTTGGTCACAGCGCCCGGACaa                              <  1:35134/149‑1 (MQ=255)                                                                                                                  gcgcgcCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAg                             >  5:32372/1‑148 (MQ=255)                                                                                                                      gcCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATg                        <  6:20558/149‑1 (MQ=255)                                                                                                                      cgccaggCGCTCGACCGTGACCGGGAT‑CATGTTCCCGGTCAC‑GTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATg                        >  4:45687/6‑149 (MQ=37)                                                                                                                       ccTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                       <  8:268578/149‑1 (MQ=255)                                                                                                                        cTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATggg                      <  6:32372/149‑1 (MQ=255)                                                                                                                          gggCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGAGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGt                     >  3:126177/1‑148 (MQ=255)                                                                                                                              gTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                          >  7:175760/1‑123 (MQ=255)                                                                                                                              gTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTcaca                                          <  8:175760/123‑1 (MQ=255)                                                                                                                                tCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCGGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAg              >  5:196107/1‑149 (MQ=255)                                                                                                                                 cGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAg              <  1:299392/148‑1 (MQ=255)                                                                                                                                   a‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAt           <  8:221480/149‑1 (MQ=255)                                                                                                                                   a‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAt           >  8:3871/1‑149 (MQ=255)                                                                                                                                   a‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGAt           <  4:269936/149‑1 (MQ=255)                                                                                                                                       tGACCGGGAT‑CATGTTCCCGGTCAC‑GTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCa       <  3:45687/149‑1 (MQ=255)                                                                                                                                             ggAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCAt  <  2:37519/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                   GTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAA‑CATGATCCCGGTCACGGTCGAGCGCCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCAT  >  NC_002947/3762800‑3763087

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