Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,286,116	C→G	6.3%	G419A (GGC→GCC)	acs ←	acetyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,286,116	0	C	G	6.3%	318.9 / 10.4	94	G419A (GGC→GCC)	acs	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (48/40);  new base G (3/3);  total (51/43)
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.43e-01

TGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTT  >  NC_000913/4285973‑4286263                                                                                                                                                |                                                                                                                                                     tGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGtt                                                                                                                                                 <  4:802784/149‑1 (MQ=255) tGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGtt                                                                                                                                                 >  2:697222/1‑149 (MQ=255)  gATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGttt                                                                                                                                                <  4:634849/149‑1 (MQ=255)  gATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGttt                                                                                                                                                <  4:414063/149‑1 (MQ=255)     tACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGGCGGTTTcat                                                                                                                                             <  2:980860/149‑3 (MQ=255)      aCCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGGCGGTTTcatg                                                                                                                                            <  2:471296/149‑4 (MQ=255)        cAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCt                                                                                                                                          <  1:397208/149‑1 (MQ=255)        cAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCt                                                                                                                                          <  4:116489/149‑1 (MQ=255)        cAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCt                                                                                                                                          >  2:781626/1‑149 (MQ=255)            cTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGcca                                                                                                                                      >  1:355958/1‑149 (MQ=255)              aCCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCTCCGGTTTCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGccacc                                                                                                                                    >  2:970591/1‑149 (MQ=255)                cTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAg                                                                                                                                  >  4:1003012/1‑149 (MQ=255)                 ttCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAgg                                                                                                                                 >  1:684312/1‑149 (MQ=255)                 ttCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAgg                                                                                                                                 >  4:522935/1‑149 (MQ=255)                  tCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGt                                                                                                                                <  3:120367/149‑1 (MQ=255)                    ggTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTAt                                                                                                                              <  1:697222/149‑1 (MQ=255)                     gTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCAGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATc                                                                                                                             >  4:976145/1‑149 (MQ=255)                        gCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGac                                                                                                                          <  4:863711/149‑1 (MQ=255)                         cccccTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGacc                                                                                                                         >  4:308559/1‑149 (MQ=255)                          ccccTCCAGCGGGTTACCTTCGTTATCGACCATCGCCGGTTGCACGCCGAAGAACGGACGTGTTACACAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCGTGAAACCGCCGGTTTCGGTCTGCCACCACGTATCGaccg                                                                                                                        >  2:69044/1‑148 (MQ=255)                          ccccTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGacca                                                                                                                        >  1:246181/1‑149 (MQ=255)                           cccTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGGCGGTTTcatga                                                                                                                                           <  4:854532/129‑5 (MQ=255)                           cccTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGGCGGTTTcatga                                                                                                                                           >  3:854532/1‑125 (MQ=255)                                      gTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTtct                                                                                                            >  1:976159/1‑149 (MQ=255)                                      gTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTtct                                                                                                            <  3:285566/149‑1 (MQ=255)                                         aCCTTCGTTATCGACCAGCGCGGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTGAGCGAGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCAACCACCGGACATTTCTCGt                                                                                                         <  4:807911/149‑1 (MQ=255)                                          ccTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGtt                                                                                                        <  2:702029/149‑1 (MQ=255)                                           cTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTg                                                                                                       >  3:237155/1‑149 (MQ=255)                                            ttCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGc                                                                                                      >  3:435566/1‑149 (MQ=255)                                            ttCGTTATCGACCAGCGCAGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGc                                                                                                      >  2:661702/1‑149 (MQ=255)                                                ttATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑t                                                                                                 >  2:104452/1‑149 (MQ=255)                                                  aTCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCGCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTCCGGTCTGCCACCAGGTATCGACCACCGAACATTTCTCGTTGCCGA‑ttt                                                                                               >  4:690050/1‑149 (MQ=255)                                                    cGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑ttttt                                                                                             <  2:306673/149‑1 (MQ=255)                                                      aCCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATAGACCACCGGACATTTCTCGTTGCCGA‑ttttttt                                                                                           <  3:665119/149‑1 (MQ=255)                                                       ccAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTc                                                                                          <  4:489585/149‑1 (MQ=255)                                                       ccAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTc                                                                                          >  2:175703/1‑149 (MQ=255)                                                        cAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTcc                                                                                         <  2:890792/149‑1 (MQ=255)                                                          gcgcCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAg                                                                                       <  3:185936/149‑1 (MQ=255)                                                           cgcCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGtt                                                                                                        >  2:420962/1‑132 (MQ=255)                                                           cgcCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGtt                                                                                                        <  1:420962/132‑1 (MQ=255)                                                           cgcCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGt                                                                                      <  3:448928/149‑1 (MQ=255)                                                                   gCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTcc                                                                              >  4:24103/1‑149 (MQ=255)                                                                     aCGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑ttttt                                                                                             >  2:533549/1‑132 (MQ=255)                                                                     aCGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑ttttt                                                                                             <  1:533549/132‑1 (MQ=255)                                                                     aCGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCa                                                                            <  4:400644/149‑1 (MQ=255)                                                                      cGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCAc                                                                           <  4:720502/149‑1 (MQ=255)                                                                       gCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACg                                                                          >  4:196787/1‑149 (MQ=255)                                                                       gCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACg                                                                          >  2:392925/1‑149 (MQ=255)                                                                           aagaaCGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTc                                                                      >  3:917913/1‑149 (MQ=255)                                                                             gaaCGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCg                                                                    >  1:1042877/1‑149 (MQ=255)                                                                              aaCGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCAGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCgg                                                                   <  3:976145/149‑1 (MQ=255)                                                                                cGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGt                                                                 <  2:741701/149‑1 (MQ=255)                                                                                cGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGt                                                                 <  2:515237/149‑1 (MQ=255)                                                                                  gACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑ttttt                                                                                             <  2:887823/119‑1 (MQ=255)                                                                                  gACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑ttttt                                                                                             >  1:887823/1‑119 (MQ=255)                                                                                     gtgtTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAAtt                                                            <  3:24103/149‑1 (MQ=255)                                                                                       gtTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCATTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTgg                                                          <  4:852601/149‑1 (MQ=255)                                                                                             gAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGcc                                                    >  3:390215/1‑149 (MQ=255)                                                                                                  ggCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACgg                                               >  3:255427/1‑149 (MQ=255)                                                                                                    cTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGaa                                             <  3:934320/149‑1 (MQ=255)                                                                                                    cTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGaa                                             >  4:536946/1‑149 (MQ=255)                                                                                                     tttCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAAc                                            >  3:596632/1‑149 (MQ=255)                                                                                                     tttCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAAc                                            >  3:89237/1‑149 (MQ=255)                                                                                                     tttCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAAc                                            >  2:291023/1‑149 (MQ=255)                                                                                                        cAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCga                                         <  1:852545/149‑1 (MQ=255)                                                                                                          gCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCgaga                                       >  1:568417/1‑149 (MQ=255)                                                                                                          gCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCgaga                                       >  4:472187/1‑149 (MQ=255)                                                                                                          gCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCgaga                                       <  1:164202/149‑1 (MQ=255)                                                                                                           cTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATc                                                                                                                             <  1:876075/63‑1 (MQ=255)                                                                                                           cTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATc                                                                                                                             >  2:876075/1‑63 (MQ=255)                                                                                                                      cAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGcgac                           >  2:666079/1‑149 (MQ=255)                                                                                                                      cAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGcgac                           <  1:812153/149‑1 (MQ=255)                                                                                                                      cAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGcgac                           >  4:548013/1‑149 (MQ=255)                                                                                                                        ggCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGcgacga                         <  4:900310/149‑1 (MQ=255)                                                                                                                         gCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGaa                        >  1:1008686/1‑149 (MQ=255)                                                                                                                         gCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGaa                        >  2:550540/1‑149 (MQ=255)                                                                                                                           aGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAAcg                      <  2:818920/149‑1 (MQ=255)                                                                                                                               gggTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAAcggtcg                  <  4:713919/149‑1 (MQ=255)                                                                                                                                  tGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATAGACCACCTGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCTACGCTTCCGCGTTAATTTGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTg               >  4:429458/1‑149 (MQ=255)                                                                                                                                   gATCATGAGACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGc              <  4:169373/149‑1 (MQ=255)                                                                                                                                   gATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTg               <  3:87752/149‑1 (MQ=255)                                                                                                                                   gATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGc              <  1:486256/149‑1 (MQ=255)                                                                                                                                    aTCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGcc             <  1:392925/149‑1 (MQ=255)                                                                                                                                    aTCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGcc             >  3:317386/1‑149 (MQ=255)                                                                                                                                     gaccgAAACCGGCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCt            >  2:467658/5‑149 (MQ=255)                                                                                                                                       aTGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTc          >  4:887462/1‑149 (MQ=255)                                                                                                                                       aTGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTc          >  3:612379/1‑149 (MQ=255)                                                                                                                                       aTGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTc          >  3:222737/1‑149 (MQ=255)                                                                                                                                       aTGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTc          >  4:752075/1‑149 (MQ=255)                                                                                                                                       aTGAAACCGCCAGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTc          >  1:557551/1‑149 (MQ=255)                                                                                                                                        tGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCg         >  2:131360/1‑149 (MQ=255)                                                                                                                                        cgAAACCGGCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCg         >  1:919760/2‑149 (MQ=255)                                                                                                                                         gAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGa        <  1:964803/149‑1 (MQ=255)                                                                                                                                          aaaCCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGAt       >  2:94853/1‑149 (MQ=255)                                                                                                                                             ccgccgGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGc    <  3:752075/149‑1 (MQ=255)                                                                                                                                              cgccgGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCt   >  4:665178/1‑149 (MQ=255)                                                                                                                                              cgccgGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCt   <  1:970591/149‑1 (MQ=255)                                                                                                                                              cgccgGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCt   <  2:246181/149‑1 (MQ=255)                                                                                                                                               gccgGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCtt  >  2:224108/1‑149 (MQ=255)                                                                                                                                               gccgGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCtt  >  4:119493/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     TGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGA‑TTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCTT  >  NC_000913/4285973‑4286263

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑