Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	733,768	G→A	5.7%	V1396M (GTG→ATG)	rhsC →	rhs element protein RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	733,768	0	G	A	5.7%	191.3 / 2.2	70	V1396M (GTG→ATG)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base G (36/30);  new base A (2/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

TTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTATTAGACTCCACTCGTGATGTGAGCGAGTTAATTAGTCTCA  >  NC_000913/733630‑733902                                                                                                                                           |                                                                                                                                       ttGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGAt                                                                                                                                    >  2:458303/1‑143 (MQ=255)    ccTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCGAAGGGTGGTAGAACGATTCTTTCGCGGAAGTAAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTa                                                                                                                                 >  1:218059/1‑143 (MQ=255)    ccTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAATGATGTa                                                                                                                                 >  2:272768/1‑143 (MQ=255)    ccTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAATGATGTa                                                                                                                                 >  2:888778/1‑143 (MQ=255)     cTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAGCATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTaa                                                                                                                                <  1:698563/143‑1 (MQ=255)            atTCACTCCGATTAATTTAGGGCTATTTAAAATCTGTCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGACGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAata                                                                                                                         <  2:259415/143‑1 (MQ=255)             ttCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAatat                                                                                                                        >  2:749332/1‑143 (MQ=255)              tCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATg                                                                                                                                   >  1:553771/1‑131 (MQ=255)              tCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATg                                                                                                                                   <  2:553771/131‑1 (MQ=255)               cacaCCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTTTTTAGAAACATGTAAGATGAAGTGCTCCGGAACTAACATTGGACGTTTTATCAGCTTTGTATTCACCGGAGTGATGTAATGAAatatat                                                                                                                      >  1:31179/1‑143 (MQ=255)                acacCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGACATAtatt                                                                                                                     >  2:475689/1‑143 (MQ=255)                     tataAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGCCGTTTTATTAGATTTGTATTCACCGGAGTGATGTACTGAAATATATTAtttt                                                                                                                >  2:429154/1‑143 (MQ=255)                     tataAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTAtttt                                                                                                                <  1:601989/143‑1 (MQ=255)                     tataAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTAtttt                                                                                                                >  1:594647/1‑143 (MQ=255)                       ttaaTTTAGGGCTATTTAAAAACTGCCGAGAGTGGTCTAATGATTCTTTAGAAGCTTGTAAGAGGAAGTGCTGCGGGACGAACATCGGACGTTTTATTAGATTTGTAGTCACCGGAGTGATGTAATGAAATATATTAtttttt                                                                                                              <  2:871987/141‑1 (MQ=255)                           tttAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTAtt                                                                                                          <  1:126825/143‑1 (MQ=255)                               gggCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCaga                                                                                                      <  2:160104/143‑1 (MQ=255)                                ggCTATTTAAGAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCagag                                                                                                     <  1:867370/143‑1 (MQ=255)                                ggCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCagag                                                                                                     >  1:79149/1‑143 (MQ=255)                                  cTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCt                                                                                                   >  1:198160/1‑143 (MQ=255)                                   tATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTACTCCCCGGAGCGATGTAATGAAATATATTATTTTTTTATTCAGAGCta                                                                                                  >  2:413295/1‑143 (MQ=255)                                      tttaaagCTGCCCAGACGGGTCTAATGATTCTGGAGAAACAGGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTatat                                                                                               <  2:374192/140‑1 (MQ=255)                                       tAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATAtg                                                                                              >  1:598111/1‑143 (MQ=255)                                        aaaaaCTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATAtgg                                                                                             <  1:180245/143‑1 (MQ=255)                                           aaCTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCGATAtggctg                                                                                          >  2:123714/1‑143 (MQ=255)                                            aCTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATAtggctgg                                                                                         <  2:716947/143‑1 (MQ=255)                                               gacAAGACTGGTTTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCtt                                                                                      <  1:483662/141‑1 (MQ=255)                                                 cAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGCTTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTa                                                                                    >  2:96616/1‑143 (MQ=255)                                                  aaGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTAt                                                                                   >  1:633932/1‑143 (MQ=255)                                                   taaaacgTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGTTTTATc                                                                                  >  2:621293/7‑143 (MQ=255)                                                     aCTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGATCTATATGTCTGGCTTTATCgt                                                                                >  2:170534/1‑142 (MQ=255)                                                        ggTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCttt                                                                             >  1:29214/1‑143 (MQ=255)                                                         gTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCtttt                                                                            <  1:98496/143‑1 (MQ=255)                                                         gTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCt                                                                               <  1:692504/140‑1 (MQ=255)                                                         gTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCt                                                                               >  2:692504/1‑140 (MQ=255)                                                          gctgaTGATTCTTGAGAGACATGTAAGGTGAAGTGCTCCGGAAAGAACATTGGACGTTTTATTGGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTa                                                                           <  2:373238/139‑1 (MQ=255)                                                           cTAATGATTCTTTAGAAACATGGAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTaa                                                                          <  1:799714/143‑1 (MQ=255)                                                           cTAATGATTCTGTAGAAACATGTAAGAGGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTaa                                                                          <  1:380639/143‑1 (MQ=255)                                                            tAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGCGCTATATGGCTGGCTGTATCCCTTTTCAt                                                                         >  1:147028/1‑143 (MQ=255)                                                             aaTGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAATGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATa                                                                        <  1:272768/143‑1 (MQ=255)                                                             aaTGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAATGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATa                                                                        <  1:888778/143‑1 (MQ=255)                                                                gATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCTGAGCTATATGGCTGTCTTTATCGTTTTTATTGCTa                                                                     >  2:823718/1‑143 (MQ=255)                                                                gATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGAGGTGATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTTTCGTTTTTAATACTa                                                                     >  2:564693/1‑143 (MQ=255)                                                                    cTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCt                                                                 >  1:274639/1‑143 (MQ=255)                                                                        aGAAACATGTAAGATGAAGTGCTCCGGCAATAACATTGGACGTTTTATTATATTTGTATTCACCGCAGTGATTTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTc                                                             <  1:95713/143‑1 (MQ=255)                                                                            agaTGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATg                                                         <  1:519366/141‑1 (MQ=255)                                                                                   aGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGAc                                                  <  1:258027/143‑1 (MQ=255)                                                                                    gATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCGTGCATAGACt                                                 >  1:717417/1‑143 (MQ=255)                                                                                        aaGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCa                                             <  2:308373/143‑1 (MQ=255)                                                                                           tGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTAACAGGAGAGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTa                                          >  1:849190/1‑143 (MQ=255)                                                                                            tgtCCAGAAATAACTTTTGAGTTTGTATTGTTTTTTTATTCACCGGAGCTATCTAATTAAATATATTATTTTTTTATTCATATCTCTATTGCTGTCTTTATCGCTTTTAATACTATTCTTTTCCATGCATCGACTTTCACTAt                                         <  2:688332/141‑1 (MQ=255)                                                                                             cTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAGTGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTAtt                                        >  1:290214/1‑143 (MQ=255)                                                                                             cTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTAtt                                        <  1:93353/143‑1 (MQ=255)                                                                                             cTCCGCAAATAACATTGGACGTTTTATTAGATTTGTATTCACTGGAGTGATGTAATGAAATATTTTGCTGTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTAtt                                        <  2:451087/143‑1 (MQ=255)                                                                                              tCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGGGTGATGTAACGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGCCTTTCACTATTa                                       >  1:489443/1‑143 (MQ=255)                                                                                               ccGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGTTATATGGCTGGCTTTATCGCTTTTAATACTCTTCTTTTTCATGCATAGAATTTGattatttt                                      >  2:236205/1‑141 (MQ=255)                                                                                                  gAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTATTAGACt                                   >  1:450925/1‑143 (MQ=255)                                                                                                       aaCATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTATTAGACTCCACt                              >  1:99599/1‑143 (MQ=255)                                                                                                          aTTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTATTAGACTCCACTCGt                           <  1:747501/143‑1 (MQ=255)                                                                                                               aCGTTTTATTAGAGTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCGTTTCCATGCATAGACTTTCACTATTAGACTCCACTCGTGAtgt                      <  2:943770/143‑1 (MQ=255)                                                                                                                cGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTATTAGACTCCACTCGTGAtgtg                     >  1:471027/1‑143 (MQ=255)                                                                                                                 gTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTATTAGACTCCACTCGTGATGTGa                    >  1:551546/1‑143 (MQ=255)                                                                                                                  ttttATTAGATTTGTATTCACCGGAGTGATGTAATTAAAGATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTATTAGACTCCACTCGTGATGTGAg                   >  1:898980/1‑143 (MQ=255)                                                                                                                     tattaGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGCCTTTCTCGCTTTTAATACTATTTTTTTACATGCATAGTCTTTCTCTATTAGACTCCACTCGTGATTTGAGCGa                >  2:505483/1‑143 (MQ=255)                                                                                                                       ttaGATTTGTATTGACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCGGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTATTAGACTCCACTCGTGATGTGAGCGAGt              <  2:424807/143‑1 (MQ=255)                                                                                                                       ttaGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTATTAGACTCCACTCGTGATGTGAGCGAGt              <  2:29214/143‑1 (MQ=255)                                                                                                                                tagtCACCTGAGTGATGTGATGAAATATATGATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTATTAGACTCCACTCGTGATGTGAGCGAGTTAATTAGtc     <  2:890466/140‑1 (MQ=255)                                                                                                                                tATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGCGCTATTTGGCTGGCTTTATCGCTTTTAATACTATTCTTTTTCATGTCTATACTTTCACTATTCCACTCCTCTCGTGATGTCAGTGTGTTAATTAGta     >  2:106933/1‑142 (MQ=255)                                                                                                                                  ttCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTATTAGACTCCACTCGTGATGTGAGCGAGTTAATTAGtctc   <  2:594647/143‑1 (MQ=255)                                                                                                                                   tCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTGTTAATACTATTCTTTTCCATGCATAGACTTTCACTATTAGACTCCACTCGTGATGTGAGCGAGTTAATTAGTCTCa  >  1:352598/1‑143 (MQ=255)                                                                                                                                           |                                                                                                                                       TTGCCTGATATATTCACACCTATAAATTTAGGGCTATTTAAAAACTGCCAAGACTGGTCTAATGATTCTTTAGAAACATGTAAGATGAAGTGCTCCGGAAATAACATTGGACGTTTTATTAGATTTGTATTCACCGGAGTGATGTAATGAAATATATTATTTTTTTATTCAGAGCTATATGGCTGGCTTTATCGCTTTTAATACTATTCTTTTCCATGCATAGACTTTCACTATTAGACTCCACTCGTGATGTGAGCGAGTTAATTAGTCTCA  >  NC_000913/733630‑733902

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