Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,851,204	A→C	7.2%	T115T (ACA→ACC)	ansA →	L‑asparaginase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,851,204	0	A	C	7.2%	252.8 / 8.4	81	T115T (ACA→ACC)	ansA	L‑asparaginase 1
Reads supporting (aligned to +/- strand):  ref base A (34/41);  new base C (3/3);  total (37/44)
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

GAAGACTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGC  >  NC_000913/1851055‑1851350                                                                                                                                                      |                                                                                                                                                   gAAGACTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACa                                                                                                                                                    >  1:723672/1‑150 (MQ=255) gAAGACTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACa                                                                                                                                                    <  2:873957/150‑1 (MQ=255) gAAGACTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACa                                                                                                                                                    >  2:1117681/1‑150 (MQ=255)   aGACTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAgg                                                                                                                                                  >  2:660570/1‑150 (MQ=255)     aCTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGt                                                                                                                                                <  1:819434/150‑1 (MQ=255)      cTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTc                                                                                                                                               <  2:709860/150‑1 (MQ=255)         gcaggaTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTTTCATTCTGCACGGAACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTcaca                                                                                                                                            <  2:810105/150‑1 (MQ=255)          cagcaTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACaa                                                                                                                                           >  1:531850/1‑150 (MQ=255)            gcaTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACCCTGTCACaatg                                                                                                                                         <  1:666510/150‑3 (MQ=255)            gcaTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACCCTGTCACaatg                                                                                                                                         <  1:1177070/150‑3 (MQ=255)            gcaTATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAAt                                                                                                                                         >  1:412337/1‑150 (MQ=255)                atTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACCCTGTCACaatgaccg                                                                                                                                     <  1:86681/150‑7 (MQ=255)                 ttGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgc                                                                                                                                    >  1:524885/1‑150 (MQ=255)                  tGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgct                                                                                                                                   <  1:55670/150‑1 (MQ=255)                   gCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctg                                                                                                                                  >  1:1347955/1‑150 (MQ=255)                     tGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggc                                                                                                                                <  2:1293255/150‑1 (MQ=255)                     tGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggc                                                                                                                                >  1:1246334/1‑150 (MQ=255)                       aaGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggctg                                                                                                                              <  2:179662/150‑1 (MQ=255)                       aaGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCgctggctg                                                                                                                              >  2:642508/1‑150 (MQ=255)                       aaGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCAAGGAACCGACAAGATGGCGTATACCGCCTCTGAGCTGTCGTTCATGAACGAGAATCTCGGAAAACCGGTCAATGTGACAGGGTCAAAAATCCCgctggctg                                                                                                                              >  2:1410878/1‑150 (MQ=255)                         gATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAg                                                                                                                            <  1:162889/150‑1 (MQ=255)                           tatTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGtt                                                                                                                          >  1:983437/1‑150 (MQ=255)                                 aGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGctc                                                                                                                    <  2:154337/150‑1 (MQ=255)                                      aCTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACg                                                                                                               <  2:1170570/150‑1 (MQ=255)                                      aCTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACg                                                                                                               <  2:1165003/150‑1 (MQ=255)                                                  atgatgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATc                                                                                                   <  1:1392729/150‑1 (MQ=255)                                                      tgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATctgct                                                                                               <  1:657384/150‑1 (MQ=255)                                                      tgGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATctgct                                                                                               >  1:37168/1‑150 (MQ=255)                                                       ggTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATctgctg                                                                                              >  2:1459469/1‑150 (MQ=255)                                                       ggTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATctgctg                                                                                              <  2:531850/150‑1 (MQ=255)                                                       ggTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATctgctg                                                                                              <  1:962110/150‑1 (MQ=255)                                                        gTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGa                                                                                             <  2:802533/150‑1 (MQ=255)                                                         ttttGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGaa                                                                                            <  1:1002752/150‑1 (MQ=255)                                                             gTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGGTCATGCTCGAGAAGCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTCCGCTCGGACGGACAAATTAATCTGCTGAATGCg                                                                                        >  2:1486950/1‑150 (MQ=255)                                                                aTTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTg                                                                                     <  1:482117/150‑1 (MQ=255)                                                                aTTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTg                                                                                     <  1:65944/150‑1 (MQ=255)                                                                aTTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTg                                                                                     >  1:1424365/1‑150 (MQ=255)                                                                aTTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTg                                                                                     <  1:737352/150‑1 (MQ=255)                                                                  tCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTa                                                                                   <  1:837939/150‑1 (MQ=255)                                                                    tGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACg                                                                                 >  2:1074532/1‑150 (MQ=255)                                                                    tGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACg                                                                                 >  1:865122/1‑150 (MQ=255)                                                                    tGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACg                                                                                 <  1:1052650/150‑1 (MQ=255)                                                                      cACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGtt                                                                               >  1:520077/1‑150 (MQ=255)                                                                       aCGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTg                                                                              <  1:457271/150‑1 (MQ=255)                                                                           cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTgcggc                                                                          >  1:653764/1‑150 (MQ=255)                                                                           cACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTgcggc                                                                          <  1:1074532/150‑1 (MQ=255)                                                                               gACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAAt                                                                      <  1:1424523/150‑1 (MQ=255)                                                                               gACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAAt                                                                      <  2:685480/150‑1 (MQ=255)                                                                               gACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAAt                                                                      >  1:1460674/1‑150 (MQ=255)                                                                                 cacGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTa                                                                    <  1:1117681/150‑1 (MQ=255)                                                                                    gATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATcc                                                                 <  2:1408728/150‑1 (MQ=255)                                                                                          gTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCaa                                                           <  1:656424/150‑1 (MQ=255)                                                                                            ataCCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACg                                                         <  1:1459469/150‑1 (MQ=255)                                                                                               ccGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAg                                                      <  2:291018/150‑1 (MQ=255)                                                                                               ccGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAg                                                      >  2:1414415/1‑150 (MQ=255)                                                                                                cGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGt                                                     >  1:650190/1‑150 (MQ=255)                                                                                                cGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGt                                                     >  1:1523518/1‑150 (MQ=255)                                                                                                           tGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCtttt                                          <  1:1533042/150‑1 (MQ=255)                                                                                                            gTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCttttt                                         >  1:587669/1‑150 (MQ=255)                                                                                                              ggTTCATGCTCGAGAATCTCGGTATACCGGTCATTTTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGATGAATGCGTTGTACGTTGCGGCTAATTATCCGATCAACGAAGTAACGCTCTTTTTCa                                       <  2:412337/149‑1 (MQ=255)                                                                                                              cGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCa                                       <  1:1270243/150‑1 (MQ=255)                                                                                                                ttCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAAt                                     >  1:921196/1‑150 (MQ=255)                                                                                                                  cATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTATTTCAATaa                                   <  2:983437/150‑1 (MQ=255)                                                                                                                    tGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAAcc                                 >  1:1481253/1‑150 (MQ=255)                                                                                                                     gCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCg                                >  1:344791/1‑150 (MQ=255)                                                                                                                      cTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGa                               >  1:173072/1‑150 (MQ=255)                                                                                                                       tCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATATCCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGAt                              >  2:1463958/1‑150 (MQ=255)                                                                                                                              tctcGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGAt                              >  2:51396/1‑143 (MQ=255)                                                                                                                              tctcGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGAt                              <  1:51396/143‑1 (MQ=255)                                                                                                                                tcGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATcgcg                     >  1:907129/1‑150 (MQ=255)                                                                                                                                       aCCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCg              >  2:934575/1‑150 (MQ=255)                                                                                                                                       aCCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCg              <  2:1080913/150‑1 (MQ=255)                                                                                                                                          ggTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCAc           >  2:60833/1‑150 (MQ=255)                                                                                                                                          ggTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCAc           >  1:1442690/1‑150 (MQ=255)                                                                                                                                           ggatttGTGACCCTGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACt          >  1:216205/5‑150 (MQ=255)                                                                                                                                           ggatttGTGACCCTGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACt          >  2:286465/5‑150 (MQ=255)                                                                                                                                           ggatttGTGACCCTGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACt          >  1:3636/5‑150 (MQ=255)                                                                                                                                           gTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACt          <  1:642508/150‑1 (MQ=255)                                                                                                                                              aTTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTAcc       <  2:1460674/150‑1 (MQ=255)                                                                                                                                               ttGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACAAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCa      <  2:524885/150‑1 (MQ=255)                                                                                                                                                 gtgACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCaaa    >  2:1468082/1‑150 (MQ=255)                                                                                                                                                  tgACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAg   <  2:482118/150‑1 (MQ=255)                                                                                                                                                   gACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGc  <  2:921196/150‑1 (MQ=255)                                                                                                                                                   gACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGc  >  1:449739/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                   GAAGACTGGCAGCATATTGCTGAAGATATTAAAGCGCACTATGACGACTATGATGGTTTTGTCATTCTGCACGGCACCGACACGATGGCGTATACCGCCTCTGCGCTGTCGTTCATGCTCGAGAATCTCGGTAAACCGGTCATTGTGACAGGGTCACAAATCCCGCTGGCTGAGTTACGCTCTGACGGACAAATTAATCTGCTGAATGCGTTGTACGTTGCGGCGAATTATCCGATCAACGAAGTAACGCTCTTTTTCAATAACCGATTGTATCGCGGCAACCGCACTACCAAAGC  >  NC_000913/1851055‑1851350

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