Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,678,456	G→A	6.4%	S10S (TCG→TCA)	ydgH →	DUF1471 family periplasmic protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,678,456	0	G	A	6.4%	203.9 / 2.9	65	S10S (TCG→TCA)	ydgH	DUF1471 family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base G (26/34);  new base A (2/2);  total (28/37)
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.03e-01

GGAATGTGATATGAAAAATGACGCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTC  >  NC_000913/1678310‑1678598                                                                                                                                                   |                                                                                                                                               ggAATGTGATATGAAAAATGACGCAGACAGTTACACCGTTTAAATGTAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCgg                                                                                                                                               >  1:284560/1‑148 (MQ=255) ggAATGTGATATGAAAAATGACGCAGACAGTTACACCGTTTAAATGTAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGc                                                                                                                                              >  5:190834/1‑149 (MQ=255)     tgtgATATGAAAAATGACGCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTg                                                                                                                                          >  1:254989/1‑149 (MQ=255)     tgagATATGAAAAATGACGCAGACAGTTACACCGTTTAAATTCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGAGCAAGGCGAAGGATTATTTTTATGAACCTTAAGAACACCCTACTGGCGTCGGCACtc                                                                                                                                          >  1:254993/4‑148 (MQ=255)         atatGAAAAATGACGCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCttt                                                                                                                                      <  2:333000/149‑1 (MQ=255)         atatGAAAAATGACGCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCttt                                                                                                                                      <  8:244190/149‑1 (MQ=255)          tatGAAAAATGACGCAGACAGTTACACCGGTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTc                                                                                                                                     <  4:221313/149‑1 (MQ=255)            tGAAAAATGACGCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTg                                                                                                                                   <  1:67187/149‑1 (MQ=255)               aaaaTGACGCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAATGATTATTTTTATGAAGCTTAAGAACACCCTCCTTGCGTCCGCACTGCTTTCTGCt                                                                                                                                 <  6:123924/148‑1 (MQ=255)                       gCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATtc                                                                                                                        <  3:193015/149‑1 (MQ=255)                          gACAGTTACACCGTTTACATGCACTAATCAGCCACGTTTCTCGTTAATAACACTACCAGTACCTGGGTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTGCTGCTATGGCATTCTcc                                                                                                                     <  5:232392/149‑1 (MQ=255)                          gACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTcc                                                                                                                     >  8:111811/1‑149 (MQ=255)                               ttACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTaa                                                                                                                >  4:36992/1‑149 (MQ=255)                                  cacCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCAGCACTGCTTTCTGCTATGTCATTCTCCGTTAACg                                                                                                              >  3:61143/1‑148 (MQ=255)                                  cacCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCAGCACTGCTTTCTGCTATGGCATTCTCCGTTAACg                                                                                                              >  3:57000/1‑148 (MQ=255)                                   acCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACg                                                                                                              >  4:47034/1‑147 (MQ=255)                                   acCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACg                                                                                                              >  6:152479/1‑147 (MQ=255)                                    ccGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACgcag                                                                                                           <  2:107481/149‑1 (MQ=255)                                     cGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACgcagc                                                                                                          <  3:99646/149‑1 (MQ=255)                                     cGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCAATGGCATTCTCCGTTAACgcag                                                                                                           >  8:312745/1‑148 (MQ=255)                                     cGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACAATTACATGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGTCATTCTCCGTTAACgcagc                                                                                                          <  4:47906/149‑1 (MQ=255)                                       tttAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTTTGGCATTCTCCGTTAACGCAGCaa                                                                                                        <  2:284579/149‑1 (MQ=255)                                           aaTGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAg                                                                                                     <  2:156327/148‑1 (MQ=255)                                             tGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAAc                                                                                                  <  8:230087/149‑1 (MQ=255)                                              gCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACt                                                                                                 >  7:109411/1‑149 (MQ=255)                                                   aatCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGacac                                                                                            <  8:295980/149‑1 (MQ=255)                                                    atCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGacac                                                                                            <  3:47034/148‑1 (MQ=255)                                                              tttCTCGTTAATAACAATACCAGTACCTGGTTTGCGCATGGCGAAGGATTATTTTTATGAAGATTAAGAACACACTCCTGGAGTCGGCAATGCTTTATGATATGGAATTATCCGTTAAAGCAGCAACATAACTGACACCGGAGCAAGCg                                                                                 >  7:326708/1‑149 (MQ=255)                                                                 ctcGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGa                                                                                               >  6:180543/1‑132 (MQ=255)                                                                 ctcGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGa                                                                                               <  5:180543/132‑1 (MQ=255)                                                                          aaCAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTaaa                                                                     >  4:192191/1‑149 (MQ=255)                                                                          aaCAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACGGGTTTCTGCTTTGGCATTCTCCGTTAACGCAGCAACAGAAAGGAACCCGCAGCAAGGGGCAGCGGTTaaa                                                                     >  2:89004/1‑149 (MQ=255)                                                                           aCAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTaaa                                                                     <  6:190834/148‑1 (MQ=255)                                                                           aCAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTaaa                                                                     <  2:284560/148‑1 (MQ=255)                                                                               tACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCttt                                                                >  5:255507/1‑149 (MQ=255)                                                                                 ccAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCtttt                                                               <  3:192191/148‑1 (MQ=255)                                                                                 ccAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCtttt                                                               <  1:187509/148‑1 (MQ=255)                                                                                 ccAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGGCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCtttt                                                               <  5:227043/148‑1 (MQ=255)                                                                                 ccAGTACCTGGTTTGCGCAAGGCGAAGGATTAATTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCtttt                                                               <  5:100362/148‑1 (MQ=255)                                                                                     tACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCg                                                          <  8:66272/149‑1 (MQ=255)                                                                                     tACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCg                                                          >  1:163706/1‑149 (MQ=255)                                                                                       ccTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCgtg                                                        <  2:163706/149‑1 (MQ=255)                                                                                        cTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCttt                                                                >  7:64326/1‑140 (MQ=255)                                                                                        cTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCttt                                                                <  8:64326/140‑1 (MQ=255)                                                                                            tttGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTTCTTTCTGCTATGGAATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGtt                                                   <  5:99015/149‑2 (MQ=255)                                                                                               gcgcAAGGCGAAGGAATATTTTTATGAAGCTTAAGAACACACTCCTGGCGTAGGCACTGCTTTCAGCTATATAATTCTCCCTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAAACTTTTGACCGTGTAGAGGt                                                 >  2:7704/1‑148 (MQ=255)                                                                                                   aaGGCGTAGGATTATTTTTATGAAGATTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCg                                            <  4:271177/149‑1 (MQ=255)                                                                                                     ggCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCCGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGt                                                 <  5:272942/142‑1 (MQ=255)                                                                                                     ggCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCCGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGt                                                 >  6:272942/1‑142 (MQ=255)                                                                                                        gAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGt                                       >  5:170958/1‑149 (MQ=255)                                                                                                         aaGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGtt                                      <  5:152479/149‑1 (MQ=255)                                                                                                          aGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCAGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGttt                                     <  4:61143/149‑1 (MQ=255)                                                                                                           ggATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCAGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGttt                                     <  4:57000/148‑1 (MQ=255)                                                                                                               tattTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATg                                >  3:296975/1‑149 (MQ=255)                                                                                                               tattTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATg                                >  7:3592/1‑149 (MQ=255)                                                                                                                attTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATg                                <  7:111811/148‑1 (MQ=255)                                                                                                                 tttttATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGc                               >  3:272312/1‑148 (MQ=255)                                                                                                                      aTGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTgg                         >  1:161837/1‑149 (MQ=255)                                                                                                                       tGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTgg                         <  6:255507/148‑1 (MQ=255)                                                                                                                               taGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGt                >  2:520/2‑149 (MQ=255)                                                                                                                                   acacCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGaaa            <  3:36992/149‑1 (MQ=255)                                                                                                                                           ctGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGa                                                                                        >  6:43701/1‑65 (MQ=255)                                                                                                                                           ctGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGa                                                                                        <  5:43701/65‑1 (MQ=255)                                                                                                                                           ctGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTtc    <  6:238544/149‑1 (MQ=255)                                                                                                                                              gCGTCGTCACTGCTTTCAGCTATGTCATTCTCCGTTAACGCAGCAACAGAACTGACACCTGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTtctc  <  7:308489/148‑1 (MQ=255)                                                                                                                                              gCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTtctc  >  6:50224/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                               GGAATGTGATATGAAAAATGACGCAGACAGTTACACCGTTTAAATGCAATAATCAGCCACGTTTCTCGTTAATAACAATACCAGTACCTGGTTTGCGCAAGGCGAAGGATTATTTTTATGAAGCTTAAGAACACCCTCCTGGCGTCGGCACTGCTTTCTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTC  >  NC_000913/1678310‑1678598

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