Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,927,702	A→T	5.5%	S380T (TCT→ACT)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,927,702	0	A	T	5.5%	240.6 / 3.0	73	S380T (TCT→ACT)	ptrB	protease II
Reads supporting (aligned to +/- strand):  ref base A (32/37);  new base T (2/2);  total (34/39)
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.61e-01

TTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTC  >  NC_000913/1927556‑1927849                                                                                                                                                   |                                                                                                                                                    ttCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGgaagaa                                                                                                                                                    >  7:414922/1‑148 (MQ=255)    gACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATaa                                                                                                                                                 <  6:224310/148‑1 (MQ=255)    gACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAAc                                                                                                                                                <  7:175045/149‑1 (MQ=255)     aCGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAAc                                                                                                                                                >  7:244462/1‑148 (MQ=255)         cATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATaa                                                                                                                                           <  1:220657/149‑1 (MQ=255)         cATCACGCGCGACTATCCACATGTGTTCAATGCGGTAATTCGCCGCATAAAAACCAGTAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATaa                                                                                                                                           <  2:193329/149‑1 (MQ=255)          aTCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAAc                                                                                                                                          >  3:71122/1‑149 (MQ=255)          aTCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAAc                                                                                                                                          <  6:176512/149‑1 (MQ=255)           tCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACg                                                                                                                                         >  3:319726/1‑149 (MQ=255)               gggCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCa                                                                                                                                              >  1:15821/1‑140 (MQ=255)               gggCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCa                                                                                                                                              <  2:15821/140‑1 (MQ=255)                   gACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATcgcg                                                                                                                                 >  1:110120/1‑149 (MQ=255)                      tATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGTATAACCATAACGCAATCGCGCgg                                                                                                                              >  2:165413/1‑149 (MQ=255)                      tATCCACAGGTGATCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGTATAACCATAACGCAATCGCGCgg                                                                                                                              >  2:159388/1‑149 (MQ=255)                        tCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGtt                                                                                                                            >  1:386556/1‑149 (MQ=255)                          cacaGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATcgcgc                                                                                                                                >  3:322158/1‑143 (MQ=255)                             aGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGGTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAgg                                                                                                                       >  3:224584/1‑149 (MQ=255)                                 gtTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCt                                                                                                                   >  3:13064/1‑149 (MQ=255)                                     aCTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCTGTATCCATATCCAGTTCAAACAAACTGTCTTGTGTAGTCATGCAAGAATAACCATAACGCAATCGAGCGGTATCACGTTCTGGAt                                                                                                               >  3:82302/1‑149 (MQ=255)                                     aCTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGAt                                                                                                               >  8:163222/1‑149 (MQ=255)                                     aCTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGGCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGa                                                                                                                <  1:333179/148‑1 (MQ=255)                                        gCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAAt                                                                                                                                                   >  1:79089/1‑110 (MQ=255)                                        gCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAAt                                                                                                                                                   <  2:79089/110‑1 (MQ=255)                                         cGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGt                                                                                                            <  2:120118/148‑1 (MQ=255)                                         cGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTa                                                                                                           <  1:404270/149‑1 (MQ=255)                                         cGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTa                                                                                                           <  7:294982/149‑1 (MQ=255)                                         cGGTAAATCGCCGCATAAAAACCATGAACTTCCGTTTGTTTTTATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTa                                                                                                           <  5:219237/149‑1 (MQ=255)                                          ggTAATTCGCCTCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAg                                                                                                          <  2:280373/149‑1 (MQ=255)                                            tAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGTAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAgg                                                                                                         >  6:223150/1‑148 (MQ=255)                                              aTTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCaa                                                                                                      <  4:319726/149‑1 (MQ=255)                                              aTTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAAACATAACGCAATCGCGCGGTTTCAGGTTCTGCAATGTAGGcca                                                                                                      >  4:28975/1‑147 (MQ=255)                                               ttCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAAt                                                                                                     <  3:241814/149‑1 (MQ=255)                                                        tAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCAAGGTca                                                                                             >  5:409243/1‑148 (MQ=255)                                                           aaaCCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATa                                                                                         <  8:424038/149‑1 (MQ=255)                                                            aaCCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAg                                                                                        <  2:110120/149‑1 (MQ=255)                                                                  gAACTTCCGTTTGTTTTAATACTAGATGCTAACCGGTATCCATATCTAGTTCAAACATAGTGTCTGTTGTATTCATTGAATAATAACCATTAGGCTATCTCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACTTATGCCGGa                                                                                  <  4:71122/149‑1 (MQ=255)                                                                          gtttgtttTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCaaa                                                                          <  7:257349/149‑1 (MQ=255)                                                                          gtttgtttTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCaaa                                                                          <  2:386556/149‑1 (MQ=255)                                                                           tttgtttTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAg                                                                         <  8:414922/149‑1 (MQ=255)                                                                              gtttTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGc                                                                       <  7:163222/148‑1 (MQ=255)                                                                                tttAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTAATGGAAGAATAACCATAAAGCAATCGCGCGGTTTAAGGTTCTGGATTTTAGGGAATACAGGTCACATAGTCCGGATCATCAAAGGCaa                                                                     >  5:289487/1‑148 (MQ=255)                                                                                  tAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATa                                                                   <  4:307847/148‑1 (MQ=255)                                                                                    aTACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAgg                                                                                                         <  4:131437/108‑1 (MQ=255)                                                                                    aTACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAgg                                                                                                         >  3:131437/1‑108 (MQ=255)                                                                                     tACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAACAATAACCATAAAGCAATCGCGCGGTTTCAGGTTCTGGAGTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCaa                                                               <  5:223150/149‑1 (MQ=255)                                                                                       cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGTATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAAt                                                              <  1:159388/148‑1 (MQ=255)                                                                                       cacGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAAt                                                              <  1:51608/148‑1 (MQ=255)                                                                                       cacGACGCTCACCGGTATCCATATCCAGTTAAAACAAAGTGTCTGGTGTAGTCATGGAAGTATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAAt                                                              <  1:165413/148‑1 (MQ=255)                                                                                        acgacgCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAAAGa                                                            >  8:350856/1‑149 (MQ=255)                                                                                         cgacgCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGAc                                                           <  2:56020/149‑1 (MQ=255)                                                                                          gacgCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATAc                                                                  <  1:291862/141‑1 (MQ=255)                                                                                          gacgCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATAc                                                                  >  2:291862/1‑141 (MQ=255)                                                                                             gCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTcc                                                       <  6:259176/149‑1 (MQ=255)                                                                                               tCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCg                                                     >  5:62473/1‑149 (MQ=255)                                                                                                          catatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAAt                                                              >  8:325890/1‑129 (MQ=255)                                                                                                           atatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGcc                                                                                     <  1:420846/105‑1 (MQ=255)                                                                                                           atatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGcc                                                                                     >  2:420846/1‑105 (MQ=255)                                                                                                           atatccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGt                                          >  8:76086/1‑148 (MQ=255)                                                                                                              tccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTaa                                       <  7:76086/148‑1 (MQ=255)                                                                                                              tccaGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAAt                                      >  2:39599/1‑149 (MQ=255)                                                                                                                  gTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTgg                                  <  5:68256/149‑1 (MQ=255)                                                                                                                  gTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTgg                                  <  5:409487/149‑1 (MQ=255)                                                                                                                  gTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGACTTCCGGTTAATTTgg                                  >  2:88055/1‑149 (MQ=255)                                                                                                                   ttCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGgc                                 <  4:45692/149‑1 (MQ=255)                                                                                                                   ttCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGgc                                 <  4:224584/149‑1 (MQ=255)                                                                                                                    tCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGgcg                                >  4:46667/1‑149 (MQ=255)                                                                                                                           aaGTGTCTGGTGTAGTCATTGACGAATAACCATAACTCAATCGCGCGGTTTCAGGGTCTTGATTGTCGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTTGCGCAAACTg                         <  6:289487/149‑1 (MQ=255)                                                                                                                             gtgtCTGGTGTAGTCATGGAAGAATAACCATAACGAAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACAAATGACTTCCAGGGTCTTGCGGTTAATTTGGCTAAAACTGGt                       >  1:386219/1‑149 (MQ=255)                                                                                                                                   ggTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAAccc                 >  2:81755/1‑149 (MQ=255)                                                                                                                                   ggTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAAccc                 >  4:377882/1‑149 (MQ=255)                                                                                                                                     tgtAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCcgc               <  3:377882/149‑1 (MQ=255)                                                                                                                                          tCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGa          >  1:193886/1‑149 (MQ=255)                                                                                                                                            aTGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACg        >  4:251950/1‑149 (MQ=255)                                                                                                                                             tGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGc       >  8:273979/1‑149 (MQ=255)                                                                                                                                                 agaaTAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCtt   >  5:130377/1‑149 (MQ=255)                                                                                                                                                  gaaTAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATAGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTc  <  7:169931/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    TTCGACGCCATCACGGGCGACTATCCACAGGTGTTCACTGCGGTAATTCGCCGCATAAAAACCAGGAACTTCCGTTTGTTTTAATACACGACGCTCACCGGTATCCATATCCAGTTCAAACAAAGTGTCTGGTGTAGTCATGGAAGAATAACCATAACGCAATCGCGCGGTTTCAGGTTCTGGATTGTAGGCAATCCAGGTCACATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTC  >  NC_000913/1927556‑1927849

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