Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	617,560	T→C	5.3%	L1135P (CTG→CCG)	entF →	enterobactin synthase multienzyme complex component, ATP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	617,560	0	T	C	5.3%	230.1 / 3.0	75	L1135P (CTG→CCG)	entF	enterobactin synthase multienzyme complex component, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (36/35);  new base C (2/2);  total (38/37)
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.29e-01

GTTATCTCGATCCACAATGGTCGATTATCGGCATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAACTGGCAGGAAAAAGAAGCTAATGGTCTG  >  NC_000913/617413‑617702                                                                                                                                                    |                                                                                                                                               gTTAGCTCGATCCACAATGGTCGATAATCGGCATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGGCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACTGCATGGCCCTTATTACctgctg                                                                                                                                               <  1:154787/149‑1 (MQ=255)       tcGATCCACAATGGTCGATTATCGGCATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCCGGGGTAt                                                                                                                                         >  7:174076/1‑149 (MQ=255)         gATCCACAATGGTCGATTATCGGCATTCAGTCACCGCGCCGCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTAtt                                                                                                                                        <  8:395833/148‑1 (MQ=255)         gATCCACAATGGTCGATTATCGGCATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTa                                                                                                                                          <  6:167255/146‑1 (MQ=255)         gATCCACAATGGTCGATTATCGGCATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTAtt                                                                                                                                        <  6:258033/148‑1 (MQ=255)          aTCCACAATGGTCGATTATCGGCATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTcc                                                                                                                                      >  1:105890/1‑149 (MQ=255)                aaTGGTCGATTATCGGCATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGc                                                                                                                                <  4:384900/149‑1 (MQ=255)                 aTGGTCGATTATCGGCATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACAGCATGTCCCTTATTACCTGCTGGGGTATTCCCTTGGCg                                                                                                                               <  8:310749/149‑1 (MQ=255)                       gATTATCGGCATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGc                                                                                                                         >  2:117276/1‑149 (MQ=255)                           aTCGGCATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTgg                                                                                                                      <  2:105890/148‑1 (MQ=255)                            tCGGCATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGgcg                                                                                                                    <  8:434637/149‑1 (MQ=255)                            tCGGCATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGgcg                                                                                                                    >  8:193754/1‑149 (MQ=255)                            tCGGCATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGgcg                                                                                                                    >  3:277250/1‑149 (MQ=255)                               gCATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAg                                                                                                                 >  4:388648/1‑149 (MQ=255)                                cATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAgg                                                                                                                <  1:117276/149‑1 (MQ=255)                                cATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAgg                                                                                                                <  7:193754/149‑1 (MQ=255)                                 aTTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAggg                                                                                                               <  1:438672/149‑1 (MQ=255)                                  ttCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTCTTCCATTGGCGGTACGCTGGCGCAgat                                                                                                               >  2:228985/1‑146 (MQ=255)                                         aCCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGAGCAGGGTATTgcgg                                                                                                       >  6:39967/1‑149 (MQ=255)                                             cgcACCAATGGCCCCATGCAGACGGCGGCTAACCTTTATGAAGTCTGCTAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGTCCCTTATTACCTGCTGGGGTATTCCCTTGGCgg                                                                                                                              <  8:403564/122‑1 (MQ=255)                                                  cAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGc                                                                                              <  6:190388/149‑1 (MQ=255)                                                          ccATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTgcggc                                                                                                      <  7:344802/133‑1 (MQ=255)                                                          ccATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTgcggc                                                                                                      >  8:344802/1‑133 (MQ=255)                                                          ccATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGt                                                                                            >  7:82678/1‑143 (MQ=255)                                                          ccATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGt                                                                                            <  8:82678/143‑1 (MQ=255)                                                           cATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTg                                                                                     >  3:64898/1‑149 (MQ=255)                                                             tGCAGACGGCGGCAAACCTGGCTGAAGTCTGCGAAGCGCAGCTGGCAACGGTACATGAACAACAACCGCATGGCCCTTATTCCCTGCTGGGGTATTCCCTTGGCGGTACGCGGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGGGGc                                                                                   <  7:124304/149‑1 (MQ=255)                                                             tGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGAc                                                                                                 <  1:107249/135‑1 (MQ=255)                                                             tGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGAc                                                                                                 >  2:107249/1‑135 (MQ=255)                                                                aGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGGAACGTTACTTGAACCACAACCGCATGACCCTTATTACCTGCTGGGGTATCCCCTTGGCGGTACGCTGGCGCAGTGTATTGCGGCGCGACTGCGTGCCCGTGGCGaa                                                                                >  5:342045/1‑149 (MQ=255)                                                                aGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTGCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTg                                                                                           <  4:148053/138‑1 (MQ=255)                                                                aGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTg                                                                                           >  3:148053/1‑138 (MQ=255)                                                                aGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGaa                                                                                >  5:109446/1‑149 (MQ=255)                                                                  aCGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCAGGCAACGGTACATGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTc                                                                                                                                       >  8:364060/1‑92 (MQ=255)                                                                     gcggcAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGt                                                                           <  2:370510/149‑1 (MQ=255)                                                                     gcggcAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCCGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGt                                                                           <  8:174076/149‑1 (MQ=255)                                                                      cggcAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTg                                                                          <  6:109446/149‑1 (MQ=255)                                                                      cggcAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTg                                                                          <  4:89430/149‑1 (MQ=255)                                                                           aaCCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCAtt                                                                     >  7:35180/1‑149 (MQ=255)                                                                               tGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCtt                                                                 >  4:267335/1‑149 (MQ=255)                                                                                ggATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGccc                                                                                        <  4:355874/125‑1 (MQ=255)                                                                                ggATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGccc                                                                                        >  3:355874/1‑125 (MQ=255)                                                                                 gATGAAGTCTTCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTg                                                                <  8:134156/148‑1 (MQ=255)                                                                                       gTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCt                                                         <  3:388648/149‑1 (MQ=255)                                                                                             gAAGAGAATCTGGCAACGTTACTTGAACAACAAACGCATGAACCTTATTAACTGCTGGTGTATCCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATa                                                    <  6:342045/148‑1 (MQ=255)                                                                                               aGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCt                                                 >  8:193287/1‑149 (MQ=255)                                                                                                gcgcATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTg                                                <  1:109964/149‑1 (MQ=255)                                                                                                  gcATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCCGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCt                                                 <  2:42949/146‑1 (MQ=255)                                                                                                  gcATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCCGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCt                                                 >  1:42949/1‑146 (MQ=255)                                                                                                    aTCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGgccg                                            >  4:3161/1‑149 (MQ=255)                                                                                                       tGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGgccgc                                           >  2:288293/1‑147 (MQ=255)                                                                                                        ggCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCAt                                                                      <  2:175759/119‑1 (MQ=255)                                                                                                        ggCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCAt                                                                      >  1:175759/1‑119 (MQ=255)                                                                                                           aaCGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGaaa                                     >  7:36203/1‑149 (MQ=255)                                                                                                           aaCGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGaaa                                     >  4:210752/1‑149 (MQ=255)                                                                                                            aCGTTACTTGAACAACAACCGCATGGCCCTTATTACCTTCTGGGGTATTCCCTTGGCGTTACGCTGGCGCAGGTTATTGCGGCTCGACTGCGTTACCGGGTCGAACAGTTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAAc                                    <  8:275531/149‑1 (MQ=255)                                                                                                            aCGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAAc                                    <  7:193287/149‑1 (MQ=255)                                                                                                                 aCTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCaa                                >  8:108278/1‑148 (MQ=255)                                                                                                                   ttGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAAc                             >  6:264916/1‑149 (MQ=255)                                                                                                                   ttGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAAc                             >  8:227649/1‑149 (MQ=255)                                                                                                                   ttGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAAc                             >  8:268473/1‑149 (MQ=255)                                                                                                                   ttGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGACTTGc                                                          >  3:315380/1‑120 (MQ=255)                                                                                                                   ttGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGACTTGc                                                          <  4:315380/120‑1 (MQ=255)                                                                                                                   ttGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCAGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAAc                             >  1:200667/1‑149 (MQ=255)                                                                                                                        caacaacCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAACTgg                          >  8:292784/1‑147 (MQ=255)                                                                                                                        caacaacCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAACTGGCa                        <  4:64898/149‑1 (MQ=255)                                                                                                                             acCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAACTGGCAGGaaa                   <  7:268473/149‑1 (MQ=255)                                                                                                                              ccGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAACTGGCAGGaaa                   >  8:392449/1‑148 (MQ=255)                                                                                                                                gCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTTCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGTATACCTGGCCGCCAGAAACGCAAAACTGGCAGGAAAaag                <  5:39967/149‑1 (MQ=255)                                                                                                                                  aTGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAACTGGCAGGAAAaag                >  6:72789/1‑147 (MQ=255)                                                                                                                                       ccTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAACTGGCAGGAAAAAGAAGCTaa         <  6:81413/149‑1 (MQ=255)                                                                                                                                       ccTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAACTGGCAGGAAAAAGAAGCTaa         >  4:82878/1‑149 (MQ=255)                                                                                                                                         ttattaCCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAACTGGCAGGAAAAAGAAGCTAATg       >  3:55844/1‑149 (MQ=255)                                                                                                                                            ttaCCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAACTGGCAGGAAAAAGAAGCTAATGGTc    <  3:3161/149‑1 (MQ=255)                                                                                                                                               ccTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTg                                                                <  6:151491/86‑1 (MQ=255)                                                                                                                                               ccTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTg                                                                >  5:151491/1‑86 (MQ=255)                                                                                                                                                 tgctgGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAACTGGCAGGAAAAAGAAGCTAATGGTCTg  >  3:50669/1‑146 (MQ=255)                                                                                                                                                 tgctgGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAACTGGCAGGAAAAAGAAGCTAATGGTCTg  <  4:50669/146‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                               GTTATCTCGATCCACAATGGTCGATTATCGGCATTCAGTCACCGCGCCCCAATGGCCCCATGCAGACGGCGGCAAACCTGGATGAAGTCTGCGAAGCGCATCTGGCAACGTTACTTGAACAACAACCGCATGGCCCTTATTACCTGCTGGGGTATTCCCTTGGCGGTACGCTGGCGCAGGGTATTGCGGCGCGACTGCGTGCCCGTGGCGAACAGGTGGCATTTCTTGGCTTGCTGGATACCTGGCCGCCAGAAACGCAAAACTGGCAGGAAAAAGAAGCTAATGGTCTG  >  NC_000913/617413‑617702

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