Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	961,757	G→A	6.0%	R186H (CGT→CAT)	cmk →	cytidylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	961,757	0	G	A	6.0%	213.3 / 3.0	68	R186H (CGT→CAT)	cmk	cytidylate kinase
Reads supporting (aligned to +/- strand):  ref base G (28/35);  new base A (2/2);  total (31/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.70e-01

CCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAG  >  NC_000913/961614‑961900                                                                                                                                                |                                                                                                                                                cccTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGAGCGTAAc                                                                                                                                             <  6:332971/148‑1 (MQ=255) cccTGATGCACCAGTGAAAATTTTCATTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTTGCCGAGAACAAAGAACGCGACTACCGAGAACCTAAcc                                                                                                                                            >  6:335043/1‑149 (MQ=255)  ccTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCATAAcc                                                                                                                                            >  3:1843/1‑148 (MQ=255)     gATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGc                                                                                                                                        >  6:280962/1‑149 (MQ=255)      aTGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAg                                                                                                                                         <  1:36447/147‑1 (MQ=255)                        tCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTcc                                                                                                                      <  4:102730/148‑1 (MQ=255)                        tCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCg                                                                                                                     >  2:256140/1‑149 (MQ=255)                         ccTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCg                                                                                                                     <  2:8491/148‑1 (MQ=255)                          cTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCgg                                                                                                                    <  5:280962/148‑1 (MQ=255)                             gACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCa                                                                                                                  <  2:107429/147‑1 (MQ=255)                             gACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGc                                                                                                                >  5:117922/1‑149 (MQ=255)                                gCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGa                                                                                                             >  8:211023/1‑149 (MQ=255)                                gCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGa                                                                                                             <  1:256140/149‑1 (MQ=255)                                gCCTCCTAGGAAGAACGTGCGCAACGCCGCATGCTACAGTTGCAGTAGAAGGGATTTAGTGTTAAAATTGAGCGCCTATTGGCCGAGATCAAAGAACGCGACGCCCGCGATCGTAACCTAGAGGTAGCGCCACTGGTTCCGGCAGCCGa                                                                                                             >  8:145506/1‑149 (MQ=255)                                   tcctcGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGc                                                                                                          >  4:325543/1‑149 (MQ=255)                                     ctcGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGACCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCtt                                                                                                        <  3:186012/149‑1 (MQ=255)                                      tcGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCtt                                                                                                        >  8:322651/1‑148 (MQ=255)                                        ggAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTa                                                                                                      <  6:230266/148‑1 (MQ=255)                                        ggAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTa                                                                                                      <  3:313156/148‑1 (MQ=255)                                        ggAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTa                                                                                                      <  8:3923/148‑1 (MQ=255)                                         gaagaaCGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAgt                                                                                                    >  7:30303/1‑149 (MQ=255)                                            gaaCGTGCTAATCGCCGTATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATAAAATAACGCGACGACCGGGTTGGTAACAGAGCGGTAGCGCCACTGGTTCCGGCAGCCGTTGCTTTAGTGtt                                                                                                 <  7:145506/149‑1 (MQ=255)                                                  gcgcATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAgcgc                                                                                                                                >  6:80982/1‑112 (MQ=255)                                                  gcgcATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAgcgc                                                                                                                                <  5:80982/112‑1 (MQ=255)                                                   cgcATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTcc                                                                                          <  8:301610/149‑1 (MQ=255)                                                   cgcATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCATAACCGAGCGGTAGCGCCACTGGTTCCGGCATCCGATGCTTTAGTGTTGGATTc                                                                                           <  4:1843/148‑1 (MQ=255)                                                    gcATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGACCGCGACGACCGCGATCGTACCCGAGCGGTAGCGCCACTTGTTCCGGCAGCCGATGCTTTAGTGTTGGATTcca                                                                                         >  2:290656/1‑149 (MQ=255)                                                    gcATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTcca                                                                                         <  2:156722/149‑1 (MQ=255)                                                      aTCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTccac                                                                                        >  8:115115/1‑148 (MQ=255)                                                       tCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGCGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTccacca                                                                                      <  2:104890/149‑1 (MQ=255)                                                           cgcATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGc                                                                                                                >  7:62931/1‑119 (MQ=255)                                                           cgcATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGc                                                                                                                <  8:62931/119‑1 (MQ=255)                                                             cATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTaa                                                                                <  5:53740/149‑1 (MQ=255)                                                                  tACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGa                                                                                              >  6:147656/1‑130 (MQ=255)                                                                  tACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGa                                                                                              <  5:147656/130‑1 (MQ=255)                                                                  tACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCAtt                                                                           >  2:306493/1‑149 (MQ=255)                                                                  tACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCATTTTGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGAGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTAAAGCAtt                                                                           >  1:317697/1‑149 (MQ=255)                                                                  tACAGTTGCAGGAGAAGGGCTTTAGTGTTAAATTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGCCCGCGATCGTAACCGAGCGGTAGCGCCACGGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTGAAGCAtt                                                                           >  4:218604/1‑149 (MQ=255)                                                                     aGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAg                                                                        <  8:37748/149‑1 (MQ=255)                                                                     aGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAg                                                                        >  7:200000/1‑149 (MQ=255)                                                                         gCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCaa                                                                     <  7:322651/148‑1 (MQ=255)                                                                                         aGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAgcgc                                                     >  8:172484/1‑148 (MQ=255)                                                                                         aGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAgcgc                                                     >  8:167664/1‑148 (MQ=255)                                                                                          gtgtTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTa                                                   <  7:115115/149‑1 (MQ=255)                                                                                           tgtTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTa                                                   >  6:221213/1‑148 (MQ=255)                                                                                             ttAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCTGCAGCCGATGCTTTAGTGTTGGATTCAACCACCTAAAGCATTTAGCAAGTGATTGAAAAAGCTCTACa                                                 >  8:160078/1‑148 (MQ=255)                                                                                             ttAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACaa                                                >  3:187715/1‑149 (MQ=255)                                                                                             ttAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACaa                                                <  7:29222/149‑1 (MQ=255)                                                                                               aaCTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATa                                              >  7:156843/1‑149 (MQ=255)                                                                                                    tGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATAcgc                                           <  8:200000/147‑1 (MQ=255)                                                                                                     gAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATAcgcgc                                         <  6:117922/148‑1 (MQ=255)                                                                                                               tGGCCGAGATCAAAGACCGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCGTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAAtt                               <  7:332150/148‑1 (MQ=255)                                                                                                                ggCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTg                              <  1:306493/148‑1 (MQ=255)                                                                                                                 gCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTg                                                                          >  3:19966/1‑103 (MQ=255)                                                                                                                 gCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTg                                                                          <  4:19966/103‑1 (MQ=255)                                                                                                                     agaTCAAAGAACGCGACGACCGCGATCGTAACCGATCGGTAGCGCCACTGGTTCCGGCAGCCGATGCGTTAGTGTTGGATTCCACCACCTTAAGCAATGAGCAAGTGATTGAAAAAGCGCTACAATAAGCGCGCCAGAAATTGGCtctc                        <  1:220048/149‑1 (MQ=255)                                                                                                                       aTCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCg                       >  6:100030/1‑148 (MQ=255)                                                                                                                           aaGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCAt                    >  4:202261/1‑147 (MQ=255)                                                                                                                           aaGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATaa                  <  3:325543/149‑1 (MQ=255)                                                                                                                                 gcgACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACgcgcgc                                       >  4:152454/1‑122 (MQ=255)                                                                                                                                 gcgACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACgcgcgc                                       <  3:152454/122‑1 (MQ=255)                                                                                                                                 gcgACGACCGCGATCATAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGAcc            >  8:81638/1‑149 (MQ=255)                                                                                                                                     cgacCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAAt        <  7:160078/149‑1 (MQ=255)                                                                                                                                      gacCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAAtt       >  5:208452/1‑149 (MQ=255)                                                                                                                                        ccGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTg     <  7:172484/149‑1 (MQ=255)                                                                                                                                        ccGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTg     <  7:167664/149‑1 (MQ=255)                                                                                                                                          gcgATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCa   <  7:211023/149‑1 (MQ=255)                                                                                                                                           cgATCATAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAg  <  7:81638/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                CCCTGATGCACCAGTGAAAATTTTCCTTGACGCCTCCTCGGAAGAACGTGCGCATCGCCGCATGCTACAGTTGCAGGAGAAGGGCTTTAGTGTTAACTTTGAGCGCCTTTTGGCCGAGATCAAAGAACGCGACGACCGCGATCGTAACCGAGCGGTAGCGCCACTGGTTCCGGCAGCCGATGCTTTAGTGTTGGATTCCACCACCTTAAGCATTGAGCAAGTGATTGAAAAAGCGCTACAATACGCGCGCCAGAAATTGGCTCTCGCATAAGCGACCGAATTTGCAG  >  NC_000913/961614‑961900

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