Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,576,077	0	C	A	44.8%	103.2 / 36.1	67	D515Y (GAT→TAT)	yddB	putative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base C (30/7);  new base A (0/30);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 1.06e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCATTAAT  >  NC_000913/1575930‑1576225                                                                                                                                                     |                                                                                                                                                     gcgcAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCa                                                                                                                                                     >  1:274071/1‑149 (MQ=255)  cgcAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCa                                                                                                                                                     >  2:967248/1‑148 (MQ=255)    cAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCAGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGaa                                                                                                                                                 <  1:898821/150‑1 (MQ=255)        aaCGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATAtt                                                                                                                                             <  2:981162/150‑1 (MQ=255)           gTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCg                                                                                                                                         >  2:722115/1‑151 (MQ=255)           gTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCg                                                                                                                                         <  2:1081868/151‑1 (MQ=255)               ttACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCAt                                                                                                                                     >  2:123278/1‑151 (MQ=255)                    gATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTa                                                                                                                                 >  2:523289/1‑150 (MQ=255)                        ttCTGCTGCAACCCCATTGCAATTTCATCGTAATAGGGCTTTTTAAATTCCTAATAACGCGTAAGATTTTTATTCCCTGATACCGATTCCGCCCAGATATGGCGAATATCACGTAATCCC‑ATAAAAAAAATATTCCCGCCATAGTAACgg                                                                                                                             <  1:929661/150‑1 (MQ=255)                         tCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGtt                                                                                                                           >  1:909112/1‑151 (MQ=255)                                 aaccccattcccattTCATCTTTAAAAGGCTTTTTAAAATCCTAATAACGCGTCAGATTTTTTTTCCCTAACACCGATCCCGCCCAGCTATTGCGAATATCACAAAACCCC‑AAATAAAGAATTTTCCCGCCATAGTAACGGTTATaa                                                                                                                       <  1:459385/144‑1 (MQ=255)                                 aaccccTTTGCCTGTTCATCTTAAAAGGGCTTTTTAAAACCCTAAAAACGCGTCAGATTTTTATTCCCTAATACCGATCCCCCCCAGCTATTGCGAATATCACAAAACCCA‑ATATAAAAAATTTTCCCGCCATAGTAACGGTTATaa                                                                                                                       <  1:631069/144‑1 (MQ=255)                                 aaccccATTGCCTTTTTATCTTTATAGGTTTTTTACAACTCCTGATAACGCGTCAGATTTTTTTTCCCTAATACCGATCCCGCCCAGCTATTGCGAATATCACATAACCCC‑ATAAAAAAAATATTCCCGCCATAGTAACGGTtata                                                                                                                        <  1:13665/143‑1 (MQ=255)                                 aaccccATGGCCAGTTCATCGTTATAAGGCGTTTTAAAATCCTGATAACGCGTAAGATTTTTATTCCCTAATACCGATCCCGCCCAGCTATTGCGGATATCACGTAACCCC‑ATATAAAGAATATTCCCGCCATAGTAACGGTTATAACct                                                                                                                    <  2:274094/147‑1 (MQ=255)                                    cccATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACctgctg                                                                                                                >  2:478673/1‑151 (MQ=255)                                     ccATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGt                                                                                                               >  1:998548/1‑151 (MQ=255)                                           ccttttcttcTTTATAAGGCTTTTTCAAATCTTGAAACCGCGTAAGATTTTTATTCCCGGATCCCGATCCCGCCCTGTTATTGCGAATATCACTTAACCCA‑ATACCAAAATTATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATc                                                                                                           <  1:1103324/145‑1 (MQ=255)                                           caatttcttcGTTATAAGTCTTTTTAAAATCCTGATAACGCGTAAGAGTTTTATTCCCTGATCCCGTTCCCGTCCAGTTATGGCGAATACCACGTAATCCA‑ATATAAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAAt                                                                                                            <  2:284592/144‑1 (MQ=255)                                               ttCATGTTTAAAGGGCGTTTAAAACTCTTGATAACGGTTAAGATTTTTTTTCCCTGATCCCTTTCCCGCCCAGCTATTGGGGATATCCCGTAACCCC‑ATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATc                                                                                                           <  2:36627/145‑1 (MQ=255)                                               ttCATCTTTATAAGGCTTTTTCAAATCCTGAAACCGCTTCAGTTTTTTATTACCTAATACCGATTCCGTCCAGTTTTTGGGGATATCACGTAACCCC‑AAATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCAt                                                                                                         <  2:681594/147‑1 (MQ=255)                                               ttCATCTTTAAAGGGTTTTTTAAAATCCTGATACCGCGTAAGTTTTTTTTTCCCTGATCCCGATTCCGCCCAGCTTTTGGGAATACCACGTAACCCC‑ATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGa                                                                                                      <  2:985181/150‑1 (MQ=255)                                                tCATCTTTATAGGGGTTTTTAAAATCTTAATAACGCGTAAGATTTTTTTTCCCTGATCCCGTTCCCGCCCAGTTATTGGGAATATCACTAAACCCA‑ATAAAAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGa                                                                                                      <  2:699043/149‑1 (MQ=255)                                                  aTCGTTATAAGGCTTTTTAAAACCTTGATAACGCGTCAGAGTTTTTTTACCTGAACCCGATTCCGTCCAGCTTTGGCGGATATCACTTAACCCC‑ATACAAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAgg                                                                                                    <  2:900047/149‑1 (MQ=255)                                                    cGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTg                                                                                                >  1:261169/1‑151 (MQ=255)                                                       tataAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGAtt                                                                                             >  2:1115777/1‑151 (MQ=255)                                                        ataAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGAtt                                                                                             >  2:405998/1‑150 (MQ=255)                                                          aaggttttttAAAATCCTGATAACGCGTCAGTTTTTTTTTCCCTGATCCCGATCCCGTCCAGCTATTGCGGATATCACGAAACCCA‑ATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTa                                                                                            <  2:531236/143‑1 (MQ=255)                                                          aaGGCTTTTTCAAACCTTGATACCGGGTAAGATTTTTATTCCCTGATCCCGATCCCGCCCAGTTATTGCGAATATCACTTACCCCC‑ATAAAAAAAATTTCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGc                                                                                          <  2:163088/151‑1 (MQ=255)                                                          aaGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGc                                                                                          >  2:70503/1‑151 (MQ=255)                                                           aGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGc                                                                                          >  1:778642/1‑150 (MQ=255)                                                              cGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaa                                                                                       >  2:90233/1‑150 (MQ=255)                                                                 tttCAAATCCTGATAACGCGTAAGAGTTTTATTACCTGATCCCGATCCCGTCCAGCTATTGCGGATATCACGTAACCCA‑ATAAAAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAat                                                                                    <  2:1035472/150‑1 (MQ=255)                                                                    cAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATcc                                                                                >  1:878406/1‑151 (MQ=255)                                                                     aaaTCCTGATAACGCGTCAGATTTTTATTCCCTGATCCCGTTTCCGCCCAGCTATTGCGAATATCACTTAACCCA‑ATAAAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATccc                                                                               <  2:963066/151‑1 (MQ=255)                                                                       aTCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCAt                                                                             >  1:823506/1‑151 (MQ=255)                                                                         ccgaaTACCGCGTAAGATTTTTTTTACCTAACCCCATTCCCGTCCAGTAATTGCGAATATCAATTAACCCA‑ATAAAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTc                                                                           <  1:866624/147‑1 (MQ=255)                                                                           tGATAACGCTTAAAATTTTTATTACCTGATCCCTTTTCCTTCCTGTTTTGGCGAATATCACGAAACCCA‑ATAACAAAATTTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTc                                                                           <  1:468894/149‑1 (MQ=255)                                                                             aaaaCGCGTTGGTGTTTTTTTCCCGGATCCCGTTTCGGCCCGGCTATTGGGGATACCACAAACTCCC‑ATAAAAAGATTATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCg                                                                         <  2:249199/147‑1 (MQ=255)                                                                             aTACCGCGTAAGTGTTTTATTCCCTGACCCCTTTTCCGCCCAGCTATGGCGAATATCACTTAACCCA‑ATAAAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCa                                                                              <  1:751175/144‑1 (MQ=255)                                                                             aTAACGCGTCAGATTTTTATTCCCTGATACCTATCCCGTCCAGCTATTGGGGATATCACGAAACCCC‑ATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTc                                                                       <  2:247982/151‑1 (MQ=255)                                                                                  gcgtaagttttttATTCCCTGATCCCAATCCCTCCCAGCTATGGCGAATATAACAAAACCCC‑ATAAAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATa                                                                    <  1:209471/140‑1 (MQ=255)                                                                                         agTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCggg                                                           >  1:897610/1‑151 (MQ=255)                                                                                                ttccttattactgtttCCGTCCGTTTATGACGAATATCACATACTCCC‑ATTTAAAAAATTTGCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGgag                                                        <  2:1069981/133‑1 (MQ=255)                                                                                                   ccTGATCCCAATCCCTTCCACCTATTGCGAATACCACATAACCCC‑ATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGatat                                                    <  2:1004460/148‑1 (MQ=255)                                                                                                          cggttcccGTCAAGTTATTGCGAATATCACGTAATCCC‑ATAAAAAGAATTTCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTg                                           <  1:1129781/144‑1 (MQ=255)                                                                                                           ctaTTCCGCCCAGCTATTGCGGATATCACGTAATCCC‑ATATAAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTg                                           <  2:907495/147‑1 (MQ=255)                                                                                                              ttCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGAc                                         >  2:128399/1‑148 (MQ=255)                                                                                                                ccGCCCAGTTATTGCGGTTATCACGTAATCCA‑ATATAAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGAt                                     <  2:148421/150‑1 (MQ=255)                                                                                                                    ccactttttGCGAATACCCCATAACCCC‑ATAAAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGt                                  <  2:422724/142‑1 (MQ=255)                                                                                                                     cAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgt                               >  1:438187/1‑151 (MQ=255)                                                                                                                      aGCTATTGCGAATATCACAAAATCCC‑ATAACAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTg                                <  1:692558/149‑1 (MQ=255)                                                                                                                       gCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtc                              >  1:239182/1‑150 (MQ=255)                                                                                                                           ttGGGGATACCACGAAATCCA‑AAATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtc                              <  2:331925/146‑1 (MQ=255)                                                                                                                             gCGGATATCACTTAACCCC‑ATATAAAGAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtca                             <  2:1028496/145‑1 (MQ=255)                                                                                                                             gCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcatagtc                        >  1:392596/1‑150 (MQ=255)                                                                                                                              gggattatacagtaTCCCAATATAAAGATTTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcatagt                         <  1:875936/136‑1 (MQ=255)                                                                                                                              cGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcatagtcat                      >  1:89310/1‑151 (MQ=255)                                                                                                                                gATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcatagtcata                     >  2:702448/1‑150 (MQ=255)                                                                                                                                gATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATAc                    >  2:127603/1‑151 (MQ=255)                                                                                                                                     cACGTAACCCA‑AAAAAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCa                <  2:38035/150‑1 (MQ=255)                                                                                                                                      aCGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACg              >  2:843879/1‑151 (MQ=255)                                                                                                                                          aaTCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCgg          >  2:469191/1‑151 (MQ=255)                                                                                                                                             ccc‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGc        >  1:902884/1‑150 (MQ=255)                                                                                                                                             ccc‑ATAGAAAGAATTTCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCa       <  1:1106465/151‑1 (MQ=255)                                                                                                                                                tatatCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGCTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCAtt     >  1:1030104/2‑151 (MQ=255)                                                                                                                                                  tatCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCATTa    >  1:845310/1‑150 (MQ=255)                                                                                                                                                   atCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCATTAAt  >  2:376592/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCGCAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCC‑ATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCATTAAT  >  NC_000913/1575930‑1576225

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