Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	3,003,619	G→A	6.0%	D311N (GAC→AAC)	bcsE →	Cyclic di‑GMP binding protein BcsE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,003,619	0	G	A	6.0%	222.8 / 2.2	66	D311N (GAC→AAC)	bcsE	Cyclic di‑GMP binding protein BcsE
Reads supporting (aligned to +/- strand):  ref base G (31/31);  new base A (2/2);  total (33/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.00e+00

GCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACAGCAACCTACCGCAGCCACCCCGCAGCGACACGGCGGCGCA  >  ACB122/3003486‑3003752                                                                                                                                      |                                                                                                                                      gcctggcctgACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTgc                                                                                                                                      >  4:154259/1‑134 (MQ=255) gcctggcctgACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGa                                                                                                                                      <  5:173477/135‑1 (MQ=255)   ctggcctgACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGAcc                                                                                                                                    <  4:231717/135‑1 (MQ=255)   ctggcctgACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCACTCCGAAGGGGTCGAGCAAGCGCTTGAcc                                                                                                                                    <  2:118886/135‑1 (MQ=255)           aCCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTc                                                                                                                            <  8:69511/135‑1 (MQ=255)           aCCCCGCTGCAGGCTTGCTCGAGTACGGTATTCCGCCTTGATGTCGACATTGCCTGCGTGGTTGGCGAAGTG‑‑‑‑‑‑CTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCt                                                                                                                       <  8:208802/134‑1 (MQ=255)             cccGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTc                                                                                                                            <  5:288147/133‑1 (MQ=255)             cccGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTc                                                                                                                            >  6:288147/1‑133 (MQ=255)             cccGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAg                                                                                                                         >  7:140518/1‑136 (MQ=255)                gCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCt                                                                                                                       <  5:306205/135‑1 (MQ=255)                   gCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTc                                                                                                                    >  2:29389/1‑135 (MQ=255)                    cAGGCTTGCTCGCGTACGGTATTTCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCa                                                                                                                   >  3:12398/1‑135 (MQ=255)                        cTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGa                                                                                                               <  2:361224/135‑1 (MQ=255)                          tGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCTCTTAACCACATCATCCATCTCTCATgg                                                                                                                >  2:120420/1‑132 (MQ=255)                          tGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTAACCACATCTTCCAGCTGTCATgg                                                                                                                >  2:115057/1‑132 (MQ=255)                                  tACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGTGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGtt                                                                                                     <  3:274801/135‑1 (MQ=255)                                  tACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGtt                                                                                                     <  2:331240/135‑1 (MQ=255)                                      gTATTCCGCCGTGATGGCGACATTGCCTGCGTGGGTGGCGAAGTGCTGTTCCTGTTCCTTTTCGCCTGCCGCTCCGAAGGGGGCGAGCAATCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTtcatca                                                                                                <  8:171227/136‑1 (MQ=255)                                          tCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGcc                                                                                             <  5:116113/135‑1 (MQ=255)                                            cgccgTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGcc                                                                                             >  7:236102/1‑133 (MQ=255)                                            cgccgTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGcc                                                                                             <  8:236102/133‑1 (MQ=255)                                            cgccgTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCa                                                                                            >  6:202905/1‑134 (MQ=255)                                            cgccgTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCa                                                                                            >  5:93591/1‑134 (MQ=255)                                              ccgTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCaggag                                                                                        >  8:71976/1‑136 (MQ=255)                                                  gATGGCGACATTGCCTGCGTGGTTGGAGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTa                                                                                     >  5:154669/1‑135 (MQ=255)                                                    tGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGcc                                                                                             >  3:293156/1‑125 (MQ=255)                                                    tGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGcc                                                                                             <  4:293156/125‑1 (MQ=255)                                                    tGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACt                                                                                   >  5:82797/1‑135 (MQ=255)                                                    tGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACt                                                                                   >  5:138163/1‑135 (MQ=255)                                                      gCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTACAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTCCTTg                                                                                 >  3:16018/1‑135 (MQ=255)                                                       cGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTg                                                                                 <  5:202905/134‑1 (MQ=255)                                                            ttGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCg                                                                           <  4:164929/135‑1 (MQ=255)                                                            ttGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCg                                                                           <  2:201812/135‑1 (MQ=255)                                                               ccTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTc                                                                         <  4:123372/134‑1 (MQ=255)                                                                  gCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACa                                                                     >  5:306216/1‑135 (MQ=255)                                                                     tGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGACGGGGTCGAGCAAGCGCTTAACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGcc                                                                  <  1:120420/135‑1 (MQ=255)                                                                     tGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGt                                                                          >  1:372531/1‑127 (MQ=255)                                                                     tGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGt                                                                          <  2:372531/127‑1 (MQ=255)                                                                     tGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTAACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGcc                                                                  <  1:115057/135‑1 (MQ=255)                                                                      ggTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCt                                                                 >  4:331156/1‑135 (MQ=255)                                                                      ggTTGGCGAAGTGCTGTTCCTGTTCCTGTGCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCt                                                                 <  3:65970/135‑1 (MQ=255)                                                                              aaGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTcc                                                        <  7:148229/136‑1 (MQ=255)                                                                              aaGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCGGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGGCGACAGCCTGGCCACTc                                                         <  1:224120/135‑1 (MQ=255)                                                                              aaGTGCGGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTc                                                         <  2:139400/135‑1 (MQ=255)                                                                                gTGCTGTTCCTGTTCCTGTTCGCCTGCAGCTCCGTAGTGGTAGAGCAAGCGCTTGACCACATCTACCATCTGTCATGGAATTAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCTTCGACAGCCTGGCCACTCCg                                                       <  5:14676/135‑1 (MQ=255)                                                                                       tcctgttcctgttcGCCTGCCGCTCCGAAGGGGGCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTc                                                 <  2:82029/134‑1 (MQ=255)                                                                                         ctgttcctgttcGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTc                                              >  8:191106/1‑135 (MQ=255)                                                                                         ctgttcctgttcGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTc                                              >  7:393078/1‑135 (MQ=255)                                                                                              cctgttcGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACa                                          >  3:62110/1‑134 (MQ=255)                                                                                                 gttcGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACAGCaa                                      >  5:66356/1‑135 (MQ=255)                                                                                                   tcGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTc                                                         <  6:7425/114‑1 (MQ=255)                                                                                                   tcGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTc                                                         >  5:7425/1‑114 (MQ=255)                                                                                                   tcGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACAGCAAc                                     >  1:95825/1‑134 (MQ=255)                                                                                                     gCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCAt                                                                                                                  <  3:138108/55‑1 (MQ=255)                                                                                                     gCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCAt                                                                                                                  >  4:138108/1‑55 (MQ=255)                                                                                                     gCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACAGCAACCTa                                  >  3:132314/1‑135 (MQ=255)                                                                                                     gCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACAGCAACCTAc                                 >  7:52284/1‑136 (MQ=255)                                                                                                           cGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACAGCAACCTACCGCAg                            <  6:82797/135‑1 (MQ=255)                                                                                                             cTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACAGCAACCTACCGCAGcc                          >  1:139050/1‑135 (MQ=255)                                                                                                             cTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACAGCAACCTACCGCAGcc                          <  4:62110/135‑1 (MQ=255)                                                                                                             cTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCAGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACAGCAACCTACCGCAGCCa                         >  7:311353/1‑136 (MQ=255)                                                                                                                cGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACAGCAACCTACCGCAGCCAcc                       <  6:183489/135‑1 (MQ=255)                                                                                                                           aGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACAGCAACCTACCGCAGCCACCCCGCAGcccca            >  3:31334/1‑131 (MQ=255)                                                                                                                           aGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACAGCAACCTACCGCAGCCACCCCGCAGCGacac           >  8:326959/1‑136 (MQ=255)                                                                                                                           aGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACAGCAACCTACCGCAGCCACCCCGCAGCGaca            <  4:44637/135‑1 (MQ=255)                                                                                                                                    ttGACCACATCTTCCAGCTGTCATGGAAAGCGCTGTTCATCAGCCAGGAGGTCCTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTGCCTCGCCAGCCACCTACCGCAGCCCCCCCGCAGCGCCACGGCGGCGCa  >  7:298338/1‑136 (MQ=255)                                                                                                                                    ttGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACAGCAACCTACCGCAGCCACCCCGCAGCGACACGgcggcg    >  6:280048/1‑134 (MQ=255)                                                                                                                                      |                                                                                                                                      GCCTGGCCTGACCCCGCTGCAGGCTTGCTCGCGTACGGTATTCCGCCGTGATGGCGACATTGCCTGCGTGGTTGGCGAAGTGCTGTTCCTGTTCCTGTTCGCCTGCCGCTCCGAAGGGGTCGAGCAAGCGCTTGACCACATCTTCCAGCTGTCATGGAAAGAGCTGTTCATCAGCCAGGAGGTACTTGCCGGCGTCGACAGCCTGGCCACTCCGGCCTTCCTCGACAGCAACCTACCGCAGCCACCCCGCAGCGACACGGCGGCGCA  >  ACB122/3003486‑3003752

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