Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	606,975	C→T	6.1%	V137I (GTC→ATC)	ybdK ←	weak gamma‑glutamyl:cysteine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	606,975	0	C	T	6.1%	210.8 / 2.7	67	V137I (GTC→ATC)	ybdK	weak gamma‑glutamyl:cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (27/35);  new base T (2/2);  total (29/38)
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.52e-01

AATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCACCGCCGCAAATTTCCAGATGATGGTCTGTG  >  NC_000913/606829‑607123                                                                                                                                                   |                                                                                                                                                     aaTATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACAt                                                                                                                                                    >  6:129168/1‑149 (MQ=255) aaTATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACAt                                                                                                                                                    >  4:311982/1‑149 (MQ=255)  atatTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATg                                                                                                                                                   <  7:130811/149‑1 (MQ=255)             tGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAATTCGTGACAAGCCGTGCAGCAAATAAATGGTGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCaaaaa                                                                                                                                        <  1:28239/149‑1 (MQ=255)                ggaggCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACgg                                                                                                                                     >  1:95850/1‑149 (MQ=255)                  aggCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTc                                                                                                                                   <  5:129168/149‑1 (MQ=255)                  aggCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTc                                                                                                                                   <  7:172989/149‑1 (MQ=255)                  aggCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTc                                                                                                                                   <  6:77361/149‑1 (MQ=255)                      aaaaCGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCct                                                                                                                               >  1:357087/1‑149 (MQ=255)                        aaCGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCctgc                                                                                                                             >  3:344050/1‑149 (MQ=255)                             gTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAACTCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAAt                                                                                                                        >  4:281353/1‑149 (MQ=255)                             gTATCCGTTCCCTACATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAAt                                                                                                                        <  7:291027/149‑1 (MQ=255)                               aTCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGc                                                                                                                                                  >  4:351492/1‑121 (MQ=255)                               aTCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGc                                                                                                                                                  <  3:351492/121‑1 (MQ=255)                                tCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGATATGCTGACCAAAAACGGTCGCCTGCTGAAt                                                                                                                        >  4:247180/1‑146 (MQ=255)                                  cGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGAt                                                                                                                   <  7:233024/149‑1 (MQ=255)                                        cTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCa                                                                                                              >  1:56388/1‑148 (MQ=255)                                        cTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCa                                                                                                              >  4:263208/1‑148 (MQ=255)                                         tGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCa                                                                                                              >  3:244508/1‑147 (MQ=255)                                             tatatGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGttt                                                                                                        >  6:251632/1‑149 (MQ=255)                                              atatGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGtttt                                                                                                       >  6:287166/1‑149 (MQ=255)                                                 tGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTc                                                                                                      >  8:142365/1‑147 (MQ=255)                                                   gCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAgc                                                                                                  <  3:281353/149‑1 (MQ=255)                                                        gcggcgGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAgcgtgcg                                                                                             <  6:237781/149‑1 (MQ=255)                                                         cggcggAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAgcgtgcgt                                                                                            <  2:357087/149‑1 (MQ=255)                                                          ggcggAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGtt                                                                                           <  8:287121/149‑1 (MQ=255)                                                             ggAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGAt                                                                                        <  4:131684/149‑1 (MQ=255)                                                              gAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATa                                                                                       <  3:311982/149‑1 (MQ=255)                                                              gAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATa                                                                                       <  5:287166/149‑1 (MQ=255)                                                                  gggCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGcgt                                                                                   <  5:251632/149‑1 (MQ=255)                                                                     cGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGATATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGcgttcg                                                                                <  3:247180/149‑1 (MQ=255)                                                                        tAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTa                                                                             >  2:224474/1‑149 (MQ=255)                                                                         acaGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTAt                                                                            <  6:115454/147‑1 (MQ=255)                                                                               cGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGT‑‑GTTATCGCATAc                                                                     <  7:113246/148‑1 (MQ=255)                                                                                ggCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATAc                                                                     >  8:357180/1‑149 (MQ=255)                                                                                ggCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATAc                                                                     <  2:95850/149‑1 (MQ=255)                                                                                 gCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATAcc                                                                    <  2:56388/149‑1 (MQ=255)                                                                                 gCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATAcc                                                                    >  8:36247/1‑149 (MQ=255)                                                                                  cacaAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATAcct                                                                   >  5:172914/1‑149 (MQ=255)                                                                                           tGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCACCATGGACATGCTGACCAAAAACGGTCGCCTGCGGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGGTATCGCATACCTCCTGACGCt                                                          >  1:205760/1‑149 (MQ=255)                                                                                               aaGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCa                                                      >  7:208535/1‑149 (MQ=255)                                                                                                  ccGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCAttt                                                   <  3:263208/149‑1 (MQ=255)                                                                                                   cGTGCACCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTCACCAAAAACGGACGCCTGCTGAAAGAGATAACCAAAGTATTCCAGAGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTc                                                  >  4:75243/1‑149 (MQ=255)                                                                                                    gTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCt                                                 <  7:142365/149‑1 (MQ=255)                                                                                                    gTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCt                                                 >  1:168606/1‑149 (MQ=255)                                                                                                    gTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGCTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCt                                                 <  4:344050/149‑1 (MQ=255)                                                                                                         gcaAATTAATTGCGTCATCGGCACTGGGCCAGCCAACATGGACATGCTGCCCTACAACGGTCGGCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAAc                                            <  6:179083/149‑1 (MQ=255)                                                                                                            aataaatGGCGTCATCGCCACTGGCGCAGCCAACATGGAGATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACggg                                         <  7:138113/149‑1 (MQ=255)                                                                                                             ataaatGGCGTCATCGCCACTGGCGCAGCCAACATGGATATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGgt                                        <  2:264165/149‑1 (MQ=255)                                                                                                               aaatGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCTTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGgtgc                                      >  6:160723/1‑149 (MQ=255)                                                                                                               aaatGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGgtgc                                      >  7:149530/1‑149 (MQ=255)                                                                                                                    gCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGcc                                 <  7:357180/149‑1 (MQ=255)                                                                                                                     cGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCa                                <  7:317420/149‑1 (MQ=255)                                                                                                                         aTCGCCACTGGCGCAGCCAACATGGATATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGCCGCTGCCATTTCTGAAACGGGTGCGTGCCAccgc                            >  8:390494/1‑149 (MQ=255)                                                                                                                         aTCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACGTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCAccgc                            >  2:328962/1‑149 (MQ=255)                                                                                                                          tCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCAccgcc                           <  6:172914/149‑1 (MQ=255)                                                                                                                          cctccACTGGCGCACCCAACATGGACATGCTGACCAAAAACGGGCGCCTGCTTAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCAccgcc                           <  2:183382/146‑1 (MQ=255)                                                                                                                            gCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCACcgccgc                         <  1:328962/149‑1 (MQ=255)                                                                                                                             ccACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCACcgccgc                         >  1:317992/1‑148 (MQ=255)                                                                                                                              cACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCAccgcc                           >  1:183382/1‑145 (MQ=255)                                                                                                                                    cgcAGCCAACATGGACATGCTGACCAAAACCGGTCGCCTGCTGAAGGAGATAACCAAAGTTTTCCAGCGTGCGGTGATAGCGTTCGTTCTCGCATTCCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCACCGCCGCAAATTTcc                 <  4:228572/149‑1 (MQ=255)                                                                                                                                    cgcAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCACCGCCGCAAATTTcc                 <  2:168606/149‑1 (MQ=255)                                                                                                                                     gcAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCACCGCCGCAAATTTCCa                <  4:244508/149‑1 (MQ=255)                                                                                                                                      cAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCACCGCCGCAAATTTCCAg               <  3:176356/149‑1 (MQ=255)                                                                                                                                       aGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCACCGCCGCAAATTTCCAga              <  8:208535/149‑1 (MQ=255)                                                                                                                                         ccAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCACCGCCGCAAATTTCCAgat             >  1:204820/1‑148 (MQ=255)                                                                                                                                              attgACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCTTTGATAGGGTTCGTTATCGCTTACCTACTGTCGTTGCCTTTTCTGATACGGGTGCGTGCCACCGCCGCAAATTTCCAGATGATGGt       <  3:75243/146‑1 (MQ=255)                                                                                                                                                   cATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACTGGCGCGTGCCCCCGCCGCAAATTTCCAGATGATGGTCtgtg  >  4:29473/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     AATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCTGCATATATGGCGACGCGGCGGAAAGGGCGATAAAGTGCGGCACAAATCGTGACAAGCCGTGCAGCAAATAAATGGCGTCATCGCCACTGGCGCAGCCAACATGGACATGCTGACCAAAAACGGTCGCCTGCTGAATGAGATAACCAAAGTTTTCCAGCGTGCGTTGATAGCGTTCGTTATCGCATACCTCCTGACGCTGCCATTTCTGAAACGGGTGCGTGCCACCGCCGCAAATTTCCAGATGATGGTCTGTG  >  NC_000913/606829‑607123

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