Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	615,932	A→G	6.1%	S592S (TCA→TCG)	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	615,932	0	A	G	6.1%	193.9 / 3.3	65	S592S (TCA→TCG)	entF	enterobactin synthase multienzyme complex component, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (31/30);  new base G (2/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

AAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCC  >  NC_000913/615790‑616076                                                                                                                                               |                                                                                                                                                 aaaTGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAAcc                                                                                                                                             <  7:377806/148‑1 (MQ=255)    tgatgCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCa                                                                                                                                          <  8:279785/148‑1 (MQ=255)               gATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCtt                                                                                                                               <  3:161760/148‑1 (MQ=255)               gATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCt                                                                                                                                >  1:268427/1‑147 (MQ=255)                    gcgTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATatc                                                                                                                          <  3:91210/148‑1 (MQ=255)                        ccgtcGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACATCCGCATCACACGGCTTATatca                                                                                                                         <  6:129839/145‑1 (MQ=255)                        ccgtcGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATatca                                                                                                                         >  5:129839/1‑145 (MQ=255)                               tGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACct                                                                                                                <  8:176314/147‑1 (MQ=255)                               tGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACct                                                                                                                >  7:176314/1‑147 (MQ=255)                                  tAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTgg                                                                                                            >  7:164892/1‑148 (MQ=255)                                      tACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTca                                                                                                                                                  >  8:238492/1‑106 (MQ=255)                                      tACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTca                                                                                                                                                  <  7:238492/106‑1 (MQ=255)                                      tACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCACAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTc                                                                                                         >  2:37053/1‑147 (MQ=255)                                         caccGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCAcc                                                                                                     >  2:389405/1‑148 (MQ=255)                                          accGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCg                                                                                                    >  3:130174/1‑148 (MQ=255)                                          accGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCg                                                                                                    >  7:264198/1‑148 (MQ=255)                                             gacgaTCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCa                                                                                                 <  8:264198/148‑1 (MQ=255)                                             gacgaTCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTACACACCCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCgg                                                                                                   >  6:52642/1‑146 (MQ=255)                                                    aaCTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGaa                                                                                          <  4:57830/148‑1 (MQ=255)                                                    aaCTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGaa                                                                                          <  4:130174/148‑1 (MQ=255)                                                     aCTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCACAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGaaa                                                                                         <  1:37053/148‑1 (MQ=255)                                                      cTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAg                                                                                        >  5:181192/1‑148 (MQ=255)                                                             gcTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAAt                                                                                 <  3:321283/148‑1 (MQ=255)                                                             gcTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAAt                                                                                 <  5:26206/148‑1 (MQ=255)                                                                ttAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGt                                                                              >  6:42929/1‑148 (MQ=255)                                                                    cGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCggg                                                                          >  7:84483/1‑148 (MQ=255)                                                                     gATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCAGACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGc                                                                         <  1:22452/148‑1 (MQ=255)                                                                      aTGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTAcc                                                                                                                 >  6:274040/1‑107 (MQ=255)                                                                      aTGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTAcc                                                                                                                 <  5:274040/107‑1 (MQ=255)                                                                              aaTTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTa                                                                >  2:105421/1‑148 (MQ=255)                                                                                 ttAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACACCCGCATCCCACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCCGGCCGAAAGGGGTAGTGGTCGGGCAGACGGcgttcg                                                             >  4:328453/1‑143 (MQ=255)                                                                                 ttAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCg                                                             >  7:96082/1‑148 (MQ=255)                                                                                  tAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATatcatc                                                                                                                       <  4:69032/89‑1 (MQ=255)                                                                                  tAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATatcatc                                                                                                                       >  3:69032/1‑89 (MQ=255)                                                                                     cAAGCCTTTGCTATAACGACCATATTACACCGAAGGGCAGTGCTCCGCTGCAAATTTGACAACCTAATCTCACGGCTTATATCATCTTTACCTCTGGCTCCAGCGGGTGGCTGAAAGGGGTAATTGTCGGGCAGTCGGCTATCGTGaa                                                         <  8:345865/148‑1 (MQ=255)                                                                                            ttGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATatcat                                                                                                                        <  4:365092/78‑1 (MQ=255)                                                                                            ttGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATatcat                                                                                                                        >  3:365092/1‑78 (MQ=255)                                                                                            ttGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTg                                                  >  7:147820/1‑148 (MQ=255)                                                                                               cTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCtt                                               >  4:218424/1‑148 (MQ=255)                                                                                                tataACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCCGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCttt                                              >  5:385509/1‑148 (MQ=255)                                                                                                 ataACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCGCAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAgg                                                                                       >  1:30880/1‑106 (MQ=255)                                                                                                 ataACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCGCAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAgg                                                                                       >  1:36911/1‑106 (MQ=255)                                                                                                 ataACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCGCAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAgg                                                                                       <  2:36911/106‑1 (MQ=255)                                                                                                 ataACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCGCAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAgg                                                                                       <  2:30880/106‑1 (MQ=255)                                                                                                 ataACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCTAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACAGCCTGCTTTg                                             >  5:265915/1‑148 (MQ=255)                                                                                                  taACGCCCAGCTTACAACGCAGGGCAGTTCGCCGCTGCAACTTTCACAACCTCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTgg                                            <  5:358268/148‑1 (MQ=255)                                                                                                   aaCGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGa                                           >  1:376481/1‑148 (MQ=255)                                                                                                   aaCGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGa                                           >  4:93988/1‑148 (MQ=255)                                                                                                   aaCGCCCCGCTTACAACGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAAAGGTCGGGCAGACGGCTATCGTCAACCGCCTGCGTTGGa                                           >  5:50004/1‑148 (MQ=255)                                                                                                    aCGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGa                                           <  8:164892/147‑1 (MQ=255)                                                                                                       ccccGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCggg                                                                                                  <  8:374473/89‑2 (MQ=255)                                                                                                       ccccGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGc                                                                                                  >  7:374473/1‑89 (MQ=255)                                                                                                               acacCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCTGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCAt                                <  6:385509/147‑1 (MQ=255)                                                                                                                         gCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTAc                     <  8:147820/148‑1 (MQ=255)                                                                                                                          cAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACa                    <  1:389405/148‑1 (MQ=255)                                                                                                                                gccgcTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGaa              <  2:206435/148‑1 (MQ=255)                                                                                                                                 ccgcTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAg             <  2:268427/148‑1 (MQ=255)                                                                                                                                     tGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGt         >  2:111201/1‑148 (MQ=255)                                                                                                                                      gCATCTTTCACAACGGCATCACACGGCTTATATCATCTTTACCTCTGGGAGCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGGGAAGATGTc        <  6:50004/148‑1 (MQ=255)                                                                                                                                       cAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTc        <  6:181192/147‑1 (MQ=255)                                                                                                                                       cAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCg       <  3:218424/148‑1 (MQ=255)                                                                                                                                        aaCTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGt      <  8:96082/148‑1 (MQ=255)                                                                                                                                         aCTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGtt     >  7:322523/1‑148 (MQ=255)                                                                                                                                            ttCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGcc  >  6:83192/1‑148 (MQ=255)                                                                                                                                            ttCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGcc  >  5:192620/1‑148 (MQ=255)                                                                                                                                             tCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGcc  <  2:203671/147‑1 (MQ=255)                                                                                                                                             tCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGCGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGcc  >  1:203671/1‑147 (MQ=255)                                                                                                                                               |                                                                                                                                                 AAATGATGCTGGAAGATGCGCGTCCGTCGCTGTTAATTACCACCGACGATCAACTGCCGCGCTTTAGCGATGTTCCCAATTTAACAAGCCTTTGCTATAACGCCCCGCTTACACCGCAGGGCAGTGCGCCGCTGCAACTTTCACAACCGCATCACACGGCTTATATCATCTTTACCTCTGGCTCCACCGGCAGGCCGAAAGGGGTAATGGTCGGGCAGACGGCTATCGTCAACCGCCTGCTTTGGATGCAAAATCATTATCCGCTTACAGGCGAAGATGTCGTTGCC  >  NC_000913/615790‑616076

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