Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,762,947	A→T	7.3%	I64F (ATC→TTC)	PP_3326 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,762,947	0	A	T	7.3%	194.3 / 5.1	69	I64F (ATC→TTC)	PP_3326	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base T (2/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.43e-01

TCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAG  >  NC_002947/3762801‑3763093                                                                                                                                                    |                                                                                                                                                      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                                                                                                             >  5:73703/1‑146 (MQ=255)                                          cGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGtgct                                                                                                             <  6:73703/146‑1 (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                                                                   <  3:72706/149‑1 (MQ=255)                                                                                cGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTgg                                                                   <  4:68922/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:194154/147‑1 (MQ=255)                                                                                               cAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTg                                                      <  8:198932/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                               <  6:204094/149‑1 (MQ=255)                                                                                                                    cgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAAt                               <  2:269352/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:56003/149‑1 (MQ=255)                                                                                                                      ccTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                             <  5:161886/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)                                                                                                                                                  tGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCTGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAg  <  7:195907/149‑1 (MQ=255)                                                                                                                                                   gATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCaa   >  2:4475/1‑147 (MQ=255)                                                                                                                                                   gATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCaa   <  1:4475/147‑1 (MQ=255)                                                                                                                                                   gATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAg  >  4:238045/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                      TCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCCCGGTCACG‑GTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCG‑TGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAG  >  NC_002947/3762801‑3763093

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