Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,130,192	A→T	8.1%	intergenic (‑247/+57)	gcvH‑I ← / ← PP_0990	glycine cleavage system protein H/DoxX family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,130,192	0	A	T	8.1%	174.0 / 5.7	62	intergenic (‑247/+57)	gcvH‑I/PP_0990	glycine cleavage system protein H/DoxX family protein
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base T (2/3);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.96e-01

TTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGTGGG  >  NC_002947/1130046‑1130337                                                                                                                                                   |                                                                                                                                                  tttatCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTacgacg                                                                                                                                                 <  3:3586/149‑1 (MQ=255)       aaTAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCg                                                                                                                                            >  1:223251/1‑148 (MQ=255)         tAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGaaa                                                                                                                                         <  6:181480/149‑1 (MQ=255)           gATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGAAAAc                                                                                                                                       <  6:170192/149‑1 (MQ=255)           gATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGaaaa                                                                                                                                        <  8:194261/148‑1 (MQ=255)           gATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAAc                                                                                                                                       <  5:309215/149‑1 (MQ=255)               acttAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGAACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGAt                                                                                                                                    <  7:316891/145‑1 (MQ=255)                    aaTTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATGCAtgt                                                                                                                              <  6:13979/149‑1 (MQ=255)                                  ccGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTgg                                                                                                                >  4:24833/1‑149 (MQ=255)                                    gTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACAtg                                                                                                                               >  6:327759/1‑132 (MQ=255)                                    gTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACAtg                                                                                                                               <  5:327759/132‑1 (MQ=255)                                    gTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGc                                                                                                               >  4:145093/1‑148 (MQ=255)                                           ttGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGaaaa                                                                                                                                        <  8:204338/116‑1 (MQ=255)                                           ttGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGTCGTAGCGaaaa                                                                                                                                        >  7:204338/1‑116 (MQ=255)                                              tgtgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCTAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAAt                                                                                                    <  4:172691/149‑1 (MQ=255)                                              tgtgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAAt                                                                                                    >  3:228794/1‑149 (MQ=255)                                              tgtgGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAAt                                                                                                    <  4:228794/149‑1 (MQ=255)                                                tgGAATGATATCGTTACAGGTTCCTCTGAGAGGAATAGGGTTTGAGCATCAAACCCTTGTTCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATg                                                                                                   <  5:275278/148‑1 (MQ=255)                                                  gAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGcc                                                                                                >  3:36547/1‑149 (MQ=255)                                                  gAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGcc                                                                                                >  5:246667/1‑149 (MQ=255)                                                  gAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGcc                                                                                                >  1:329456/1‑149 (MQ=255)                                                  gAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGc                                                                                                 <  5:366189/148‑1 (MQ=255)                                                   aaTGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCt                                                                                               >  7:205542/1‑149 (MQ=255)                                                          tCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTc                                                                                         <  4:36547/148‑1 (MQ=255)                                                                cAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCa                                                                                   <  6:120607/148‑1 (MQ=255)                                                                   gTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGc                                                                                                               <  6:344212/117‑1 (MQ=255)                                                                   gTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGc                                                                                                               >  5:344212/1‑117 (MQ=255)                                                                   gTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATg                                                                                                   >  4:253677/1‑129 (MQ=255)                                                                   gTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATg                                                                                                   <  3:253677/129‑1 (MQ=255)                                                                       ctctGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTcagc                                                                           <  3:24833/149‑1 (MQ=255)                                                                          tGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTa                                                                                                                                               >  6:46895/1‑78 (MQ=255)                                                                          tGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTa                                                                                                                                               <  5:46895/78‑1 (MQ=255)                                                                           gagaGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAga                                                                        <  7:127389/148‑1 (MQ=255)                                                                           gagaGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAga                                                                        >  8:127389/1‑148 (MQ=255)                                                                               ggAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATc                                                                   <  8:205542/149‑1 (MQ=255)                                                                                gAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATc                                                                   >  5:292958/1‑148 (MQ=255)                                                                                 aaaaGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCa                                                                 <  7:46088/149‑1 (MQ=255)                                                                                  aaaGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGa                                                                                                             >  7:104123/1‑104 (MQ=255)                                                                                  aaaGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGa                                                                                                             <  8:104123/104‑1 (MQ=255)                                                                                        tttGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCTAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCg                                                                                            <  2:176193/115‑1 (MQ=255)                                                                                        tttGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCg                                                                                            >  1:176193/1‑115 (MQ=255)                                                                                            aGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAAcgc                                                       >  7:16829/1‑148 (MQ=255)                                                                                            aGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAAcgc                                                       <  8:16829/148‑1 (MQ=255)                                                                                                        ttGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCCGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCgg                                          >  2:32430/1‑149 (MQ=255)                                                                                                        ttGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCgg                                          >  2:33942/1‑149 (MQ=255)                                                                                                          gtgCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCggg                                         <  2:329456/148‑1 (MQ=255)                                                                                                          gtgCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCgg                                          <  1:335311/147‑1 (MQ=255)                                                                                                          gtgCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCgg                                          >  2:335311/1‑147 (MQ=255)                                                                                                             cACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAAc                                     <  3:145093/149‑1 (MQ=255)                                                                                                                cGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGgc                                                 >  6:200917/1‑134 (MQ=255)                                                                                                                cGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGgc                                                 <  5:200917/134‑1 (MQ=255)                                                                                                                      cGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGcc                                                                                                >  8:208562/1‑81 (MQ=255)                                                                                                                      cGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGcc                                                                                                <  7:208562/81‑1 (MQ=255)                                                                                                                       gaggagGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGagcag                            >  4:310620/1‑148 (MQ=255)                                                                                                                            ggCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAagag                      <  3:310620/149‑1 (MQ=255)                                                                                                                                    aTCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGagca                             >  6:120651/1‑134 (MQ=255)                                                                                                                                    aTCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGagca                             <  5:120651/134‑1 (MQ=255)                                                                                                                                     tCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCTGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGcc               >  5:254274/1‑147 (MQ=255)                                                                                                                                       gTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcagg           <  2:223251/149‑1 (MQ=255)                                                                                                                                        ttttCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggt          >  8:304579/1‑149 (MQ=255)                                                                                                                                        ttttCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggt          >  4:80203/1‑149 (MQ=255)                                                                                                                                            cGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggtcag       <  7:304579/148‑1 (MQ=255)                                                                                                                                             gCTACGTCGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCcaggtcaggt     >  2:147856/1‑149 (MQ=255)                                                                                                                                                 cgacgTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAggtggg  <  6:292958/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TTTATCAATAGATTACGTAAATTTACGGAGAAGCCGTATTTATTGTGTGGAATGATATCGTTACAGGTGCCTCTGAGAGGAAAAGGGTTTGAGCATCAAACCCTTGTGCACCGACGGCGAGGAGGCTGTGTATCGTTTTCGCTACGACGTAGCGAAAACGATACATGTGGTGCCTGGACTGGCGAGGTGCCAATGGCCTTCGGTCATCTCACTTTCAGCAGATGATCCAGCGAAAACGCGCCGGCGCCTCGGGCAACCAGGAGCAGCAAGAGCCCTGCCCAGGTCAGGTGGG  >  NC_002947/1130046‑1130337

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