Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,510,737	C→A	6.7%	I166I (ATC→ATA)	PP_5575 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,510,737	0	C	A	6.7%	192.9 / 4.0	60	I166I (ATC→ATA)	PP_5575	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (28/28);  new base A (2/2);  total (30/30)
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.62e-01

TTCGCCAACGTGACTGCCTCCTGGCTGGGAACCCGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTATCAGGAATTTTCATCGCATGGATTACGAATTACAGGAAGG  >  NC_002947/3510592‑3510881                                                                                                                                                     |                                                                                                                                                  ttCGCCAACGTGACTGCCTCCTGGCTGGGAACCCGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCCGGAGCGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATAtc                                                                                                                                                   >  1:108556/1‑149 (MQ=255)     ccAACGTGACTGCCTCCTGGCTGGGAACCCGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCg                                                                                                                                            >  8:1878/1‑149 (MQ=255)      cAACGTGACTGCCTCCTGGCTGGGAACCCGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGt                                                                                                                                           >  1:38697/1‑149 (MQ=255)       aaCGTGACTGCCTCCTGGCTGGGAACCCGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTg                                                                                                                                          <  6:200551/149‑1 (MQ=255)          gTGACTGCCTCCTGGCTGGGAACCCGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATc                                                                                                                                       <  5:161397/149‑1 (MQ=255)            gACTGCCTCCTGGCTGGGAACCCGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGa                                                                                                                                     >  4:18665/1‑149 (MQ=255)            gACTGCCTCCTGGCTGGGAACCCGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGa                                                                                                                                     >  1:79131/1‑149 (MQ=255)               tGCCTCCTGGCTGGGAACCCGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGAccc                                                                                                                                  >  3:165060/1‑149 (MQ=255)                 cctcctGGCTGGGAACCCGGCCAAGGTCGTGGACTTTCACTTCAATGACGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGAccc                                                                                                                                  >  6:155596/1‑147 (MQ=255)                  ctcctGGCTGGGAACCCGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGTGATATCCAGccattgca                                                                                                                                      <  1:172393/142‑8 (MQ=255)                  ctcctGGCTGGGAACCCGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGTGATATCCAGccattgca                                                                                                                                      >  2:172393/1‑135 (MQ=255)                      tggctggGAACCCGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAg                                                                                                                           >  4:28404/1‑149 (MQ=255)                           gggAACCCGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCgg                                                                                                                      >  8:104864/1‑149 (MQ=255)                                cccGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGa                                                                                                                 <  5:155596/149‑1 (MQ=255)                                cccGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGa                                                                                                                 <  4:120207/149‑1 (MQ=255)                                      cAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAg                                                                                                                               <  6:283576/129‑1 (MQ=255)                                      cAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAg                                                                                                                               >  5:283576/1‑129 (MQ=255)                                          gTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGtggct                                                                                                       <  7:104864/149‑1 (MQ=255)                                                gACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGa                                                                                                 <  3:18665/149‑1 (MQ=255)                                                 aCTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAg                                                                                                <  2:79131/149‑1 (MQ=255)                                                  cTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGt                                                                                               <  2:75389/149‑1 (MQ=255)                                                  cTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGATGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGt                                                                                               >  4:213400/1‑149 (MQ=255)                                                   tttCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGtt                                                                                              >  1:219875/1‑149 (MQ=255)                                                       aCTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGa                                                                                           <  2:108556/148‑1 (MQ=255)                                                         ttCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCg                                                                                        >  3:207698/1‑149 (MQ=255)                                                                 cGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGtggctgg                                                                                                     >  4:76322/1‑128 (MQ=255)                                                                 cGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGtggctgg                                                                                                     <  3:76322/128‑1 (MQ=255)                                                                 cGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTc                                                                                 >  4:87492/1‑148 (MQ=255)                                                                   gCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGAt                                                                              <  2:219875/149‑1 (MQ=255)                                                                     c‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATgg                                                                            >  5:106324/1‑149 (MQ=255)                                                                            aGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTc                                                                        <  3:28404/149‑1 (MQ=255)                                                                              ccTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTcg‑c                                                                     >  5:180838/1‑149 (MQ=255)                                                                                tCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTcg‑cg                                                                    >  5:195502/1‑148 (MQ=255)                                                                                  ggCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAg                                                                 <  6:180838/149‑1 (MQ=255)                                                                                   gCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGtgg                                                                                                         >  8:98817/1‑109 (MQ=255)                                                                                   gCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGtgg                                                                                                         <  7:98817/109‑1 (MQ=255)                                                                                         tCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGAc                                                          >  6:45419/1‑149 (MQ=255)                                                                                           gCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGAccc                                                        <  6:106324/149‑1 (MQ=255)                                                                                                gTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGAt                                                    <  4:165060/148‑1 (MQ=255)                                                                                                     cgcACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTc                                              <  7:1878/149‑1 (MQ=255)                                                                                                      gcACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCg                                             <  6:195502/149‑1 (MQ=255)                                                                                                       cACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCCGACCCCGATTTTTTCGc                                            <  3:213400/149‑1 (MQ=255)                                                                                                            gTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTc                                              <  7:179254/142‑1 (MQ=255)                                                                                                            gTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTc                                              >  8:179254/1‑142 (MQ=255)                                                                                                            gTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTAt                                         >  6:228099/1‑147 (MQ=255)                                                                                                            gTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTATCa                                       <  2:38697/149‑1 (MQ=255)                                                                                                                ccATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTATCAGGaa                                   <  2:178420/149‑1 (MQ=255)                                                                                                                 cATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGAccc                                                        <  3:284171/127‑1 (MQ=255)                                                                                                                 cATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGAccc                                                        >  4:284171/1‑127 (MQ=255)                                                                                                                  aTCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTATCAGGAAtt                                 >  2:83047/1‑149 (MQ=255)                                                                                                                         gCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTATCAGGAATTTTCATCg                          >  2:325292/1‑149 (MQ=255)                                                                                                                         gCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTATCAGGAATTTTCATCg                          >  5:348024/1‑149 (MQ=255)                                                                                                                           aGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTATCAGGAATTTTCATCGCa                        >  5:53253/1‑149 (MQ=255)                                                                                                                            gCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTATCAGGAATTTTCATCGCAt                       <  4:207698/149‑1 (MQ=255)                                                                                                                               aTTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTATCAGGAATTTTCATCGCATGGa                    <  3:87492/149‑1 (MQ=255)                                                                                                                                ttCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCGCCGCAGGGCTGACCCCGATTTTTTCGCTATCAGGAATTTTCATCGCATGGa                    <  1:83047/149‑1 (MQ=255)                                                                                                                                      gCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTATCAGGAATTTTCATCGCATGGATTACGaa             >  3:164797/1‑149 (MQ=255)                                                                                                                                          acgCTGGAGATATCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTATCAGGAATTTTCATCGCATGGATTACGAATTa          >  6:81653/3‑148 (MQ=255)                                                                                                                                            gCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTATCAGGAATTTTCATCGCATGGATTACGAATTACa        <  6:53253/148‑1 (MQ=255)                                                                                                                                            gCTGGAGATATCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTATCAGGAATTTTCATCGCATGGATTACGAATTACAg       <  5:81653/149‑1 (MQ=255)                                                                                                                                                 atatCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTATCAGGAATTTTCATCGCATGGATTACGAATTACAGGAAgg  >  7:274091/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TTCGCCAACGTGACTGCCTCCTGGCTGGGAACCCGGCCAAGGTCGTGGACTTTCACTTCAATGCCGGCC‑‑A‑CGAGCCTCGGCACGGTCGCCTGGTCGGCGCACTGGTGCCCATCGCCGGCAGCGATTCGCTGCAATGGCTGGATATCTCCAGCGTGATCGACCCCAGCAAGCCCGGCCTGAGCCAGTGGCTGGCCGAGTTCGACCGCATGCTCGATGGCTTCG‑CGCAGGGCTGACCCCGATTTTTTCGCTATCAGGAATTTTCATCGCATGGATTACGAATTACAGGAAGG  >  NC_002947/3510592‑3510881

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