Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,762,942	A→T	7.6%	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	7.6%	187.9 / 5.2	66	N62I (AAC→ATC)	PP_3326	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/32);  new base T (2/3);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.08e-01

ATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACC  >  NC_002947/3762797‑3763084                                                                                                                                                   |                                                                                                                                                  aTGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACa                                                                                                                                                 <  8:271473/148‑1 (MQ=255)     tCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATc                                                                                                                                            >  4:72706/1‑149 (MQ=255)       aGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATccc                                                                                                                                          >  3:178823/1‑149 (MQ=255)          gTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGATCATGTTCCCggg                                                                                                                                       <  5:77047/149‑2 (MQ=255)             aaCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTTGCAATATCGAGGCCGTGTCCGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGATCATGTTCCCGGTCa                                                                                                                                     <  5:306987/148‑1 (MQ=255)             aaCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCAc                                                                                                                                    <  8:157337/149‑1 (MQ=255)                  aGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTc                                                                                                                               >  8:20524/1‑148 (MQ=255)                           cGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCC‑AGGCGCTCGACCGTGACCGGGATCATGTTCCCGGTCAC‑‑GTCGAACG‑cccag                                                                                                                    <  1:115456/148‑4 (MQ=25)                                  tgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCa                                                                                                             >  6:324508/1‑149 (MQ=255)                                     aCCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGc                                                                                                          >  8:212868/1‑149 (MQ=255)                                          ggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGtgc                                                                                                     <  2:71037/149‑1 (MQ=255)                                             gCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGtgctgc                                                                                                  <  5:219826/149‑1 (MQ=255)                                              cGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGtgct                                                                                                    <  6:73703/146‑1 (MQ=255)                                              cGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGtgct                                                                                                    >  5:73703/1‑146 (MQ=255)                                                   tGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCAc                                                                                            >  8:19838/1‑149 (MQ=255)                                                     gCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATACCGGTCACG‑GTCGAGCG‑CCTGGGGGATCAATCCGTGCTGCAACCAca                                                                                           >  6:56003/1‑147 (MQ=255)                                                      cAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTa                                                                                         <  7:212868/149‑1 (MQ=255)                                                       aaTATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTAc                                                                                        >  4:11275/1‑149 (MQ=255)                                                          atCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACtt                                                                                          >  2:300924/1‑144 (MQ=255)                                                          atCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACtt                                                                                          <  1:300924/144‑1 (MQ=255)                                                            cGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGc                                                                                   >  2:241907/1‑149 (MQ=255)                                                                 ccGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCg                                                                              <  7:19838/149‑1 (MQ=255)                                                                  cttgtACGGCCACGGCGCCGAGTATAACGCCTACAAGGCGGAAAAGATCATGCGCGCCTGGGCGTTCGA‑CGTGACCGGTAACATGATCCCGGTCAGGTGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCg                                                                              <  1:241907/147‑1 (MQ=255)                                                                   gtgtACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCgg                                                                             >  5:168207/1‑148 (MQ=255)                                                                    tgtACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGAG‑ACTGGCGGATCATGCCGTGCAGCAACCACTTACCCCGCCTGCGGCa                                                                           >  1:220493/1‑149 (MQ=255)                                                                        cGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTg                                                                       >  5:218246/1‑149 (MQ=255)                                                                         ggCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGc                                                                      >  5:13639/1‑149 (MQ=255)                                                                         ggCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGc                                                                      >  4:65024/1‑149 (MQ=255)                                                                           ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                                    >  6:161886/1‑149 (MQ=255)                                                                            cACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCCGGAACATGATACCGGTCACG‑GTCGACCG‑CCTGGCGGATCATGCCGTGATGCAACCACATAACCCGCCTGCGGCACTTGCCG‑t                                                                  >  4:150711/1‑149 (MQ=255)                                                                                    cGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTgg                                                          <  4:68922/149‑1 (MQ=255)                                                                                    cGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTgg                                                          <  3:72706/149‑1 (MQ=255)                                                                                       ggATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGAc                                                        >  2:210503/1‑148 (MQ=255)                                                                                          tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACg                                                    <  3:65024/149‑1 (MQ=255)                                                                                            aCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACcgg                                                  >  6:231643/1‑149 (MQ=255)                                                                                             cGCCTAGAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGc                                                 <  8:107343/149‑1 (MQ=255)                                                                                             cGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGc                                                 <  6:168207/149‑1 (MQ=255)                                                                                               catACAAGGCGTAAAAGATCAAGCCCGCATTGGCGTTCGA‑CGTGACCGGGAACATGATCCCGTTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCt                                               <  8:263258/147‑1 (MQ=255)                                                                                                   cAAGGCGTAATAGATCTCGCGCGCCTGGGCGTTCGA‑CGTGACCGTTAACTTGATCCCTGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGAAACCACTTACCCCGCCTGCGTCACTTTCCG‑TGCTGCTGGACACCGGTCCTTgc                                            <  2:220493/148‑2 (MQ=255)                                                                                                   cAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTg                                             <  8:198932/147‑1 (MQ=255)                                                                                                   cAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTg                                             <  8:194154/147‑1 (MQ=255)                                                                                                   cAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTg                                             <  5:8287/147‑1 (MQ=255)                                                                                                       gCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTcaca                                       <  5:324508/149‑1 (MQ=255)                                                                                                       gCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTcaca                                       <  8:219849/149‑1 (MQ=255)                                                                                                         ggAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAgc                                     >  5:204094/1‑149 (MQ=255)                                                                                                             aaGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGccc                                 <  7:20524/149‑1 (MQ=255)                                                                                                               gATCAAGCGGGCCTGGGGGTTCGA‑CGTGTCCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGCTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCg                                <  1:210503/149‑1 (MQ=255)                                                                                                                aTCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGa                              >  1:52303/1‑149 (MQ=255)                                                                                                                   aaGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGAc                             >  5:78257/1‑147 (MQ=255)                                                                                                                        cgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAAt                      <  2:269352/149‑1 (MQ=255)                                                                                                                        cgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAAt                      <  6:204094/149‑1 (MQ=255)                                                                                                                        cgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAAACACTTACCCCGCATGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAAt                      <  5:231643/149‑1 (MQ=255)                                                                                                                          ccTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                    <  3:11275/149‑1 (MQ=255)                                                                                                                          ccTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                    <  5:161886/149‑1 (MQ=255)                                                                                                                          ccTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                    <  5:56003/149‑1 (MQ=255)                                                                                                                                  ttCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAAt              >  1:142704/1‑147 (MQ=255)                                                                                                                                     gA‑CGTGACCGGGAACATGATCCGGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACACCGGGCCTTGGACCCAGCGCCCGGACAAGCAAGGGGTCAATCAgg          >  7:203001/1‑148 (MQ=255)                                                                                                                                     gA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGa         <  4:178823/149‑1 (MQ=255)                                                                                                                                      a‑CGTGACCGGGAACATGATCACGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCAGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTAGGACACAGCGCCCGGACAAGCAATGGGACAATCAGGAt        >  8:195907/1‑149 (MQ=255)                                                                                                                                           gACCGGGATCATGTTCCCGGTCAC‑‑GTCGAACGCCCAGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCa    >  2:114362/1‑149 (MQ=38)                                                                                                                                            aCCGGGATCATGTTCCCGGTCAC‑‑GTCGAACGCCCAGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTAGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCAc   >  7:132939/1‑149 (MQ=38)                                                                                                                                             ccGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGAc                             <  3:103520/122‑1 (MQ=255)                                                                                                                                             ccGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGAc                             >  4:103520/1‑122 (MQ=255)                                                                                                                                              cGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAg                          >  3:129220/1‑124 (MQ=255)                                                                                                                                              cGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAg                          <  4:129220/124‑1 (MQ=255)                                                                                                                                              cGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCAcc  >  2:68889/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                  ATGGTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACC  >  NC_002947/3762797‑3763084

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