BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,594,566	A→G	G515G (GGT→GGC)	ydeK ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,594,566	0	A	G	98.8%	518.4 / ‑3.4	170	G515G (GGT→GGC)	ydeK	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (88/80); total (90/80)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 173 total aligned reads displayed.

ATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCA  >  NC_000913/1594476‑1594654

aTTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGcc                                                                                        <  1:3682875/93‑1 (MQ=255)

 ttCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCa                                                                                       >  1:1280612/1‑93 (MQ=255)

   ccAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCAtt                                                                                     <  2:1866494/93‑1 (MQ=255)

    cAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTa                                                                                    >  2:3167614/1‑93 (MQ=255)

      gagaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACt                                                                                  >  1:4929813/1‑93 (MQ=255)

      gagaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACt                                                                                  >  1:3751634/1‑93 (MQ=255)

      gagaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACt                                                                                  >  1:3366/1‑93 (MQ=255)

        gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                                                                <  1:6931039/93‑1 (MQ=255)

        gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                                                                >  1:4638777/1‑93 (MQ=255)

        gaCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGa                                                                                >  1:1622466/1‑93 (MQ=255)

         gcTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGAc                                                                               <  1:2426671/92‑1 (MQ=255)

         aCTGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGAc                                                                               >  2:3245713/1‑93 (MQ=255)

           tGTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAt                                                                             >  1:361825/1‑93 (MQ=255)

            gTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                                                                            <  1:7370099/93‑1 (MQ=255)

            gTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                                                                            <  2:4413134/93‑1 (MQ=255)

            gTCCGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACAtt                                                                            <  2:1335914/93‑1 (MQ=255)

              ccGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTa                                                                          >  2:3099291/1‑93 (MQ=255)

              ccGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTa                                                                          <  1:3530793/93‑1 (MQ=255)

               cGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTAc                                                                         >  1:4856855/1‑93 (MQ=255)

               cGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTAc                                                                         <  2:4224781/93‑1 (MQ=255)

               cGTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTAc                                                                         <  1:4117189/93‑1 (MQ=255)

                gTTGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACg                                                                        >  1:137039/1‑93 (MQ=255)

                 ttGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                                                                       >  2:1176853/1‑93 (MQ=255)

                 ttGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGa                                                                       <  2:217182/93‑1 (MQ=255)

                  tGAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGAt                                                                      >  1:1822348/1‑93 (MQ=255)

                   gAAATATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATa                                                                     <  1:7539859/93‑1 (MQ=255)

                     aaTATTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATAcc                                                                   <  1:5757032/93‑1 (MQ=255)

                      atatTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGCTACCg                                                                  >  1:1527794/1‑93 (MQ=255)

                       tatTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGc                                                                 >  2:1632923/1‑93 (MQ=255)

                       tatTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGc                                                                 >  2:1632922/1‑93 (MQ=255)

                        atTAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCt                                                                >  2:7174948/1‑93 (MQ=255)

                         ttAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCtc                                                               >  2:4457587/1‑92 (MQ=255)

                         ttAACTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGAGGCGCCATTACTGACATTTACGATCCCGCtc                                                               >  1:654826/1‑92 (MQ=255)

                           aaCTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTAcc                                                             >  1:4651709/1‑93 (MQ=255)

                            aCTCTCCCCTTACCGGGGGCATTGCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCa                                                            <  1:3812221/93‑1 (MQ=255)

                            aCTATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCa                                                            >  1:2586710/1‑93 (MQ=255)

                              tATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaa                                                          >  2:5738421/1‑93 (MQ=255)

                              tATCCCCTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCaaa                                                          >  1:7405782/1‑93 (MQ=255)

                                 ccccTTACCGGGGGCATTTCCTCCGATATATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAAtt                                                       >  1:3876070/1‑93 (MQ=255)

                                 ccccTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGACTTAAGGGTGGCACCATTACTGACATTTACGATACCGCTTCCCAAAtt                                                       >  1:3578356/1‑93 (MQ=255)

                                  cccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTa                                                      >  2:238024/1‑93 (MQ=255)

                                   ccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAc                                                     <  1:1139609/93‑1 (MQ=255)

                                   ccTTACCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTAc                                                     >  2:2407714/1‑93 (MQ=255)

                                       aCCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGAc                                                 <  1:5432587/93‑1 (MQ=255)

                                       aCCGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACA‑TTACGATACCGCTACCAAAATTATCGACa                                                <  1:7274640/93‑1 (MQ=255)

                                        ccGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa                                                >  1:4603665/1‑93 (MQ=255)

                                        ccGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa                                                >  1:714341/1‑93 (MQ=255)

                                        ccGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa                                                <  2:1915204/93‑1 (MQ=255)

                                        ccGGGGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACa                                                >  2:2867609/1‑93 (MQ=255)

                                           ggggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATAACGCTACCAAAATTACCGACATTg                                             >  2:4280981/3‑93 (MQ=255)

                                            gggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGa                                            >  2:4192130/2‑93 (MQ=255)

                                            gggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGa                                            <  1:7771765/92‑1 (MQ=255)

                                              gCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGagc                                          >  2:5106281/1‑91 (MQ=255)

                                              gCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATa                                          <  2:361681/93‑1 (MQ=255)

                                                aTTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAcc                                        >  1:2115227/1‑93 (MQ=255)

                                                aTTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAcc                                        >  2:249262/1‑93 (MQ=255)

                                                aTTGCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAcc                                        <  2:4803069/93‑1 (MQ=255)

                                                 tttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAccc                                       >  2:1790650/1‑93 (MQ=255)

                                                 tttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAccc                                       >  1:1076713/1‑93 (MQ=255)

                                                  ttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTATCGACATTGATACCca                                      >  1:2012100/1‑93 (MQ=255)

                                                     ctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAg                                   >  1:7091346/1‑93 (MQ=255)

                                                     ctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAg                                   >  2:1389658/1‑93 (MQ=255)

                                                      tccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGt                                  >  2:7215142/1‑93 (MQ=255)

                                                      tccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGt                                  <  1:697033/93‑1 (MQ=255)

                                                        cGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTAt                                <  1:766768/93‑1 (MQ=255)

                                                        cGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTAt                                <  2:3572294/93‑1 (MQ=255)

                                                        cGATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTAt                                <  1:5848326/93‑1 (MQ=255)

                                                         gATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATc                               >  2:1586245/1‑93 (MQ=255)

                                                         gATAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATc                               >  2:5757726/1‑93 (MQ=255)

                                                          aTAAATCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCg                              >  1:5273125/1‑93 (MQ=255)

                                                            aaaTCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGaa                            <  1:1189346/93‑1 (MQ=255)

                                                               tCCATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACga                         <  2:282509/93‑1 (MQ=255)

                                                                 cATAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatg                       >  1:7107942/1‑93 (MQ=255)

                                                                  aTAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatga                      <  1:2536970/93‑1 (MQ=255)

                                                                  aTAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatga                      <  2:5305534/93‑1 (MQ=255)

                                                                   tAGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACgatgat                     >  2:2047377/1‑93 (MQ=255)

                                                                    aGCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGatgatg                    <  1:3689992/93‑1 (MQ=255)

                                                                     gCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGt                   >  2:4930970/1‑93 (MQ=255)

                                                                     gCCCGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGt                   <  2:4903713/93‑1 (MQ=255)

                                                                     gCCCGTTGAATTAAGGGCGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGAATCGAACGATGATGt                   >  1:2056020/1‑93 (MQ=255)

                                                                      cccGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTc                  <  1:2208268/93‑1 (MQ=255)

                                                                      cccGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTc                  <  2:5590518/93‑1 (MQ=255)

                                                                        cGTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGc                >  1:2998416/1‑93 (MQ=255)

                                                                         gTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGca               <  1:6574928/93‑1 (MQ=255)

                                                                         gTTGAATTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGca               >  1:1791954/1‑93 (MQ=255)

                                                                          ttGAATTAAGGATGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcat              >  1:6184204/1‑93 (MQ=255)

                                                                            gAATTAAGGGCGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcc            >  2:429575/1‑93 (MQ=255)

                                                                             aaTTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcca           <  2:3063060/93‑1 (MQ=255)

                                                                              aTTAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatccat          >  2:6009246/1‑93 (MQ=255)

                                                                               ttAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatccatc         <  2:3838778/93‑1 (MQ=255)

                                                                                tAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCa        <  1:5237807/93‑1 (MQ=255)

                                                                                tAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCa        <  2:2613188/93‑1 (MQ=255)

                                                                                tAAGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCa        >  2:1493308/1‑93 (MQ=255)

                                                                                  aGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAgg      <  1:3533524/93‑1 (MQ=255)

                                                                                  aGGGTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAgg      <  2:6729615/93‑1 (MQ=255)

                                                                                   gggTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGt     >  1:769968/1‑93 (MQ=255)

                                                                                   aggTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGt     <  2:6782199/92‑1 (MQ=255)

                                                                                    ggTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGtt    <  1:6450448/93‑1 (MQ=255)

                                                                                    ggTGGCGCCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGtt    >  2:3804701/1‑93 (MQ=255)

                                                                                      tGGCGCCATTACTGACATGTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCa  <  2:3986781/93‑1 (MQ=255)

ATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCA  >  NC_000913/1594476‑1594654

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