BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	294,990	T→C	H444R (CAC→CGC)	intF ←	putative phage integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	294,990	0	T	C	98.9%	860.8 / ‑2.6	269	H444R (CAC→CGC)	intF	putative phage integrase
Reads supporting (aligned to +/- strand): ref base T (1/2); new base C (148/118); total (149/120)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 278 total aligned reads displayed.

TATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATA  >  NC_000913/294900‑295079

tatGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCgcgc                                                                                         >  2:4142469/1‑93 (MQ=255)

tatGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCgcgc                                                                                         <  2:489930/93‑1 (MQ=255)

 atGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGcc                                                                                        >  1:3252442/1‑93 (MQ=255)

 atGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGcc                                                                                        >  2:4125711/1‑93 (MQ=255)

  tGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCa                                                                                       >  1:2532268/1‑93 (MQ=255)

  tGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCa                                                                                       >  1:2159652/1‑93 (MQ=255)

     tGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCAttt                                                                                    <  2:4449074/93‑1 (MQ=255)

     tGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCAttt                                                                                    <  2:2582573/93‑1 (MQ=255)

      gATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCAtttt                                                                                   >  2:4533767/1‑93 (MQ=255)

       aTTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTc                                                                                  <  1:3608590/93‑1 (MQ=255)

       aTTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTc                                                                                  <  2:2416391/93‑1 (MQ=255)

         ttttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGt                                                                                <  1:673078/93‑1 (MQ=255)

          tttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTa                                                                               >  2:4026088/1‑93 (MQ=255)

           ttAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTaa                                                                              <  2:2058190/93‑1 (MQ=255)

            tAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAAc                                                                             <  2:2845536/93‑1 (MQ=255)

              aTGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAg                                                                           >  2:815691/1‑93 (MQ=255)

               tGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGa                                                                          <  1:599950/93‑1 (MQ=255)

               tGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGa                                                                          <  1:1042971/93‑1 (MQ=255)

                ggAATCAACGCATATCAACGTTGTTTTTTATTGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCTCCATTTTCGTAACAGAt                                                                         <  1:1691791/93‑1 (MQ=255)

                ggAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGAt                                                                         >  2:1387131/1‑93 (MQ=255)

                 gAATCAACTCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATc                                                                        <  1:1139977/93‑1 (MQ=255)

                 gAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATc                                                                        <  2:2205233/93‑1 (MQ=255)

                   aTCACCGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTa                                                                      <  2:2997861/93‑1 (MQ=255)

                    tCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTaa                                                                     <  2:1960660/93‑1 (MQ=255)

                     cAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAAc                                                                    <  1:4726990/93‑1 (MQ=255)

                       aCGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACgg                                                                  >  1:583601/1‑93 (MQ=255)

                        cGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGga                                                                 <  2:1141709/93‑1 (MQ=255)

                           atatCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGgacga                                                              >  1:1467147/1‑93 (MQ=255)

                            tatCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGgacgac                                                             >  1:3814905/1‑93 (MQ=255)

                            tatCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGgacgac                                                             >  2:2345332/1‑93 (MQ=255)

                             atCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGacgacg                                                            >  2:3443729/1‑93 (MQ=255)

                             atCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGacgacg                                                            >  2:2523490/1‑93 (MQ=255)

                               cAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGc                                                          >  1:724906/1‑93 (MQ=255)

                                aaCGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCg                                                         >  1:2340865/1‑93 (MQ=255)

                                 aCGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGa                                                        >  1:1014988/1‑93 (MQ=255)

                                  cGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGAt                                                       <  2:3128789/93‑1 (MQ=255)

                                  cGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGAt                                                       <  1:13304/93‑1 (MQ=255)

                                   gttgttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATa                                                      >  2:118062/1‑93 (MQ=255)

                                     tgttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAgt                                                    >  2:3553017/1‑93 (MQ=255)

                                     tgttTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATATGACGGACGACGGCGATAgt                                                    >  1:2602812/1‑93 (MQ=255)

                                       ttttttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAgtgt                                                  >  2:3833576/1‑93 (MQ=255)

                                       ttttttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAgtgt                                                  >  2:721053/1‑93 (MQ=255)

                                        tttttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGtt                                                 <  2:4602786/93‑1 (MQ=255)

                                        tttttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGtt                                                 <  1:395624/93‑1 (MQ=255)

                                         ttttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGttt                                                <  1:862172/93‑1 (MQ=255)

                                          tttATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGtttt                                               >  1:4326440/1‑93 (MQ=255)

                                             aTGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATGGTGTTTTTCg                                            >  2:2083594/1‑93 (MQ=255)

                                              tGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCgg                                           >  2:2762233/1‑93 (MQ=255)

                                                  aaTACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAag                                       >  2:4055027/1‑93 (MQ=255)

                                                   aTACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAaga                                      <  2:4490552/93‑1 (MQ=255)

                                                    tACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAagag                                     >  1:2066456/1‑93 (MQ=255)

                                                     aCCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgc                                    <  2:1220253/93‑1 (MQ=255)

                                                     aCCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgc                                    >  2:3310269/1‑93 (MQ=255)

                                                      ccTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcg                                   >  2:1280336/1‑93 (MQ=255)

                                                      ccTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcg                                   >  1:1209839/1‑93 (MQ=255)

                                                       cTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgcgc                                  <  1:4137383/93‑1 (MQ=255)

                                                       cTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAgagc                                  >  1:2806164/1‑93 (MQ=255)

                                                        tGCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGAGCt                                 >  1:4070976/1‑93 (MQ=255)

                                                         gCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTg                                <  2:1814261/93‑1 (MQ=255)

                                                         gCTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTg                                >  2:4675417/1‑93 (MQ=255)

                                                          cTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTgg                               <  2:1909105/93‑1 (MQ=255)

                                                          cTTCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTgg                               <  2:2935620/93‑1 (MQ=255)

                                                           ttCATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTggg                              >  1:4400821/1‑93 (MQ=255)

                                                             cATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGtt                            >  1:64338/1‑93 (MQ=255)

                                                             cATTTAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGtt                            <  1:2907560/93‑1 (MQ=255)

                                                               tttAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTtg                          >  1:263377/1‑93 (MQ=255)

                                                                ttAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTtgt                         <  2:1767873/93‑1 (MQ=255)

                                                                ttAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTtgt                         <  2:88176/93‑1 (MQ=255)

                                                                ttAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTtgt                         <  2:574367/93‑1 (MQ=255)

                                                                ttAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCCGGGTTtgt                         >  1:3149030/1‑93 (MQ=255)

                                                                 tAAGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTtgtg                        <  1:460000/93‑1 (MQ=255)

                                                                  aaGATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTgtgt                       <  1:4712926/93‑1 (MQ=255)

                                                                    gATCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTcc                     >  2:4315282/1‑93 (MQ=255)

                                                                     aTCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTccc                    <  2:2154225/93‑1 (MQ=255)

                                                                     aTCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTccc                    >  1:128756/1‑93 (MQ=255)

                                                                      tctaTGTCTCAAATTTCTCGCGCCATCTTCGTAACAGATCTAACGGACGACGGCGACAGCGATTTTCGGCAAGAGCGTTGGGTAGGTGTCCCa                   >  2:4218436/4‑93 (MQ=255)

                                                                      tCCATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCa                   >  2:913312/1‑93 (MQ=255)

                                                                       ccATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCAt                  >  2:4730020/1‑93 (MQ=255)

                                                                       ccATGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGAGCTGGGTTTGTGTTCCAt                  >  2:1481970/1‑93 (MQ=255)

                                                                       ccATGTCTCAATTTTCTCGCGCCATTTTCGGAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCAt                  <  1:3439586/93‑1 (MQ=255)

                                                                         aTGTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATaa                >  2:1255694/1‑93 (MQ=255)

                                                                           gTCTCAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAtt              <  1:776382/93‑1 (MQ=255)

                                                                            tctcAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAAttt             <  1:3306146/93‑1 (MQ=255)

                                                                              tcAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGa           <  2:3253417/93‑1 (MQ=255)

                                                                              tcAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGa           >  2:102276/1‑93 (MQ=255)

                                                                              tcAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGa           >  1:4383469/1‑93 (MQ=255)

                                                                              tcAATTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGa           <  2:2484969/93‑1 (MQ=255)

                                                                                aaTTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGc         >  2:4509452/1‑93 (MQ=255)

                                                                                 aTTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGca        >  2:3702160/1‑93 (MQ=255)

                                                                                 aTTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGca        >  1:1595181/1‑93 (MQ=255)

                                                                                 aTTTTCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGca        >  1:3572601/1‑93 (MQ=255)

                                                                                  ttttCTCGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCACGAGCGCTGGGTTTGTGTCCCATAATTTGAGaaa       >  2:1852829/1‑93 (MQ=255)

                                                                                     tctcGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaacaa    <  1:190666/93‑1 (MQ=255)

                                                                                     tctcGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaacaa    <  2:2624851/93‑1 (MQ=255)

                                                                                     tctcGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaacaa    <  1:3134045/93‑1 (MQ=255)

                                                                                     tctcGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGcaacaa    <  1:4005131/93‑1 (MQ=255)

                                                                                       tcGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATa  >  2:4592747/1‑93 (MQ=255)

                                                                                       tcGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATa  <  1:94642/93‑1 (MQ=255)

                                                                                       tcGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATa  <  1:904042/93‑1 (MQ=255)

                                                                                       tcGCGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATa  >  2:1653112/1‑93 (MQ=255)

TATGTTGATTTTTAATGGAATCAACGCATATCAACGTTGTTTTTTATGGTAATACCTGCTTCATTTAAGATCCATGTCTCAATTTTCTCGTGCCATTTTCGTAACAGATCTAACGGACGACGGCGATAGTGTTTTTCGGCAAGAGCGCTGGGTTTGTGTCCCATAATTTGAGCAACAATA  >  NC_000913/294900‑295079

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: `‑`