Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,439,556	G→A	6.0%	intergenic (‑294/‑31)	ytfJ ← / → ytfK	putative transcriptional regulator/DUF1107 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,439,556	0	G	A	6.0%	264.8 / 4.1	82	intergenic (‑294/‑31)	ytfJ/ytfK	putative transcriptional regulator/DUF1107 family protein
Reads supporting (aligned to +/- strand):  ref base G (35/42);  new base A (3/2);  total (38/44)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GACCGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAATTCGATAAGGGTAAGATT  >  NC_000913/4439410‑4439703                                                                                                                                                   |                                                                                                                                                    gACCGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGc                                                                                                                                                    <  5:309623/148‑1 (MQ=255)  aCCGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGca                                                                                                                                                   <  2:350523/148‑1 (MQ=255)   cctcgcTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTATGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  <  4:99275/145‑1 (MQ=255)   ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATATTAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  <  6:63251/148‑1 (MQ=255)   ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  <  2:82899/148‑1 (MQ=255)   ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  <  7:101295/148‑1 (MQ=255)   ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  <  6:177931/148‑1 (MQ=255)   ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGca                                                                                                                                                   <  8:96672/147‑1 (MQ=255)   ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGca                                                                                                                                                   >  5:334424/1‑147 (MQ=255)   ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGca                                                                                                                                                   >  6:104261/1‑147 (MQ=255)   ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGca                                                                                                                                                   >  1:357851/1‑147 (MQ=255)   ccGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGCGGCACCGcac                                                                                                                                                  <  1:306863/148‑1 (MQ=255)           cGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTa                                                                                                                                           <  7:250495/147‑1 (MQ=255)            ggTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  >  6:166690/1‑139 (MQ=255)            ggTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  <  5:166690/139‑1 (MQ=255)            ggTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGGCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAAt                                                                                                                                         <  6:337877/148‑1 (MQ=255)              tAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  <  6:358564/137‑1 (MQ=255)              tAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGcac                                                                                                                                                  >  5:358564/1‑137 (MQ=255)                          tCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGgaga                                                                                                                            >  1:332702/1‑147 (MQ=255)                                 aatTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCCCAGGCGTAATACTCAGCAGGAGATAACAAt                                                                                                                     >  5:380186/1‑147 (MQ=255)                                      tAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAAt                                                                                                               <  1:137158/148‑1 (MQ=255)                                      tAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAAt                                                                                                               <  3:320401/148‑1 (MQ=255)                                            ttAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCACACAGGCGTAATACTcagc                                                                                                                                  >  6:301280/1‑123 (MQ=255)                                            ttAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGCTGGGCCGATTGTGGCACCACACAGGCGTAATACTcagc                                                                                                                                  <  5:301280/123‑1 (MQ=255)                                             tAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCa                                                                                                         <  7:206947/147‑1 (MQ=255)                                               gttAAAAAAAGTTATACTCGGGGGAAACATTGCCCGTATAGTCTATAGTCACTAAGCATTAAAATTTTCGCCTCATAATAGTTGGGCCGATTGTGGCAACGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAAc                                                                                                       <  5:45117/145‑1 (MQ=255)                                                    aaaaaGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACa                                                                                                 <  5:104261/148‑1 (MQ=255)                                                    aaaaaGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGACTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACa                                                                                                 <  2:332702/148‑1 (MQ=255)                                                                       aaaCATTGCCCGGATAGTCTATAGTCACTAAGCAATAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCg                                                                                >  8:198258/1‑146 (MQ=255)                                                                        aaCATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAAt                                                                                                               <  4:131184/114‑1 (MQ=255)                                                                        aaCATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAAt                                                                                                               >  3:131184/1‑114 (MQ=255)                                                                                  ggATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCTCAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAGCGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGa                                                                   <  6:380186/148‑1 (MQ=255)                                                                                    aTAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATc                                                                 >  8:154330/1‑148 (MQ=255)                                                                                      aGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCt                                                               >  1:187024/1‑148 (MQ=255)                                                                                         cTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCACACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGtt                                                            >  8:123853/1‑148 (MQ=255)                                                                                           ataGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTcc                                                          >  5:223235/1‑148 (MQ=255)                                                                                              gTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGAt                                                                  <  2:275788/137‑1 (MQ=255)                                                                                              gTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGAt                                                                  >  1:275788/1‑137 (MQ=255)                                                                                              gTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGt                                                        <  2:187024/147‑1 (MQ=255)                                                                                                cACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACtt                                                                                         >  2:351438/1‑112 (MQ=255)                                                                                                cACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACtt                                                                                         <  1:351438/112‑1 (MQ=255)                                                                                                    aaGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACaa                                                                                                >  8:265604/1‑101 (MQ=255)                                                                                                    aaGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACaa                                                                                                <  7:265604/101‑1 (MQ=255)                                                                                                    aaGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGt                                                 >  4:121233/1‑148 (MQ=255)                                                                                                    aaGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGt                                                 >  4:4337/1‑148 (MQ=255)                                                                                                       cATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTtat                                              <  3:4337/148‑1 (MQ=255)                                                                                                         ttAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATc                                                                 >  2:252702/1‑127 (MQ=255)                                                                                                         ttAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATc                                                                 <  1:252702/127‑1 (MQ=255)                                                                                                          tAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACa                                           >  8:310790/1‑148 (MQ=255)                                                                                                          tAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCATAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACa                                           >  3:331748/1‑148 (MQ=255)                                                                                                          tAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCACACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACa                                           >  1:207367/1‑148 (MQ=255)                                                                                                           aaaaTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAAGCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACAt                                          >  4:159906/1‑148 (MQ=255)                                                                                                           aaaaTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCg                                                         >  7:5804/1‑133 (MQ=255)                                                                                                           aaaaTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCg                                                         <  8:5804/133‑1 (MQ=255)                                                                                                           aaaaTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTAAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACAt                                          >  6:283902/1‑148 (MQ=255)                                                                                                            aaaTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGt                                                        <  8:33665/133‑1 (MQ=255)                                                                                                            aaaTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGt                                                        >  7:33665/1‑133 (MQ=255)                                                                                                              aTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCaa                                       <  2:357851/148‑1 (MQ=255)                                                                                                              aTTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCaa                                       <  6:334424/148‑1 (MQ=255)                                                                                                                       tCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTgg                              >  6:322487/1‑148 (MQ=255)                                                                                                                        cATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCa                                                                                       <  5:129837/90‑1 (MQ=255)                                                                                                                        cATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCa                                                                                       >  6:129837/1‑90 (MQ=255)                                                                                                                        cATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTtata                                             >  7:201153/1‑132 (MQ=255)                                                                                                                        cATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTtata                                             <  8:201153/132‑1 (MQ=255)                                                                                                                        cATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTg                               >  8:5731/1‑146 (MQ=255)                                                                                                                        cATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGGGTACGGTTATACATCAAGGACGTTGgc                             >  4:273729/1‑148 (MQ=255)                                                                                                                        cATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATACTGTTCCGTGGACGGTTATACATCAAGGACGTTGgc                             >  2:304492/1‑148 (MQ=255)                                                                                                                         ataataGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTgg                              <  7:5731/146‑1 (MQ=255)                                                                                                                          taataGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAAttt                                                                                                             <  1:260759/66‑1 (MQ=255)                                                                                                                          taataGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAAttt                                                                                                             >  2:260759/1‑66 (MQ=255)                                                                                                                            ataGTTGGGCCGATTGTGGCACCACACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCtt                         <  7:123853/148‑1 (MQ=255)                                                                                                                                ttGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGa                     <  7:198258/148‑1 (MQ=255)                                                                                                                                  gggCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAAGCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGaa                    <  3:159906/147‑1 (MQ=255)                                                                                                                                   ggCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAAtt                  <  7:154330/148‑1 (MQ=255)                                                                                                                                     ccGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAATTc                 >  8:61890/1‑147 (MQ=255)                                                                                                                                     ccGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAATTCg                <  7:310790/148‑1 (MQ=255)                                                                                                                                       gATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATCACAATGAAAATTTTCCAACGGTACACCCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAATTCGAt              <  5:239796/148‑1 (MQ=255)                                                                                                                                        aTTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTc                                                                                        <  8:270962/73‑1 (MQ=255)                                                                                                                                        aTTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTc                                                                                        >  7:270962/1‑73 (MQ=255)                                                                                                                                         ttGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAATTCGAt              >  1:286170/1‑146 (MQ=255)                                                                                                                                           gtgGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGgc                             >  7:318400/1‑129 (MQ=255)                                                                                                                                           gtgGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGgc                             <  8:318400/129‑1 (MQ=255)                                                                                                                                                   gCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAATTCGATAAGGGTAAGAtt  >  5:278353/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GACCGCGCTTCGGTAAAAAGATAATTCTGAATAATTGTAACCTTTAGGTAAAAAAAGTTATACGCGGTGGAAACATTGCCCGGATAGTCTATAGTCACTAAGCATTAAAATTTGCGCCTCATAATAGTTGGGCCGATTGTGGCACCGCACAGGCGTAATACTCAGCAGGAGATAACAATGAAAATTTTCCAACGGTACAACCCACTTCAGGTGGCGAAGTACGTAAAGATCCTGTTCCGTGGACGGTTATACATCAAGGACGTTGGCGCTTTTGAATTCGATAAGGGTAAGATT  >  NC_000913/4439410‑4439703

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