Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	26,245	A→T	8.0%	Q122L (CAG→CTG)	PP_0021 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	26,245	0	A	T	8.0%	182.9 / 5.4	62	Q122L (CAG→CTG)	PP_0021	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base T (3/2);  total (32/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 = 9.75e-01

CTAAGAGCATTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATGGATCTTTGC  >  NC_002947/26101‑26392                                                                                                                                                 |                                                                                                                                                    ctaAGACCATTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGAAACGATTGACCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTc                                                                                                                                                 <  8:284144/149‑1 (MQ=255)     gagCATTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATc                                                                                                                                             <  4:112710/149‑1 (MQ=255)        cATTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATcgcg                                                                                                                                          <  4:2440/149‑1 (MQ=255)         aTTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGt                                                                                                                                         <  6:235597/149‑1 (MQ=255)           tCGCGGCATCGCTCCATTTGTTATTGAGGCCCGTAACTTTTTCCCGTACTTCCAGCCTCGGGAGATCAAGACTACAGAACGTATATCCATGTTTGACTTGGTCCCACGGATGATTCTGAAGCTATTGATCGTACat                                                                                                                                                    <  4:220723/136‑2 (MQ=255)           tCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAg                                                                                                                                                    >  3:220723/1‑136 (MQ=255)                  aTCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACAt                                                                                                                                <  6:237737/149‑1 (MQ=255)                     gCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCt                                                                                                                             >  6:52105/1‑149 (MQ=255)                       tCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACTGTTCGATCGCGTGGTTTACATCCTg                                                                                                                            >  4:303061/1‑148 (MQ=255)                         gaTTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTGGCCAGCACTAGAAAACAGAACGAAATACCTGGTTTGACGGGGTCCCACGGAAGATGCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGga                                                                                                                          <  3:239207/147‑1 (MQ=255)                         cATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGga                                                                                                                          <  5:52105/148‑1 (MQ=255)                                  gTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGAcc                                                                                                                <  3:121313/149‑1 (MQ=255)                                      ggCCCGTAACGATTGCCCTTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAAc                                                                                                            <  7:215615/149‑1 (MQ=255)                                      ggCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAAc                                                                                                            <  8:289145/149‑1 (MQ=255)                                         ccGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTgg                                                                                                         >  3:66112/1‑149 (MQ=255)                                         ccGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTgg                                                                                                         >  1:283811/1‑149 (MQ=255)                                                gATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCACACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTc                                                                                                  >  7:166063/1‑149 (MQ=255)                                                            ttCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAag                                                                                                                     <  5:74455/118‑1 (MQ=255)                                                            ttCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAag                                                                                                                     >  6:74455/1‑118 (MQ=255)                                                            ttCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGtt                                                                                      >  5:196252/1‑149 (MQ=255)                                                             tCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGt                                                                                       <  8:20664/147‑1 (MQ=255)                                                             tCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGt                                                                                       >  7:20664/1‑147 (MQ=255)                                                              ccAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTgt                                                                                                                           >  5:319832/1‑109 (MQ=255)                                                              ccAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTgg                                                                                                                           <  6:319832/110‑1 (MQ=255)                                                                aGCCTCGCCAGCACAAGAAAACTGAACGAAAAACCAGGTTTGCCGGCGTCCCACGGAAGATTCTCAAGCTAATGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTgaga                                                                                  <  8:267812/149‑1 (MQ=255)                                                                              aaGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCg                                                                    >  5:65121/1‑149 (MQ=255)                                                                              aaGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCg                                                                    >  3:31283/1‑149 (MQ=255)                                                                                gAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGtt                                                                  <  4:66112/149‑1 (MQ=255)                                                                                gAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGt                                                                   >  8:238309/1‑148 (MQ=255)                                                                                gAAAACAGAACGAAAAACCAGGTTTGACCGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCGTTACGAGAACGCGTTGATAGGCATGCGGCCGCGt                                                                   >  4:244074/1‑148 (MQ=255)                                                                                      aGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAgg                                                            <  8:166063/149‑1 (MQ=255)                                                                                      aGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAgg                                                            <  4:31283/149‑1 (MQ=255)                                                                                       gAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAggg                                                           <  2:104418/149‑1 (MQ=255)                                                                                        aaCGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGga                                                          <  6:196252/149‑1 (MQ=255)                                                                                           gAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGagag                                                                                 >  6:304874/1‑123 (MQ=255)                                                                                           gAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGagag                                                                                 <  5:304874/123‑1 (MQ=255)                                                                                             aaaaCCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACAACCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCg                                                     >  3:3953/1‑149 (MQ=255)                                                                                               aaCCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACTGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTg                                                   <  3:303061/149‑1 (MQ=255)                                                                                                 ccAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGc                                                 <  6:65121/149‑1 (MQ=255)                                                                                                            ggCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATAc                                      >  1:102776/1‑149 (MQ=255)                                                                                                                tCCCACGGAAGATTCTCAAGCTATTGATCGAACAGGTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCACCTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGgc                                  >  5:92056/1‑149 (MQ=255)                                                                                                                     cGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGAc                             <  2:283811/149‑1 (MQ=255)                                                                                                                          gATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTc                        >  3:102619/1‑149 (MQ=255)                                                                                                                               tcAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAAt                   >  1:116775/1‑149 (MQ=255)                                                                                                                                 aaGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAAtgt                 >  4:221269/1‑149 (MQ=255)                                                                                                                                 aaGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAAtgt                 <  2:102776/149‑1 (MQ=255)                                                                                                                                 aaGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAAtgt                 >  2:145941/1‑149 (MQ=255)                                                                                                                                     tATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCATCAGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATTCCGGCGCGACCCTTCGCAATGTGCCa             >  1:183723/1‑149 (MQ=255)                                                                                                                                      cgcgATCGAACTGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCAt            >  2:306740/4‑149 (MQ=255)                                                                                                                                      cgcgATCGAACTGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCAt            <  5:281341/146‑1 (MQ=255)                                                                                                                                      cgcgATCGAACTGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCAt            >  6:281341/4‑149 (MQ=255)                                                                                                                                       ttGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATg           >  3:67269/1‑149 (MQ=255)                                                                                                                                          aTCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATGGAt        >  6:150378/1‑149 (MQ=255)                                                                                                                                          aTCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATGGAt        >  1:249973/1‑149 (MQ=255)                                                                                                                                             gAACAGTTCGATCGCGTGGTTTTCATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATGGATCt      <  7:238309/148‑1 (MQ=255)                                                                                                                                             gAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGTGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCAt            <  5:87845/142‑1 (MQ=255)                                                                                                                                             gAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGTGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCAATACAGGCGCGACGCTTCGCAATGAGCCAt            >  6:87845/1‑142 (MQ=255)                                                                                                                                             gAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAAtgt                 >  7:230700/1‑137 (MQ=255)                                                                                                                                             gAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAAtgt                 <  8:230700/137‑1 (MQ=255)                                                                                                                                             gAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATg           <  4:300390/143‑1 (MQ=255)                                                                                                                                             gAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATg           >  3:300390/1‑143 (MQ=255)                                                                                                                                               aCAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGTGAGAGCGTGTCTACCTCTATACCGGCGCGACGCTTCGCAATGAGCAATGGAGCTTTg   >  2:249560/1‑149 (MQ=255)                                                                                                                                                cAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATGGATCTTTGc  <  5:150378/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CTAAGAGCATTCGCGGCATCGCTCCATTTGATAGTGAGGCCCGTAACGATTGCCCGTACTTCCAGCCTCGCCAGCACAAGAAAACAGAACGAAAAACCAGGTTTGACGGCGTCCCACGGAAGATTCTCAAGCTATTGATCGAACAGTTCGATCGCGTGGTTTACATCCTGGAGAAAGAGACCCAACTGGTCCTTTCCGAGAACGCGTTGAGAGGCATGCTGCAGCGTTACAAGGGAGAGCGTGGCTACCTCTATACCGGCGCGACGCTTCGCAATGTGCCATGGATCTTTGC  >  NC_002947/26101‑26392

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