BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	509,575	T→C	P343P (CCA→CCG)	ybaR ←	putative ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	509,575	0	T	C	82.9%	532.7 / 114.3	252	P343P (CCA→CCG)	ybaR	putative ATPase
Reads supporting (aligned to +/- strand): ref base T (24/19); new base C (100/108); total (125/127)
Fisher's exact test for biased strand distribution p-value = 4.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 253 total aligned reads displayed.

GTCAGCATCGCTTCATCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCTGGCTGCACTTCTGCCAGAGGCACGCTTTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTAACTTTTCCAGC  >  NC_000913/509483‑509665

  cAGCATCGCTTCATCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCTgg                                                                                          <  2:6753770/93‑1 (MQ=255)

  cAGCATCGCTTCATCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCgg                                                                                          >  2:3007684/1‑93 (MQ=255)

   aGCATCGCTTCATCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGc                                                                                         <  2:944568/93‑1 (MQ=255)

   aGCATCGCTTCATCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGc                                                                                         >  2:3260940/1‑93 (MQ=255)

   aGCATCGCTTCATCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGc                                                                                         >  1:3179079/1‑93 (MQ=255)

   aGCATCGCTTCATCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGc                                                                                         <  2:4630605/93‑1 (MQ=255)

      aTCGCTTCATCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGc                                                                                      <  2:256114/93‑1 (MQ=255)

      aTCGCTTCATCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGc                                                                                      >  2:2714908/1‑93 (MQ=255)

         gCTTCATCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACt                                                                                   >  1:810638/1‑93 (MQ=255)

         gCTTCATCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACt                                                                                   >  1:1273138/1‑93 (MQ=255)

           ttCATCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTc                                                                                 <  2:5354599/93‑1 (MQ=255)

           ttCATCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTc                                                                                 <  1:3447163/93‑1 (MQ=255)

            tCATCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCt                                                                                <  2:5018631/93‑1 (MQ=255)

              aTCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCTGGCTGCACTTCTGc                                                                              >  1:1618582/1‑93 (MQ=255)

               tCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCTGGCTGCACTTCTGcc                                                                             >  2:1284018/1‑93 (MQ=255)

                  aGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCaga                                                                          <  1:6839066/93‑1 (MQ=255)

                  aGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCaga                                                                          <  2:5982879/93‑1 (MQ=255)

                  aGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCaga                                                                          <  1:6501387/93‑1 (MQ=255)

                   gCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCagag                                                                         >  1:1849578/1‑93 (MQ=255)

                   gCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCagag                                                                         <  1:4510924/93‑1 (MQ=255)

                   gCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTAAGGCGCAGCAACATACCCGGCTGCACTTCTGCCagag                                                                         >  1:3985013/1‑93 (MQ=255)

                     cATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGc                                                                       <  2:37317/93‑1 (MQ=255)

                     cATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGc                                                                       >  1:1391892/1‑93 (MQ=255)

                       tGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCAc                                                                     >  2:5791495/1‑93 (MQ=255)

                       tGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCAc                                                                     >  2:7181343/1‑93 (MQ=255)

                         cTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGc                                                                   >  2:3709787/1‑93 (MQ=255)

                              cccTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGCttttt                                                              >  1:4441449/1‑93 (MQ=255)

                               ccTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGCtttttc                                                             >  1:5619565/1‑92 (MQ=255)

                               ccTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGCtttttc                                                             >  1:6123099/1‑92 (MQ=255)

                                 tGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGCtttttcac                                                           <  1:78904/93‑4 (MQ=255)

                                  gggTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGCtttttcacc                                                          <  2:5330415/93‑5 (MQ=255)

                                   ggTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCTGGCTGCACTTCTGCCAGAGGCACGATTTTTTCAcc                                                         <  2:4021886/93‑1 (MQ=255)

                                   ggTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCt                                                        <  1:3946909/93‑1 (MQ=255)

                                   ggTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCt                                                        <  1:4183001/93‑1 (MQ=255)

                                     tAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCTGGCTGCACTTCTGCCAGAGGCACGCTTTTTTCACCtt                                                       <  1:2402399/93‑1 (MQ=255)

                                     tAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCTGGCTGCACTTCTGCCAGAGGCACGCTTTTTTCACCtt                                                       <  1:5502078/93‑1 (MQ=255)

                                     tAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTc                                                      >  2:1246977/1‑93 (MQ=255)

                                     tAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTc                                                      <  2:6291413/93‑1 (MQ=255)

                                          tcGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAg                                                 >  2:3208262/1‑93 (MQ=255)

                                           cGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGt                                                <  2:5132740/93‑1 (MQ=255)

                                            gCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTa                                               >  2:3086232/1‑93 (MQ=255)

                                            gCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTa                                               >  2:5452464/1‑93 (MQ=255)

                                             ccATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTaa                                              <  1:5212072/93‑1 (MQ=255)

                                              cATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCTGGCTGCACTTCTGCCAGAGGCACGCTTTTTTCACCTTCGTCAGTaa                                              <  1:5884499/93‑1 (MQ=255)

                                              cATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCTGGCTGCACTTCTGCCAGAGGCACGCTTTTTTCACCTTCGTCAGTaa                                              <  1:3850226/93‑1 (MQ=255)

                                              cATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAAc                                             >  2:6630608/1‑93 (MQ=255)

                                              cATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAAc                                             <  1:2345839/93‑1 (MQ=255)

                                              cATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAAc                                             >  2:2148329/1‑93 (MQ=255)

                                               aTCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAAcc                                            >  1:6336599/1‑93 (MQ=255)

                                               aTCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAAcc                                            <  2:1147023/93‑1 (MQ=255)

                                                tCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCa                                           >  2:3146666/1‑93 (MQ=255)

                                                tCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCa                                           >  1:5128100/1‑93 (MQ=255)

                                                tCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCa                                           <  2:5202061/93‑1 (MQ=255)

                                                 cTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAg                                          <  1:3502253/93‑1 (MQ=255)

                                                 cTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAg                                          <  2:557601/93‑1 (MQ=255)

                                                  tACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAgg                                         >  2:4510083/1‑93 (MQ=255)

                                                  tACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAgg                                         <  2:2091033/93‑1 (MQ=255)

                                                    ccGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCTGGCTGCACTTCTGCCAGAGGCACGCTTTTTTCACCTTCGTCAGTAACCAGGc                                        >  2:480210/1‑93 (MQ=255)

                                                     cGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGt                                      <  2:665536/93‑1 (MQ=255)

                                                      ggCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTg                                     >  2:1050996/1‑93 (MQ=255)

                                                        cACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGcc                                   <  2:6463213/93‑1 (MQ=255)

                                                         aCGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCg                                  <  1:3249095/93‑1 (MQ=255)

                                                         aCGCGATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCg                                  >  2:6075426/1‑93 (MQ=255)

                                                             gATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTcgg                              >  1:3554127/1‑93 (MQ=255)

                                                             gATCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTcgg                              >  2:1354747/1‑93 (MQ=255)

                                                              aTCGCCGGTCGTCAGGCGCAGCAACATACCTGGCTGCACTTCTGCCAGAGGCACGCTTTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTcgg                              >  1:4283960/1‑93 (MQ=255)

                                                              aTCGCCGGTCGTCAGGCGCAGCAACATACCTGGCTGCACTTCTGCCAGAGGCACGCTTTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTcgg                              >  1:2654204/1‑93 (MQ=255)

                                                               tCGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTcggcg                            >  1:1057801/1‑93 (MQ=255)

                                                                cGCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTcggcgg                           <  2:6845515/93‑1 (MQ=255)

                                                                 gCCGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCggg                          <  1:2233127/93‑1 (MQ=255)

                                                                  ctggTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCgggg                         >  1:6599058/3‑93 (MQ=255)

                                                                  ccGGTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCgggg                         <  1:4702921/93‑1 (MQ=255)

                                                                   gggTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGt                        <  1:3721467/92‑1 (MQ=255)

                                                                    ggTCGTCAGGCGCAGCAACATACCTGGCTGCACTTCTGCCAGAGGCACGCTTTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGt                        >  1:1110088/1‑93 (MQ=255)

                                                                    ggTCGTCAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGtt                       >  1:2272286/1‑93 (MQ=255)

                                                                       cgtcAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTaaa                    <  2:1780137/93‑1 (MQ=255)

                                                                       cgtcAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTaaa                    <  2:4931032/93‑1 (MQ=255)

                                                                        gtcAGGCGCAGCAACATACCTGGCTGCACTTCTGCCAGAGGCACGCTTTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTaaa                    <  1:2198134/93‑1 (MQ=255)

                                                                        gtcAGGCGCAGCAACATACCTGGCTGCACTTCTGCCAGAGGCACGCTTTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTaaa                    >  1:4631742/1‑93 (MQ=255)

                                                                        gtcAGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGGTAAAt                   >  1:718497/1‑93 (MQ=255)

                                                                          aaGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGGGAAATCg                 >  1:1704750/2‑93 (MQ=255)

                                                                           aGGCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGa                <  1:623713/93‑1 (MQ=255)

                                                                            ggCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAg               <  2:713752/93‑1 (MQ=255)

                                                                            ggCGCAGCAACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAg               <  1:7216661/93‑1 (MQ=255)

                                                                               gcagcaACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTaa            >  2:2417285/1‑93 (MQ=255)

                                                                               gcagcaACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTaa            <  2:4368295/93‑1 (MQ=255)

                                                                                  gcaACATACCTGGCTGCACTTCTGCCAGAGGCACGCTTTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTAACt          >  2:1898624/1‑93 (MQ=255)

                                                                                  gcaACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTAACtt         <  1:1978151/93‑1 (MQ=255)

                                                                                  gcaACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTAACtt         >  2:1481406/1‑93 (MQ=255)

                                                                                   caACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTAACttt        >  2:3489156/1‑93 (MQ=255)

                                                                                   caACATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTAACttt        <  2:3946290/93‑1 (MQ=255)

                                                                                    aaCATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTAACtttt       >  2:2338265/1‑93 (MQ=255)

                                                                                     aCATACCTGGCTGCACTTCTGCCAGAGGCACGCTTTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTAACtttt       >  2:4295094/1‑93 (MQ=255)

                                                                                     aCATACCTGGCTGCACTTCTGCCAGAGGCACGCTTTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTAACtttt       >  1:2743768/1‑93 (MQ=255)

                                                                                     aCATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTAACTTTTc      <  2:3388188/93‑1 (MQ=255)

                                                                                      cATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTAACTTTTcc     <  1:6932380/93‑1 (MQ=255)

                                                                                      cATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTAACTTTTcc     >  2:1591501/1‑93 (MQ=255)

                                                                                      cATACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGGTAAATCGAGTAACTTTTcc     >  2:4749383/1‑93 (MQ=255)

                                                                                       aTACCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTAACTTTTCCa    >  1:5076188/1‑93 (MQ=255)

                                                                                         aCCCGGCTGCACTTCTGCCAGAGGCACGC‑TTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTAACTTTTCCAgc  <  1:215736/93‑1 (MQ=255)

GTCAGCATCGCTTCATCCAGCCATGCTTCGCCCTGGGTAATCTCGCCATCTACCGGCACGCGATCGCCGGTCGTCAGGCGCAGCAACATACCTGGCTGCACTTCTGCCAGAGGCACGCTTTTTTCACCTTCGTCAGTAACCAGGCGTGCCGTCGGCGGGGTTAAATCGAGTAACTTTTCCAGC  >  NC_000913/509483‑509665

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