BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,035,428	A→G	noncoding (330/2905 nt)	rrlA →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,035,428	0	A	G	99.0%	677.4 / ‑1.9	289	noncoding (330/2905 nt)	rrlA	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (1/2); new base G (145/141); total (146/143)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

Only 100 of 297 total aligned reads displayed.

CCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG  >  NC_000913/4035337‑4035506

  aGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAg                                                                             >  2:128797/1‑93 (MQ=18)

   gTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGTAAAGGCGTGCGATACAGGGTGGCAGc                                                                            >  1:5681891/1‑93 (MQ=16)

   gTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGc                                                                            >  1:1263884/1‑93 (MQ=18)

   gTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGc                                                                            >  1:6775225/1‑93 (MQ=18)

   gTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGc                                                                            <  2:6647771/93‑1 (MQ=18)

   gTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGGGGCAGc                                                                            >  1:329599/1‑93 (MQ=16)

    tAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGcc                                                                           <  1:1701053/93‑1 (MQ=18)

    tAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGcc                                                                           <  2:5170808/93‑1 (MQ=18)

     aGCGGCGAGCGAACGGCGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGccc                                                                          <  2:4241885/93‑1 (MQ=255)

      gcggcgAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGGTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGcccc                                                                         >  2:3283889/1‑93 (MQ=16)

             gcgaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTAcaca                                                                  >  1:4308540/1‑93 (MQ=18)

             gcgaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTAcaca                                                                  <  2:1171362/93‑1 (MQ=21)

              cgaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACaa                                                                 <  1:5085456/93‑1 (MQ=21)

               gaACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACaaa                                                                <  1:2789440/93‑1 (MQ=21)

                 aCGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGTCCCGTACACaaaaa                                                              >  2:6537096/1‑93 (MQ=16)

                  cGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAAt                                                             >  2:7176643/1‑93 (MQ=18)

                   ggggAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCTTCTGGAAAGGAGTGCGATCCAGGGTGGCAGCCCCGCACACAAAAAtc                                                            >  1:802805/1‑92 (MQ=11)

                    gggAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGc                                                           >  2:6404607/1‑93 (MQ=18)

                    gggAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGc                                                           >  1:4708670/1‑93 (MQ=18)

                    gggAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGAGGCAGCCCCGTACACAAAAATGc                                                           >  2:7582623/1‑93 (MQ=16)

                     ggAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGca                                                          >  1:5828239/1‑93 (MQ=18)

                      gAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGcac                                                         <  2:1954656/93‑1 (MQ=21)

                       ggcagcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGcaca                                                        <  2:6967885/92‑1 (MQ=18)

                       agcagcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGcaca                                                        <  1:2382110/93‑1 (MQ=21)

                         cagcCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATg                                                      <  1:2537885/93‑1 (MQ=18)

                            cccAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCtg                                                   <  1:3723200/93‑1 (MQ=18)

                              cAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCtgtg                                                 >  2:188511/1‑93 (MQ=18)

                               agagCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGa                                                >  2:4348601/1‑93 (MQ=18)

                               agagCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGa                                                <  2:3312296/93‑1 (MQ=18)

                                gagCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAg                                               <  1:3494633/93‑1 (MQ=21)

                                gagCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAg                                               >  2:5418767/1‑93 (MQ=18)

                                gagCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAg                                               >  2:7499316/1‑93 (MQ=18)

                                 agCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGc                                              >  2:3345321/1‑93 (MQ=18)

                                 agCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGc                                              <  2:6266190/93‑1 (MQ=18)

                                    cTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCg                                           >  2:7851970/1‑93 (MQ=18)

                                    cTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCg                                           >  2:4879602/1‑93 (MQ=18)

                                    cTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCg                                           >  1:1029623/1‑93 (MQ=18)

                                     tGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGa                                          >  1:2360938/1‑93 (MQ=18)

                                     tGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGa                                          <  2:4030434/93‑1 (MQ=18)

                                     tGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGa                                          >  2:2027385/1‑93 (MQ=18)

                                      gAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGAt                                         >  1:6079675/1‑93 (MQ=18)

                                      gAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGAt                                         <  1:2290539/93‑1 (MQ=18)

                                      gAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGAt                                         >  2:2259158/1‑93 (MQ=18)

                                      gAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGAt                                         <  1:4392269/93‑1 (MQ=18)

                                       caTCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATg                                        >  2:4793989/2‑93 (MQ=18)

                                       aaTCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATg                                        >  2:7153516/1‑93 (MQ=18)

                                       aaTCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATg                                        >  1:1967427/1‑93 (MQ=18)

                                        aTCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTAGATGa                                       >  2:4151187/1‑93 (MQ=16)

                                         tCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAg                                      >  2:5787575/1‑93 (MQ=18)

                                         tCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATAAAGGGTGGCAGCCCCGGACACAAAAATGAACATGCTGTGAGCTCGATGAg                                      >  1:2012918/1‑93 (MQ=11)

                                          cAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGt                                     >  2:5935660/1‑93 (MQ=18)

                                           aGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTa                                    <  1:1776038/93‑1 (MQ=18)

                                           aGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTa                                    <  1:3941430/93‑1 (MQ=18)

                                           aGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTa                                    <  1:4156012/93‑1 (MQ=18)

                                           aGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTa                                    <  1:6964469/93‑1 (MQ=18)

                                           aGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACTAAAATGCACATGCTGTGAGCTCGATGAGTa                                    >  2:7729879/1‑93 (MQ=21)

                                            gtgtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAg                                   <  1:2623135/93‑1 (MQ=18)

                                              gtgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAggg                                 <  1:1722977/93‑1 (MQ=18)

                                               tgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGc                                >  1:7265646/1‑93 (MQ=18)

                                               tgtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGc                                <  1:7223320/93‑1 (MQ=18)

                                                gtgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCg                               >  2:5631358/1‑93 (MQ=18)

                                                 tgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCgg                              <  1:2077638/93‑1 (MQ=18)

                                                 tgtgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCgg                              <  1:576305/93‑1 (MQ=18)

                                                   tgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAgggaggga                            >  2:1514399/1‑93 (MQ=16)

                                                   tgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGa                            >  1:3228875/1‑93 (MQ=18)

                                                   tgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGa                            <  1:7408628/93‑1 (MQ=18)

                                                   tgtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGa                            <  2:7809763/93‑1 (MQ=18)

                                                    gtTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGac                           >  1:2934066/1‑93 (MQ=18)

                                                      tAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGacac                         <  1:4964081/93‑1 (MQ=18)

                                                      tAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGacac                         <  1:7831948/93‑1 (MQ=18)

                                                       aGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGCGTAGGGCGGGACACg                        <  1:7706989/93‑1 (MQ=16)

                                                       aGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACg                        <  1:701081/93‑1 (MQ=18)

                                                        gTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGt                       >  2:6736770/1‑93 (MQ=18)

                                                          ggAAGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTgg                     <  1:6212473/93‑1 (MQ=18)

                                                            aaGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTa                   <  2:174701/93‑1 (MQ=18)

                                                             aGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTAt                  <  2:3418683/93‑1 (MQ=18)

                                                             aGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTAt                  >  2:445468/1‑93 (MQ=21)

                                                             aGCGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTAt                  <  1:456357/93‑1 (MQ=18)

                                                               cGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATcc                <  2:2814857/93‑1 (MQ=18)

                                                               cGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATcc                <  1:5578024/93‑1 (MQ=18)

                                                               cGTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGGAGGGCGGGACACGTGGTATcc                >  1:1445266/1‑93 (MQ=16)

                                                                gTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAACATGCACATGCTGTGAGCTCGATGAGTAGGGGGGGAAACGTGGTATCCt               >  1:4063825/1‑93 (MQ=11)

                                                                gTCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCt               >  2:1007664/1‑93 (MQ=21)

                                                                 tCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTg              >  2:3834460/1‑93 (MQ=18)

                                                                 tCTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTg              >  2:741649/1‑93 (MQ=18)

                                                                  cTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGt             >  1:2240734/1‑93 (MQ=21)

                                                                  cTGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGt             <  1:6863351/93‑1 (MQ=18)

                                                                   tGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTc            >  2:6322750/1‑93 (MQ=18)

                                                                   tGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTc            >  2:5501975/1‑93 (MQ=18)

                                                                   tGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTc            <  1:3353329/93‑1 (MQ=18)

                                                                   tGGAAAGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTc            >  2:2994027/1‑93 (MQ=18)

                                                                       aaGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGaa        <  2:3647371/93‑1 (MQ=18)

                                                                       aaGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGaa        <  2:686145/93‑1 (MQ=18)

                                                                       aaGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGaa        >  2:2550216/1‑93 (MQ=18)

                                                                       aaGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGaa        <  1:7134939/93‑1 (MQ=18)

                                                                       aaGGCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGaa        >  1:571679/1‑93 (MQ=18)

                                                                         ggCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAata      <  2:7364416/93‑1 (MQ=18)

                                                                          gCGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAatat     <  1:3840503/93‑1 (MQ=18)

                                                                           cGTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATg    <  2:3042909/93‑1 (MQ=18)

                                                                            gTGCGATACAGGGTGGCAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgg   <  1:1063208/93‑1 (MQ=17)

CCAGTAGCGGCGAGCGAACGGGGAGCAGCCCAGAGCCTGAATCAGTGTGTGTGTTAGTGGAAGCGTCTGGAAAGGCGTGCGATACAGGGTGACAGCCCCGTACACAAAAATGCACATGCTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG  >  NC_000913/4035337‑4035506

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