Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,652,855:1	(G)6→7	100%	coding (362/846 nt)	sseA →	3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,652,849	1	.	G	96.4%	284.6 / 3.0	84	coding (356/846 nt)	sseA	3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (37/44);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAA  >  NC_000913/2652754‑2652942                                                                                                  |                                                                                                tAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑‑gggg                                                                                            <  2:1535741/100‑1 (MQ=255)  aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑‑ggggg                                                                                           <  2:254896/100‑1 (MQ=255)  aaGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑‑ggggg                                                                                           <  1:144502/100‑1 (MQ=255)       ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑‑ggggg‑ggtgg                                                                                     <  2:557399/100‑5 (MQ=255)       ccTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGAGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑‑gggcg‑g                                                                                         >  1:518856/1‑96 (MQ=255)        cTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑‑ggggg‑ggtgga                                                                                    >  2:545318/1‑95 (MQ=255)         tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGAc                                                                                    >  1:418254/1‑100 (MQ=255)         tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGAc                                                                                    >  1:978954/1‑100 (MQ=255)         tGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGAc                                                                                    >  1:1487587/1‑100 (MQ=255)          gATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑‑ggggg‑ggggggac                                                                                  >  2:2441513/1‑93 (MQ=255)            ttGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTg                                                                                 <  1:1937187/100‑1 (MQ=255)             tGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGc                                                                                <  2:1809639/100‑1 (MQ=255)              gTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCa                                                                               >  1:429404/1‑100 (MQ=255)               tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAg                                                                              >  2:1075321/1‑100 (MQ=255)               tCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAg                                                                              <  1:2576446/100‑1 (MQ=255)                cTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAgn                                                                             <  1:881813/100‑2 (MQ=255)                 tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGc                                                                            <  1:837070/100‑1 (MQ=255)                 tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGc                                                                            >  2:2471575/1‑100 (MQ=255)                 tATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGACGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGc                                                                            <  2:1149189/100‑1 (MQ=255)                  aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCt                                                                           <  2:460873/100‑1 (MQ=255)                  aTGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCt                                                                           <  2:2387531/100‑1 (MQ=255)                   tGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTg                                                                          <  2:495038/100‑1 (MQ=255)                    gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGCCTTGCAGGCTgg                                                                         >  2:1781607/1‑100 (MQ=255)                    gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTgg                                                                         >  1:2038503/1‑100 (MQ=255)                    gACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTgg                                                                         >  2:2209942/1‑100 (MQ=255)                     aCGAAGGTAATCTTGTCTCAGCCCCACGAGCATGGTAGATGCTGAGAACCTCTGGTGTAGAGAAAGTGTCGATTC‑T‑‑GGGGGAGTGGACTTGCAGGCTGGc                                                                        <  1:2169926/100‑1 (MQ=255)                         aGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGgcggcg                                                                    >  2:1773718/1‑100 (MQ=255)                           gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAgcgct                                                                  <  1:1716719/100‑2 (MQ=255)                           gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGcgcg                                                                  >  2:6884/1‑100 (MQ=255)                           gTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGcgcg                                                                  <  1:2141349/100‑1 (MQ=255)                            tAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCga                                                                 >  1:1792074/1‑100 (MQ=255)                              aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCgatg                                                               >  1:1055374/1‑100 (MQ=255)                              aTCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTGGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCgatg                                                               <  1:313946/100‑1 (MQ=255)                                cTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCgatgat                                                             >  1:1371864/1‑100 (MQ=255)                                 ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATc                                                            <  2:1116898/100‑1 (MQ=255)                                 ttttCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATc                                                            <  2:691198/100‑1 (MQ=255)                                     ctgaGCCCCACGAGCATGGTGGATGCTGCGCACCGTTGGTGTAGAGAAAGTGTCGAGTC‑TG‑GGGGG‑GTGGACTGGCAGGCTGGCAGCGCGATGATctgct                                                        <  2:808581/97‑1 (MQ=255)                                     ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                        >  2:1042438/1‑100 (MQ=255)                                     ctcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTAGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                        <  1:1441364/100‑1 (MQ=255)                                      tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAACGTGTCGAGCCGG‑‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgctg                                                       >  1:1463909/1‑100 (MQ=255)                                      tcAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TGGGGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATctgct                                                        <  1:399452/100‑1 (MQ=255)                                        aGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                     <  2:1855577/100‑1 (MQ=255)                                         gCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTa                                                    <  2:527208/100‑1 (MQ=255)                                          ccccACGAGCATGGTGGATGCTGCGTACCTTTGTTGTAGAGAAAGTGTCGATTC‑T‑‑GGGTG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGtt                                                     <  1:2397270/97‑1 (MQ=255)                                          ccccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTag                                                   >  1:2471222/1‑100 (MQ=255)                                           cccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTaga                                                  >  1:871278/1‑100 (MQ=255)                                            ccACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGGGTGGATGC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaa                                                 >  1:2548045/1‑100 (MQ=255)                                             cACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTagaag                                                <  2:2565306/100‑1 (MQ=255)                                               cGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAgaagaa                                              >  2:756591/1‑100 (MQ=255)                                                 aGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAgg                                            <  2:1494363/100‑1 (MQ=255)                                                  gCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGt                                           <  1:728817/100‑1 (MQ=255)                                                   cATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTg                                          <  2:80213/100‑1 (MQ=255)                                                   cATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTg                                          >  1:1406880/1‑100 (MQ=255)                                                    aTGTTGGATGCTGCGGACGTTGGGTGTAGGGAAAGGGTCGATTG‑TG‑GGGGG‑GTGGACTTGCAGGGTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                         <  1:2455131/100‑1 (MQ=255)                                                    aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                         >  2:845314/1‑100 (MQ=255)                                                    aTGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGc                                         <  2:1404316/100‑1 (MQ=255)                                                     tggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GGGGACTTGCAGGCTGGCAGCGCGATGATCGGCTGTTAGAAGAAGGTGc                                         >  2:1677235/1‑99 (MQ=255)                                                      ggtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAg                                       <  2:1052677/100‑1 (MQ=255)                                                       gtggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGt                                      <  2:1436035/100‑1 (MQ=255)                                                         ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                    <  1:1440240/100‑1 (MQ=255)                                                         ggATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTag                                    <  1:1310716/100‑1 (MQ=255)                                                           gatctgcGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagag                                  <  1:1956716/97‑1 (MQ=255)                                                                gcgcACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTagagagatc                             >  2:1956726/1‑95 (MQ=255)                                                                 cgcACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                            <  2:733491/100‑1 (MQ=255)                                                                 cgcACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                            <  1:1927824/100‑1 (MQ=255)                                                                  gtaCCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCg                            <  1:949925/97‑1 (MQ=255)                                                                  gcACCTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCgg                           <  2:927972/100‑1 (MQ=255)                                                                     ccTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                        <  1:1014557/100‑1 (MQ=255)                                                                     ccTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAg                        <  1:823395/100‑1 (MQ=255)                                                                      cTTTGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgg                       <  1:1524954/100‑1 (MQ=255)                                                                         tGGTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGgaga                    <  2:653700/100‑1 (MQ=255)                                                                          ggTGTAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                   >  2:1502118/1‑100 (MQ=255)                                                                          ggTGTAGAGAAAGCGTCGATTC‑TG‑GGGGG‑GTAGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                   >  2:1743416/1‑100 (MQ=255)                                                                           gtgtAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGagag                   >  2:1331619/1‑99 (MQ=255)                                                                           gtgtAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGt                  >  2:718384/1‑100 (MQ=255)                                                                            tgtgtagaAAGGGTTGATTC‑TG‑GGGGG‑GTGGAGGTGCAGGCTGGAAGCGCGATGATCTGTGGTTAGAAGAAGGTGCAGTAGAGCTGGCGGAAGGAGAGtt                 <  1:679771/95‑1 (MQ=255)                                                                            tgtAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCCGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGtt                 >  1:555283/1‑100 (MQ=255)                                                                             gtAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                >  2:1868930/1‑100 (MQ=255)                                                                             gtAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                <  1:1619173/100‑1 (MQ=255)                                                                             gtAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGttt                >  2:117825/1‑100 (MQ=255)                                                                              tAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa               >  1:211490/1‑100 (MQ=255)                                                                              tAGAGAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTa               >  2:2584039/1‑100 (MQ=255)                                                                               tcttccgatctCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTaa              <  1:1160315/90‑1 (MQ=255)                                                                               agagAAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAgatcggaagag              >  1:1552426/1‑90 (MQ=255)                                                                                gagaAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGAGAGATGAACTGCTGGTAGAAGCAGCTGGAGTAGGGCTGCCGGAAGAAGAGtata               >  1:588990/1‑98 (MQ=255)                                                                                 agaAAGTGTCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAAcg            <  2:1050787/100‑1 (MQ=255)                                                                                      gtgtCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg       <  1:153170/100‑1 (MQ=255)                                                                                      gtgtCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCcgcg       >  1:1344560/1‑100 (MQ=255)                                                                                         tCGATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGttt    <  1:712872/100‑1 (MQ=255)                                                                                           gATTC‑TG‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTaa  >  2:1814772/1‑100 (MQ=255)                                                                                            aTTC‑TG‑GGGGG‑GTGGACTTGCGGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGACGCGTTTa   >  1:2189459/1‑98 (MQ=255)                                                                                                  |                                                                                                TAAGCACCTGATTGTCTATGACGAAGGTAATCTTTTCTCAGCCCCACGAGCATGGTGGATGCTGCGCACCTTTGGTGTAGAGAAAGTGTCGATTC‑T‑‑GGGGG‑GTGGACTTGCAGGCTGGCAGCGCGATGATCTGCTGTTAGAAGAAGGTGCAGTAGAGCTGCCGGAAGGAGAGTTTAACGCCGCGTTTAA  >  NC_000913/2652754‑2652942

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