Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,155,575	G→A	7.0%	R151H (CGT→CAT)	tmk →	thymidylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,155,575	0	G	A	7.0%	276.0 / 7.4	86	R151H (CGT→CAT)	tmk	thymidylate kinase
Reads supporting (aligned to +/- strand):  ref base G (44/36);  new base A (2/4);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 4.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

TCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTGATGGATG  >  NC_000913/1155429‑1155718                                                                                                                                                   |                                                                                                                                                tctctcCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTg                                                                                                                                               <  2:483722/149‑1 (MQ=255)    ctcCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTgcgc                                                                                                                                            >  3:752635/1‑149 (MQ=255)    ctcCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTgcgc                                                                                                                                            >  2:580028/1‑149 (MQ=255)     tcCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTgcgcg                                                                                                                                           >  3:315034/1‑149 (MQ=255)        aCTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGTAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCAcgcgcgcgt                                                                                                                                        <  2:96997/149‑2 (MQ=255)        aCTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCgcgg                                                                                                                                        <  3:55650/149‑1 (MQ=255)        aCTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCAcgcgcgcgt                                                                                                                                        <  4:592780/149‑2 (MQ=255)                 tATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGa                                                                                                                               <  4:315034/149‑1 (MQ=255)                  aTCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGAt                                                                                                                              >  4:508738/1‑149 (MQ=255)                     aGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGt                                                                                                                           >  2:369823/1‑149 (MQ=255)                      gggCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTATCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTa                                                                                                                          <  3:756146/149‑1 (MQ=255)                      gggCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTa                                                                                                                          <  1:583399/149‑1 (MQ=255)                       ggcggcGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTAt                                                                                                                         >  3:319522/1‑149 (MQ=255)                        gcggcgGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTAtt                                                                                                                        >  4:181464/1‑149 (MQ=255)                        gcggcgGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTAtt                                                                                                                        <  2:145623/149‑1 (MQ=255)                          ggcggACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGa                                                                                                                      <  3:557696/149‑1 (MQ=255)                           gcggACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAg                                                                                                                     <  4:488457/149‑1 (MQ=255)                              gACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCaa                                                                                                                  <  2:51508/149‑1 (MQ=255)                               aCGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAg                                                                                                                 >  3:534041/1‑149 (MQ=255)                                  tGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAAt                                                                                                              <  3:181464/149‑1 (MQ=255)                                   ggTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCAcgcgcgcgtt                                                                                                                                       <  1:321872/123‑3 (MQ=255)                                   ggTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCAcgcgcgcgtt                                                                                                                                       >  2:321872/1‑121 (MQ=255)                                    gTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCt                                                                                                            >  4:578117/1‑149 (MQ=255)                                     tATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCtt                                                                                                           >  2:570585/1‑149 (MQ=255)                                     tATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGAAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCtt                                                                                                           >  1:551667/1‑149 (MQ=255)                                      aTTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCttt                                                                                                          <  3:435420/149‑1 (MQ=255)                                      aTTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCttt                                                                                                          <  2:389722/149‑1 (MQ=255)                                       ttGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTc                                                                                                         <  1:456101/149‑1 (MQ=255)                                       ttGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTc                                                                                                         <  3:98484/149‑1 (MQ=255)                                         gACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGa                                                                                                       <  3:119488/149‑1 (MQ=255)                                         gACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGa                                                                                                       <  1:249917/149‑1 (MQ=255)                                                atatGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCttt                                                                                                <  1:291873/149‑1 (MQ=255)                                                   tGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCttt                                                                                                          >  2:754956/1‑136 (MQ=255)                                                   tGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCttt                                                                                                          <  1:754956/136‑1 (MQ=255)                                                    gCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATc                                                                                            >  3:683419/1‑149 (MQ=255)                                                     cTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCg                                                                                           >  3:510140/1‑149 (MQ=255)                                                        gCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCAc                                                                                        >  2:318659/1‑149 (MQ=255)                                                        gCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCAc                                                                                        <  4:660514/149‑1 (MQ=255)                                                             acTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGc                                                                                                                    <  3:734146/116‑1 (MQ=255)                                                             acTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGc                                                                                                                    >  4:734146/1‑116 (MQ=255)                                                             acTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCcgcg                                                                                   >  4:614837/1‑149 (MQ=255)                                                                     aTGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTAt                                                                           <  2:521170/149‑1 (MQ=255)                                                                      tGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATc                                                                          >  2:634418/1‑149 (MQ=255)                                                                           ttCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGaa                                                                     <  2:590330/149‑1 (MQ=255)                                                                           ttCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCAcgcgcgcgtt                                                                                                                                       >  2:704508/1‑81 (MQ=255)                                                                           ttCTCGGGGATTTTCGCCCCGACTTAACGATCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCAcgcgcgcgtt                                                                                                                                       <  1:704508/83‑3 (MQ=255)                                                                              tcGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTg                                                                  >  4:586093/1‑149 (MQ=255)                                                                                ggggATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGgc                                                                >  4:158631/1‑149 (MQ=255)                                                                                ggggATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGgc                                                                >  4:349223/1‑149 (MQ=255)                                                                                  ggATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGgcag                                                              >  2:771536/1‑149 (MQ=255)                                                                                     ttttCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGcac                                                           >  4:508665/1‑149 (MQ=255)                                                                                       ttCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACaa                                                         <  4:386375/149‑1 (MQ=255)                                                                                           ccccGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATa                                                     >  3:195615/1‑149 (MQ=255)                                                                                            cccGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGTTTAAAACGCGCGCGTGCGCACGGCGAGCTGGAACGTATTGAGCAAGAATCTTTAGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATaa                                                    >  4:281711/1‑149 (MQ=255)                                                                                            cccGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATaa                                                    >  3:110735/1‑149 (MQ=255)                                                                                            cccGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATaa                                                    >  3:110488/1‑149 (MQ=255)                                                                                              cGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATaaaa                                                  <  1:749720/149‑1 (MQ=255)                                                                                                  ttAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCAt                                              >  2:294140/1‑149 (MQ=255)                                                                                                      cGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCAt                                          >  2:333772/1‑149 (MQ=255)                                                                                                        ctctATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTc                                                                                                   <  1:163508/90‑1 (MQ=255)                                                                                                        ctctATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTc                                                                                                   >  2:163508/1‑90 (MQ=255)                                                                                                        ctctATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATAc                                        <  1:466299/149‑1 (MQ=255)                                                                                                               tcGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGTATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGAt                                 >  3:585264/1‑149 (MQ=255)                                                                                                               tcGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTgg                                                                 <  2:761705/117‑1 (MQ=255)                                                                                                               tcGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTgg                                                                 >  1:761705/1‑117 (MQ=255)                                                                                                               tcGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCa                                               <  1:512400/135‑1 (MQ=255)                                                                                                               tcGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCa                                               >  2:512400/1‑135 (MQ=255)                                                                                                               tcGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGAt                                 >  1:144526/1‑149 (MQ=255)                                                                                                               tcGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGAt                                 >  3:621142/1‑149 (MQ=255)                                                                                                               tcGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGAt                                 <  4:510140/149‑1 (MQ=255)                                                                                                                 gATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGc                               >  2:324309/1‑149 (MQ=255)                                                                                                                      tACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCa                                               >  4:268739/1‑128 (MQ=255)                                                                                                                      tACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCa                                               >  2:176510/1‑128 (MQ=255)                                                                                                                      tACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCa                                               <  1:176510/128‑1 (MQ=255)                                                                                                                      tACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCa                                               <  3:268739/128‑1 (MQ=255)                                                                                                                          ccGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGGTAAAAGCATTCATACCATTGATGCCACCCAGcc                      >  1:431194/1‑149 (MQ=255)                                                                                                                          ccGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGcc                      <  2:422696/149‑1 (MQ=255)                                                                                                                          ccGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGcc                      >  1:672653/1‑149 (MQ=255)                                                                                                                            ggAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGc                    <  3:586093/149‑1 (MQ=255)                                                                                                                                     cTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCg           <  4:275137/149‑1 (MQ=255)                                                                                                                                      ttAAAACTCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCATAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGt          >  2:417019/1‑149 (MQ=255)                                                                                                                                      ttAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGt          >  2:141611/1‑149 (MQ=255)                                                                                                                                        aaaaCGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCGGCTGGAGGCCGTga        <  1:580028/149‑1 (MQ=255)                                                                                                                                        aaaaCGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCAGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTga        >  1:429381/1‑149 (MQ=255)                                                                                                                                          aaCGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTgatg      >  4:259017/1‑149 (MQ=255)                                                                                                                                           aCGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTgatgg     <  3:578117/149‑1 (MQ=255)                                                                                                                                             gcgcgcgTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTgatggat   <  3:635219/149‑1 (MQ=255)                                                                                                                                             gcgcgcgTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTgatggat   <  3:164456/149‑1 (MQ=255)                                                                                                                                             gcgcgcgTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTgatggat   <  4:752635/149‑1 (MQ=255)                                                                                                                                              cgcgcgTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTgatggatg  <  2:439627/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                TCTCTCCACTCAGGCGTATCAGGGCGGCGGACGTGGTATTGACCAACATATGCTGGCAACACTGCGTGATGCTGTTCTCGGGGATTTTCGCCCCGACTTAACGCTCTATCTCGATGTTACCCCGGAAGTTGGCTTAAAACGCGCGCGTGCGCGCGGCGAGCTGGATCGTATTGAGCAAGAATCTTTCGATTTCTTTAATCGCACCCGCGCCCGCTATCTGGAACTGGCAGCACAAGATAAAAGCATTCATACCATTGATGCCACCCAGCCGCTGGAGGCCGTGATGGATG  >  NC_000913/1155429‑1155718

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