Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	440,351	T→A	18.0%	intergenic (+181/‑121)	dmdC_1 → / → CJ019_00398	3‑methylmercaptopropionyl‑CoA dehydrogenase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	440,351	0	T	A	18.0%	121.8 / 4.7	67	intergenic (+181/‑121)	dmdC_1/CJ019_00398	3‑methylmercaptopropionyl‑CoA dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (24/27);  new base A (3/9);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCTTGTTTCGTCTTAA  >  ACB122/440229‑440485                                                                                                                            |                                                                                                                                       tCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTTCCCGCGTATCGGCCGGCCCTGCTGGCCACTCTCCATCCGCCGTGCCCCGTTGGCCGCGCGTCCC‑TCTTCGCGGGTATACACCGTCCATCAAGGCTcgac                                                                                                                            >  4:258983/1‑133 (MQ=255)            gTGGTCCCTCGTAGGAGCGGGTTTACCAGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCGACAAGGCTCGCGTGGCGGCCg                                                                                                                   >  8:452886/1‑133 (MQ=255)            gTGGTCACTCGTACGAGCCGGTTTACCCGCGAAGCGGCCGGCCCTGCTGGCCAATCAAAATCAGACGTTCACCGATGCCAGCCCAACCG‑TATAAGCGGGTAAACCCCCACCAACAAAGCAAGCCTGGCCGCCCGa                                                                                                                 >  4:239500/1‑135 (MQ=12)              ggTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTTAACCCGCGACTACAAGGCTCGCGTGGCGGCCGGAAg                                                                                                               >  7:258242/1‑135 (MQ=255)                               gggttAACCGCGAAGAGGCCGGCCCTGCTGGCCAATCTGCATCAGCCGGTCACCGGTGCCAGCGCTACCC‑TATTCGCGGGGAAACCCGCGCCGACAAAGCGAGCGGGGGCGCCGCGAGGTGGGTGATCCAACCg                                                                                               >  6:203540/4‑134 (MQ=12)                                    ccccGCTAAGAGGCCGGCCCTGCTGCCAAAACACCACCCCCCGACCACGGTGGCCAGCGCAACCT‑TATTCGCGGTAAAACCCGCTCCAACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAAtt                                                                                         <  2:92857/134‑1 (MQ=255)                                     cccGCTACGCGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGAAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAAtt                                                                                         <  1:306279/134‑1 (MQ=255)                                          taaGAGGCGGGCCCTGCTGGCCAATCTCCATCAGCCGTGCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAACCCCGCGCCTACAAGGCTCGCGTGGCAGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTg                                                                                   <  6:448146/134‑1 (MQ=37)                                          gAATAGCCCGGCCCTTCTGGCAAATCTCAATCAGCGTTTCCCCGTTGCCTGCGCTACCT‑TATTCGCGGTTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTg                                                                                   <  2:306533/135‑1 (MQ=37)                                            agagGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGtt                                                                                 >  8:187381/1‑135 (MQ=255)                                                   ggccCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCCTCCAACCGCAAATTCAACTGTTCCAGTGc                                                                          >  8:152698/1‑135 (MQ=255)                                                    gccCTGCTGGCCAATCTCCAGCAGCCGTGCCCCGTTCCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCATACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCa                                                                         <  4:415077/135‑1 (MQ=37)                                                    gccCTGCTGGCAAATCTCCATCACCCGTTCACCGTGGCCCGCTCTACCT‑TATTCGCGGTAAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCa                                                                         <  1:97902/135‑1 (MQ=38)                                                    gccCTGCGGGCCAATCTCAAGCAGCCGAGCACCTTTCCCCGCGCTACCG‑TATTCGCGGGAAAACCCGCTCCAACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCa                                                                         <  7:299946/135‑1 (MQ=255)                                                     cccTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGACAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCAc                                                                        >  6:114234/1‑135 (MQ=255)                                                       ctgctgGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCGCCGACCAGGCGCGCGTGGCGGCCGGCAGTTGGGGGATGCCACCGCCAATTCAACTGGTCCCGTggcctg                                                                      >  2:458044/1‑130 (MQ=18)                                                       ctgctgGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCCCTCCATCAAGGCGCGCGTGGCCGCCGGGATGTAGTTCATCCAACCGCAAATTGAAATGATCCATAGCACTg                                                                      >  2:107969/1‑135 (MQ=18)                                                         gctgGCCAAGCTCCATCAGCCGTTCACGGTTGCCCGCGCTACGG‑TATTCGCGGGAAAACCCGCTCAAACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGcc                                                                    <  1:377640/135‑1 (MQ=255)                                                           tgGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGcc                                                                    <  7:187381/133‑1 (MQ=255)                                                           tgGCAAATCTCCTTCAGCTGTTCCCCTTTGCCAGCGCTACCG‑TTTTGGCGGTTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGa                                                                  <  3:205527/135‑1 (MQ=255)                                                           gggCAAATCACCACCACCCGGTCCCGGTTGCCCGCGCACCCG‑TATTCGCGGGTAACCCCGCACCAACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGa                                                                  <  6:42657/134‑1 (MQ=255)                                                             gCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGAt                                                                 >  7:443236/1‑134 (MQ=255)                                                                aaTCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAg                                                             <  8:372914/135‑1 (MQ=255)                                                                     ccATCAGCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAgg                                                        >  1:442208/1‑135 (MQ=255)                                                                       accAGCCGTTCCCCGTGCCCAGCGCTCCCG‑TATGCGCGGGAAAACCCGCTCCAACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGc                                                      <  5:162812/133‑1 (MQ=255)                                                                       aTCCGCCGTCCCCCGTTGCCCCCGCTACCG‑TATCCGCGGTAAAACCCGCACCAACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGc                                                      <  1:327886/135‑1 (MQ=255)                                                                       aTCAGCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGc                                                      >  6:8921/1‑135 (MQ=255)                                                                       aTCAGCCGTGCACCGTTGCCAGCGCTCCCG‑TATTCGCGGGAAAACCCGCGCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAgg                                                        <  6:383839/133‑1 (MQ=255)                                                                       aTCAGCCGGTCACCGTGGCCAGCGCTACCG‑TATTCGCGGTAAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGc                                                      <  6:103012/135‑1 (MQ=255)                                                                        cccgCCGTCCCCCGTGCCCAGCGC‑ACCGTTATTCGCGGGAAACCCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGca                                                     <  1:167018/132‑1 (MQ=255)                                                                          aGCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCCAATTCCACTGGTGCCTTGCACTGCCGATGAAGTGAGGGgaca                                                   >  2:84612/1‑135 (MQ=37)                                                                           gCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACaa                                                  <  5:8921/135‑1 (MQ=255)                                                                            ccGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCaaa                                                                                           >  8:62073/1‑93 (MQ=255)                                                                            ccGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCaaa                                                                                           <  7:62073/93‑1 (MQ=255)                                                                                tCACCGTTCCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCa                                             <  4:239367/135‑1 (MQ=255)                                                                                  aCCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACaa                                                  <  3:92850/128‑1 (MQ=255)                                                                                  aCCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACaa                                                  >  4:92850/1‑128 (MQ=255)                                                                                   ccGTTGCCAGCGCTACCG‑TATTCGCGGGTTAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCCCAATGCCatta                                          >  4:51381/1‑135 (MQ=255)                                                                                   ccGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCatta                                          >  8:429711/1‑135 (MQ=255)                                                                                     gTTGCCCGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCattatt                                        <  6:107494/135‑1 (MQ=255)                                                                                        gCCAGCGCTCCCG‑TATTCGCGGGAAAACCCGCTCCAACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGAt                                     <  1:365224/135‑1 (MQ=255)                                                                                         ccAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGGGGCCGGGAGTTCGTTCATCCAACCGCAAATTCAACTGGTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATc                                    >  6:341170/1‑135 (MQ=255)                                                                                         ccAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGAt                                     <  3:26417/134‑1 (MQ=255)                                                                                         ccAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATc                                    >  2:267/1‑135 (MQ=255)                                                                                              gcTACCG‑TATTGGCGGCTAACCCGTCGCATACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCACCCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCg                               <  7:440170/135‑1 (MQ=255)                                                                                               atACCG‑TATTCGCGGGTAAACCGGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCgg                              <  3:198737/134‑1 (MQ=255)                                                                                                 aCCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCACCCGCAAATTCAACTGTTCCAGTGCCCTGCCGATGAAGGGAGGGCACAATGCCATTGTTGATCTGCCGGGt                            >  3:324252/1‑135 (MQ=255)                                                                                                    g‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGtt                           >  7:416893/1‑133 (MQ=255)                                                                                                      tATTCGCGGGTAAACCCGCTCCTTCAAGGCTCGCGTGGCGGGCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGga                        >  3:251508/1‑135 (MQ=255)                                                                                                      tATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGga                        <  1:458044/135‑1 (MQ=255)                                                                                                              ggTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTAccc                >  5:334695/1‑135 (MQ=255)                                                                                                                tAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCtt              >  5:360057/1‑135 (MQ=255)                                                                                                                tAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGGTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGGTGGAGATTACCCtt              >  4:434149/1‑135 (MQ=255)                                                                                                                 aaccccGCTCCAACAAGGCTCGAGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCTTg             <  2:212959/132‑1 (MQ=255)                                                                                                                 aaaCCCGCTCGTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCTTg             <  6:358117/135‑1 (MQ=255)                                                                                                                 aaaCCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCTTg             <  2:442208/135‑1 (MQ=255)                                                                                                                 aaaCCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCTTg             <  5:114234/135‑1 (MQ=255)                                                                                                                 aaaCCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCTTg             >  3:205742/1‑135 (MQ=255)                                                                                                                 aaaCCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCTTg             >  1:374602/1‑135 (MQ=255)                                                                                                                 aaaCCCGCCCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCTTg             <  6:411819/135‑1 (MQ=255)                                                                                                                  accccGCGCCCACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCTTg             <  6:262960/132‑1 (MQ=255)                                                                                                                  aaCGCGCCCCAACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCTTg             <  2:121993/134‑1 (MQ=255)                                                                                                                  aaCCCGCCCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCTTg             <  7:277440/134‑1 (MQ=255)                                                                                                                   aCCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACaa                                                  <  5:364698/96‑1 (MQ=255)                                                                                                                   aCCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACaa                                                  >  6:364698/1‑96 (MQ=255)                                                                                                                   aCCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCTTGtt           >  4:434343/1‑135 (MQ=255)                                                                                                                         tCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCTTGTTTCGTCt     <  6:201619/135‑1 (MQ=255)                                                                                                                            tACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCTTGTTTCGTCTTaa  >  5:384224/1‑135 (MQ=255)                                                                                                                            |                                                                                                                                       TCAGGTTCAGCGTGGTCCCTCGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCTGCTGGCCAATCTCCATCAGCCGTTCACCGTTGCCAGCGCTACCG‑TATTCGCGGGTAAACCCGCTCCTACAAGGCTCGCGTGGCGGCCGGAAGTTCGTTCATCCAACCGCAAATTCAACTGTTCCAGTGCACTGCCGATGAAGTGAGGGCACAATGCCATTATTGATCTGCCGGGTTGGAGATTACCCTTGTTTCGTCTTAA  >  ACB122/440229‑440485

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