Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,762,947	A→T	11.1%	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	11.1%	158.2 / 10.5	63	I64F (ATC→TTC)	PP_3326	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base T (4/3);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.11e-01

GTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAG  >  NC_002947/3762800‑3763093                                                                                                                                                     |                                                                                                                                                      gTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGAt                                                                                                                                                      <  1:166914/149‑1 (MQ=255) gTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGAt                                                                                                                                                      >  6:36057/1‑149 (MQ=255)    aGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATc‑‑cc                                                                                                                                                 <  4:109928/149‑1 (MQ=255)      cGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCgg                                                                                                                                               >  3:61459/1‑149 (MQ=255)      cGTAAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGAAc‑‑cccg                                                                                                                                               >  1:344304/1‑147 (MQ=255)       gTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACCGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTC‑‑CCgg                                                                                                                                               <  4:227686/149‑1 (MQ=255)          aaCGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGCGACCGGGAACATGATC‑‑CCGGTCAc                                                                                                                                           >  7:305682/1‑149 (MQ=255)               aGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTCTACGGCCACGGCGCCGAGGATAACGACTACAAGGCGGAAAAGATCAAGCGCTCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTc                                                                                                                                       >  7:256404/1‑148 (MQ=255)                            cTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGa                                                                                                                        >  1:220716/1‑149 (MQ=255)                             tgtgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTC‑‑CCgg                                                                                                                                               >  1:326333/1‑127 (MQ=39)                             tgtgGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTC‑‑CCgg                                                                                                                                               <  2:326333/127‑1 (MQ=39)                             tgtgGACAGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGCTCGACCGTGACCGGGATCATGTTC‑‑CCgct                                                                                                                                              >  5:199724/1‑126 (MQ=38)                                 gACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATg                                                                                                                   <  2:220716/149‑1 (MQ=255)                                     gaggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGACGGATCATGCCGt                                                                                                               >  1:191266/1‑149 (MQ=255)                                      aggagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGtg                                                                                                              <  5:157355/149‑1 (MQ=255)                                        gagCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGCACATGATC‑‑CCGGTCACGGGCGAGCG‑CCTGCCGGATCATGCCGagct                                                                                                            >  2:177124/1‑149 (MQ=255)                                         agCGACGTGGAAATATCGAGGCCGGGTACGGCCTCGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGtgctg                                                                                                           <  2:344304/149‑1 (MQ=255)                                           cGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGtgctgc                                                                                                          >  1:44794/1‑148 (MQ=255)                                               gTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCa                                                                                                     <  5:36057/149‑1 (MQ=255)                                                tGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCAc                                                                                                    >  7:213594/1‑149 (MQ=255)                                                     atatCGAGGCCGTGTACGGCCTCGGCGCCGAGGTTAACGGCTAGAAGGCGGAAAAGATCAAGCGCGCCTGGGAGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGACCATGCCGTGCTGCAACGACTTAcc                                                                                               <  5:234523/149‑1 (MQ=255)                                                         cGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGc                                                                                           >  3:144085/1‑149 (MQ=255)                                                          gAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGc                                                                                           <  4:277695/148‑1 (MQ=255)                                                                      ggCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGc                                                                              <  8:124526/149‑1 (MQ=255)                                                                      ggCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGc                                                                              >  3:336397/1‑149 (MQ=255)                                                                        ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                                            >  2:23909/1‑149 (MQ=255)                                                                        ccACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                                            >  5:387725/1‑149 (MQ=255)                                                                         cACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGt                                                                           <  4:88985/149‑1 (MQ=255)                                                                          aCGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGtg                                                                          <  2:194986/149‑1 (MQ=255)                                                                                ccGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGctgctg                                                                    <  4:144085/149‑1 (MQ=255)                                                                                  gAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTg                                                                               <  5:219600/136‑1 (MQ=255)                                                                                  gAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTg                                                                               >  6:219600/1‑136 (MQ=255)                                                                                       tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                                            >  8:167677/1‑134 (MQ=255)                                                                                       tAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCg                                                                            <  7:167677/134‑1 (MQ=255)                                                                                        aaCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATCAGGCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGCCAAc                                                              >  1:113378/1‑149 (MQ=255)                                                                                         aCGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACgg                                                            <  2:44794/148‑1 (MQ=255)                                                                                              tACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                                                      >  1:295581/1‑149 (MQ=255)                                                                                                cAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTg                                                      <  7:308463/147‑1 (MQ=255)                                                                                                          aaGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGccc                                          <  4:61459/149‑1 (MQ=255)                                                                                                          aaGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGccc                                          >  2:96966/1‑149 (MQ=255)                                                                                                          aaGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGACGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGccc                                          <  2:191266/149‑1 (MQ=255)                                                                                                             aTCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGa                                       >  1:229539/1‑149 (MQ=255)                                                                                                                aaGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACaa                                    <  6:218476/149‑1 (MQ=255)                                                                                                                aaGCGCGCCTGGGCGTTCGA‑CGCGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACaa                                    <  8:305682/149‑1 (MQ=255)                                                                                                                    gcgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACg                                                             >  4:218623/1‑120 (MQ=255)                                                                                                                    gcgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACg                                                             <  3:218623/120‑1 (MQ=255)                                                                                                                    gcgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACg                                                             <  3:337707/120‑1 (MQ=255)                                                                                                                    gcgcCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACg                                                             >  4:337707/1‑120 (MQ=255)                                                                                                                       ccTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCt                                                                                         >  3:24949/1‑89 (MQ=255)                                                                                                                       ccTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCt                                                                                         <  4:24949/89‑1 (MQ=255)                                                                                                                       ccTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATgg                             <  8:213594/149‑1 (MQ=255)                                                                                                                        cTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATggg                            >  3:267441/1‑149 (MQ=255)                                                                                                                        cTGGGCGGTAGA‑CGGTACCGGGAACATGATC‑‑CCGGGCACGGGCGAGCG‑CCTGGCGGAGCATGCCGTGCTGCAACGACTTACCCCGCCTGCGGCAATTGCCGGGCGGCTGGACAACGGTCCGTGGTAACAGCGCCCGGACAATCAATggg                            <  2:194606/149‑1 (MQ=255)                                                                                                                         tGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGt                           >  4:323501/1‑149 (MQ=255)                                                                                                                                   accGTGACCGGGATCATGTTC‑‑CCGGTCAC‑GTCGAACGCCCAGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAgg                   >  8:337444/3‑148 (MQ=37)                                                                                                                                       tGACCGGGATCATGTTC‑‑CCGGTCAC‑GTCGAACG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACa                                     >  2:143440/1‑125 (MQ=39)                                                                                                                                       tGACCGGGATCATGTTC‑‑CCGGTCAC‑GTCGAACG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACa                                     <  1:143440/125‑1 (MQ=39)                                                                                                                                       tGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCg               >  7:236779/1‑148 (MQ=255)                                                                                                                                           cGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCAccc          <  4:336397/149‑1 (MQ=255)                                                                                                                                             ggAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCAt        >  1:54390/1‑149 (MQ=255)                                                                                                                                                   tGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAg  <  2:54390/149‑1 (MQ=255)                                                                                                                                                   tGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAg  <  2:113378/149‑1 (MQ=255)                                                                                                                                                    gATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAg  >  6:54240/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GTCAGCGTCAACGCAGGCGCCCACGGCCTGTGGACCGAGGAGCGACGTGGCAATATCGAGGCCGTGTACGGCCACGGCGCCGAGGATAACGCCTACAAGGCGGAAAAGATCAAGCGCGCCTGGGCGTTCGA‑CGTGACCGGGAACATGATC‑‑CCGGTCACGGTCGAGCG‑CCTGGCGGATCATGCCGTGCTGCAACCACTTACCCCGCCTGCGGCACTTGCCGTGCTGCTGGACAACGGGCCTTGGTCACAGCGCCCGGACAAGCAATGGGTCAATCAGGATCGCACCCATGTCAAG  >  NC_002947/3762800‑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: ‑