Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,286,116	C→G	8.1%	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	8.1%	202.6 / 9.5	62	G419A (GGC→GCC)	acs	acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (34/23);  new base G (2/3);  total (36/26)
Fisher's exact test for biased strand distribution p-value = 6.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

TCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCT  >  NC_000913/4285969‑4286262                                                                                                                                                    |                                                                                                                                                   tCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAAccgc                                                                                                                                                    <  1:963914/148‑1 (MQ=255) tCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAAccgc                                                                                                                                                    <  1:511233/148‑1 (MQ=255)   ggTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACTGCCgg                                                                                                                                                 >  3:571526/1‑149 (MQ=255)    gTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCgg                                                                                                                                                 >  1:335182/1‑148 (MQ=255)      gATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGtt                                                                                                                                               <  2:475552/148‑1 (MQ=255)        ttACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCg                                                                                                                                            <  4:603346/149‑1 (MQ=255)         tACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGGCGGTTTca                                                                                                                                            <  1:336038/148‑2 (MQ=255)         tACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGGCGGTTTca                                                                                                                                            <  1:891668/148‑2 (MQ=255)            cAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCt                                                                                                                                        <  3:264899/149‑1 (MQ=255)             aGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCt                                                                                                                                        >  2:513844/1‑148 (MQ=255)               gCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCACGAAACCGCCGGTTTCGGTCTGc                                                                                                                                      >  1:906318/1‑148 (MQ=255)                    cTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAg                                                                                                                                <  4:571526/149‑1 (MQ=255)                              ccccTCCAGCAGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGacca                                                                                                                      >  3:61338/1‑149 (MQ=255)                               cccTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGacca                                                                                                                      >  2:313436/1‑148 (MQ=255)                                   ccAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCgg                                                                                                                  >  2:650528/1‑148 (MQ=255)                                       cGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACAtt                                                                                                             <  3:911529/149‑1 (MQ=255)                                         ggTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACAttt                                                                                                            >  2:505171/1‑148 (MQ=255)                                                ttCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTg                                                                                                     <  2:910828/148‑1 (MQ=255)                                                ttCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTg                                                                                                     >  1:898311/1‑148 (MQ=255)                                                ttCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGc                                                                                                    <  4:629793/149‑1 (MQ=255)                                                ttCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGc                                                                                                    >  4:281291/1‑149 (MQ=255)                                                 tCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGGCGGTTTca                                                                                                                                            <  1:785766/108‑2 (MQ=255)                                                 tCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGGCGGTTTca                                                                                                                                            >  2:785766/1‑107 (MQ=255)                                                 tCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTc                                                                                                                                         <  1:483879/111‑1 (MQ=255)                                                 tCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTc                                                                                                                                         >  2:483879/1‑111 (MQ=255)                                                      aTCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGAttt                                                                                              >  3:473019/1‑149 (MQ=255)                                                           ccAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTc                                                                                         <  4:701125/149‑1 (MQ=255)                                                             aGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTcc                                                                                        >  2:382966/1‑148 (MQ=255)                                                                   ggTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTAcc                                                                                  <  2:87771/148‑1 (MQ=255)                                                                    gTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCa                                                                                 >  1:84404/1‑148 (MQ=255)                                                                         aCGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCa                                                                           >  4:575993/1‑149 (MQ=255)                                                                              gaagaaCGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCtt                                                                      >  4:480701/1‑149 (MQ=255)                                                                                agaaCGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTc                                                                     >  2:336111/1‑148 (MQ=255)                                                                                 gaaCGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTcc                                                                    <  1:688195/148‑1 (MQ=255)                                                                                  aaCGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCg                                                                   <  1:313436/148‑1 (MQ=255)                                                                                   aCGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCgg                                                                  <  2:213113/148‑1 (MQ=255)                                                                                     ggACATGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGt                                                                >  2:398577/1‑148 (MQ=255)                                                                                       aCGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTa                                                              <  2:335182/148‑1 (MQ=255)                                                                                          tgtTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTCTCGGTCTCCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAAtt                                                           <  2:690971/148‑1 (MQ=255)                                                                                             tGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCt                                                       <  4:399333/149‑1 (MQ=255)                                                                                               ccGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTc                                                      >  2:64263/1‑148 (MQ=255)                                                                                                 gAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGc                                                    <  1:579841/148‑1 (MQ=255)                                                                                                  aaCCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGcc                                                   >  1:943720/1‑148 (MQ=255)                                                                                                    ccGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCAc                                                <  3:480701/149‑1 (MQ=255)                                                                                                        cTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGa                                             >  1:480141/1‑148 (MQ=255)                                                                                                         tttCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCCGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAAc                                           >  3:407641/1‑149 (MQ=255)                                                                                                              gCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCgaga                                      >  3:628100/1‑149 (MQ=255)                                                                                                              gCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCgag                                       >  2:250024/1‑148 (MQ=255)                                                                                                               cTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCgaga                                      >  1:140102/1‑148 (MQ=255)                                                                                                               cTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCgaga                                      >  1:126660/1‑148 (MQ=255)                                                                                                                tCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAAt                                    <  4:832798/149‑1 (MQ=255)                                                                                                                  ggTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATg                                   <  1:513844/148‑1 (MQ=255)                                                                                                                   gTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATgc                                  >  1:412743/1‑148 (MQ=255)                                                                                                                           aGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGcgacg                         >  4:186731/1‑149 (MQ=255)                                                                                                                           aGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGcgacg                         >  4:81252/1‑149 (MQ=255)                                                                                                                                 cGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAAcgg                    >  2:459148/1‑148 (MQ=255)                                                                                                                                   gggTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGGGAATGCGCAGCGACGAAcggtcg                 <  4:197366/149‑1 (MQ=255)                                                                                                                                      tGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAAcggtcggt               >  2:305708/1‑148 (MQ=255)                                                                                                                                        aTCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGcc            >  4:293359/1‑149 (MQ=255)                                                                                                                                           aTGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTc         <  3:575993/149‑1 (MQ=255)                                                                                                                                           aTGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTc         >  4:430785/1‑149 (MQ=255)                                                                                                                                             gAAACCGGCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCg        >  1:768093/1‑148 (MQ=255)                                                                                                                                             gAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCg        <  1:832956/148‑1 (MQ=255)                                                                                                                                               aaCCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATc     >  4:97014/1‑149 (MQ=255)                                                                                                                                                  cgccgGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGc   <  1:725972/148‑1 (MQ=255)                                                                                                                                                  cgccgGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGc   <  2:84404/148‑1 (MQ=255)                                                                                                                                                  cgccgGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCt  <  4:427303/149‑1 (MQ=255)                                                                                                                                                  cgccgGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCt  <  3:814839/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TCGGTGATTACCAGGCTACCTTCGGTGGCCCCCTCCAGCGGGTTACCTTCGTTATCGACCAGCGCCGGTTGCACGCCGAAGAACGGACGTGTTGCCGAACCGGCTTTCAGCTCGGTAGCGCCAGGCAGCGGGGTGATCATGAAACCGCCGGTTTCGGTCTGCCACCAGGTATCGACCACCGGACATTTCTCGTTGCCGATTTTTTTCCAGTACCACTCCCACGCTTCCGGGTTAATTGGCTCGCCCACGGAACCGAGAATGCGCAGCGACGAACGGTCGGTGCCTTCGATCGCT  >  NC_000913/4285969‑4286262

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