Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,368,671	A→C	7.7%	A13A (GCT→GCG)	greA ←	transcription elongation factor GreA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,368,671	0	A	C	7.7%	179.4 / 7.2	65	A13A (GCT→GCG)	greA	transcription elongation factor GreA
Reads supporting (aligned to +/- strand):  ref base A (31/29);  new base C (2/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-01

TTCGACCATGCCCTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAGTGGTG  >  NC_002947/5368527‑5368811                                                                                                                                                  |                                                                                                                                             ttCGACCATGCCCTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAgcgc                                                                                                                                          >  1:62683/1‑149 (MQ=255)       cATGCCCTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACGAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGa                                                                                                                                    >  5:225127/1‑149 (MQ=255)       cATGCCCTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGa                                                                                                                                    <  2:79094/149‑1 (MQ=255)            ccTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGt                                                                                                                                >  6:148662/1‑148 (MQ=255)            ccTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGt                                                                                                                                >  2:24519/1‑148 (MQ=255)             cTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCTCGCGCCCTGGaagaaga                                                                                                                              <  1:262865/149‑7 (MQ=255)              tGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCAt                                                                                                                             >  8:22878/1‑149 (MQ=255)                      aCGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTAtt                                                                                                                     >  1:299078/1‑149 (MQ=255)                        gggCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTg                                                                                                                   <  7:59597/149‑1 (MQ=255)                          gcggcATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCTCGCGCCCTg                                                                                                                                      >  4:16890/1‑128 (MQ=255)                          gcggcATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCTCGCGCCCTg                                                                                                                                      <  3:16890/128‑1 (MQ=255)                            ggcATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTaa                                                                                                               <  2:62683/149‑1 (MQ=255)                                     ttCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGc                                                                                                      >  1:65900/1‑149 (MQ=255)                                       cGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGGGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCt                                                                                                     <  2:65900/148‑1 (MQ=255)                                        ggCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCttt                                                                                                   <  2:299078/149‑1 (MQ=255)                                         gCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGc                                                                                                      >  4:68557/1‑145 (MQ=255)                                         gCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGc                                                                                                      <  3:68557/145‑1 (MQ=255)                                         gCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCCGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCttt                                                                                                   >  6:213608/1‑148 (MQ=255)                                                ccTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGGTTTCAGTCCtg                                                                                           >  3:59267/1‑149 (MQ=255)                                                ccTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCt                                                                                            <  3:40236/148‑1 (MQ=255)                                                 cTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCtg                                                                                           >  6:250850/1‑148 (MQ=255)                                                  ttGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCa                                                                                                 <  3:256253/141‑1 (MQ=255)                                                  ttGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCa                                                                                                 >  4:256253/1‑141 (MQ=255)                                                    gAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCtgca                                                                                         <  8:176587/147‑1 (MQ=255)                                                     aGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCtgcatgc                                                                                      <  1:145918/149‑1 (MQ=255)                                                       ggcACCCAGCTCACGCGCTTCACCAATGGCCTTTGTCAGGGTCGGGCGGTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACTGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCGtgcatgca                                                                                     <  7:129397/146‑1 (MQ=255)                                                        tCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGa                                                                                   <  5:250850/149‑1 (MQ=255)                                                            ccAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCt                                                                               <  3:163318/149‑1 (MQ=255)                                                             cAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCttt                                                                                                   <  8:1636/128‑1 (MQ=255)                                                             cAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCttt                                                                                                   >  7:1636/1‑128 (MQ=255)                                                                 tCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTCGGGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGc                                                                             >  1:61957/1‑147 (MQ=255)                                                                  cACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGc                                                                         >  6:192401/1‑149 (MQ=255)                                                                   aCGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGc                                                                         >  6:196029/1‑148 (MQ=255)                                                                    cgcgCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGcga                                                                       <  6:225127/149‑1 (MQ=255)                                                                           cACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACCGGTGCTCTTCTTCCCGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACg                                                                 <  8:69534/148‑1 (MQ=255)                                                                           cACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGc                                                                             <  6:35126/136‑1 (MQ=255)                                                                           cACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGc                                                                             >  5:35126/1‑136 (MQ=255)                                                                            aCCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGt                                                               >  3:261996/1‑149 (MQ=255)                                                                              cAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCt                                                                                                           <  7:24915/103‑1 (MQ=255)                                                                              cAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCt                                                                                                           >  8:24915/1‑103 (MQ=255)                                                                                    ccTGACTCAGGCGCGGTCGCTC‑GTTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGTGCGCGAGCTCCCTGGACGGTCATCGGGTATTTGTTAATGCTAATTCTTTCAGTCCTTCATGCAGATCCTGCAAGAGACGAACGGTCTTTTCCg                                                       <  8:104709/149‑1 (MQ=255)                                                                                    ccTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCtgcatgc                                                                                      <  2:260781/118‑1 (MQ=255)                                                                                    ccTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCtgcatgc                                                                                      >  1:260781/1‑118 (MQ=255)                                                                                     cTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCg                                                       >  3:113411/1‑148 (MQ=255)                                                                                       gACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGAc                                                    >  5:194836/1‑149 (MQ=255)                                                                                        aCTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGAcc                                                   >  1:235214/1‑149 (MQ=255)                                                                                          tCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGcccac                                                 >  8:180965/1‑147 (MQ=255)                                                                                            aGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAAt                                               <  6:194836/149‑1 (MQ=255)                                                                                              gcgcGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAAttt                                             <  5:192401/149‑1 (MQ=255)                                                                                                    gcgcTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAgcgc                                       >  4:6369/1‑149 (MQ=255)                                                                                                      gcTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCtt                                     <  4:261996/149‑1 (MQ=255)                                                                                                        tC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCg                                   >  7:75039/1‑149 (MQ=255)                                                                                                             tCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACgg                                                                <  8:28132/116‑1 (MQ=255)                                                                                                             tCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACgg                                                                >  7:28132/1‑116 (MQ=255)                                                                                                                 gCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATgg                            <  3:6369/148‑1 (MQ=255)                                                                                                                    cAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCtt                         >  3:176126/1‑148 (MQ=255)                                                                                                                     aGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCtt                         <  2:61957/147‑1 (MQ=255)                                                                                                                            ggTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCg                <  7:22878/149‑1 (MQ=255)                                                                                                                                 tcttctTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATg           >  5:238073/1‑149 (MQ=255)                                                                                                                                     cttcCAGGGCGCGAGCGCCCTGGACGGTCATCTGTTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCATATACTGCAAGCTACGAACTGTCTTTTACGGACAGATTTTCAGCGCTTCGCAGATGGCTTCACAAGACGCAATGGTAg       >  2:220705/1‑149 (MQ=255)                                                                                                                                       tcCAGGGCGCTCGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAg                                          <  3:273549/112‑1 (MQ=255)                                                                                                                                       tcCAGGGCGCTCGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAg                                          >  4:273549/1‑112 (MQ=255)                                                                                                                                          aGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAgtggtg  <  7:261534/149‑1 (MQ=255)                                                                                                                                           gggCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAAt            >  3:65627/1‑138 (MQ=255)                                                                                                                                           gggCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAAt            <  4:65627/138‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                             TTCGACCATGCCCTGCTCCTCACGGGCGGCATGGTATTCGGCGTTTTCCTTGAGGTCACCCAGCTCACGCGCTTCACCAATGGCCTGACTCAGGCGCGGGCGCTC‑GGTCTTGCTCAGGAACAGGTGCTCTTCTTCCAGGGCGCGAGCGCCCTGGACGGTCATCGGGTATTTGGTAATGCTCATGCTTTCAGTCCTGCATGCAGATCCTGCAAGCGACGAACGGTCTTTTCCGGACCGAATTTCAGCGCTTCGCAGATGGCTTCACCAGCCGCAATGGTAGTGGTG  >  NC_002947/5368527‑5368811

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