BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	810,837	T→C	A31A (GCT→GCC)	bioC →	BioC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	810,837	0	T	C	99.6%	724.3 / NA	223	A31A (GCT→GCC)	bioC	BioC
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (125/97); minor base A (0/1); total (125/98)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 228 total aligned reads displayed.

GGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCG  >  NC_000913/810747‑810926

 gCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGAc                                                                                        <  1:9943708/93‑1 (MQ=255)

 gCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGAc                                                                                        <  2:3850904/93‑1 (MQ=255)

 gCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGAc                                                                                        <  1:537471/93‑1 (MQ=255)

   aaCGGTTAATAACCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGc                                                                                      <  2:4444929/93‑1 (MQ=255)

    aCGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGTCAGAGTGCCGACGcc                                                                                     <  2:144891/93‑1 (MQ=255)

    aCGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGcc                                                                                     >  1:4599975/1‑93 (MQ=255)

       gTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTa                                                                                  >  2:8299584/1‑93 (MQ=255)

        ttAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCTGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTAc                                                                                 >  2:746849/1‑93 (MQ=255)

        ttAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTAc                                                                                 >  1:2332880/1‑93 (MQ=255)

        ttAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTAc                                                                                 <  2:6229934/93‑1 (MQ=255)

          aataaataAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTg                                                                               <  1:1838983/93‑1 (MQ=255)

          aataaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTg                                                                               >  1:2761950/1‑93 (MQ=255)

          aataaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTg                                                                               <  2:7212457/93‑1 (MQ=255)

            taaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGc                                                                             >  2:1870408/1‑93 (MQ=255)

            taaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGc                                                                             >  1:3826070/1‑93 (MQ=255)

            taaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGc                                                                             >  1:7056000/1‑93 (MQ=255)

            taaACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGc                                                                             >  2:5407754/1‑93 (MQ=255)

             aaaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCTGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCa                                                                            <  1:6788349/93‑1 (MQ=255)

              aaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCaa                                                                           <  1:6370449/93‑1 (MQ=255)

              aaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCaa                                                                           <  2:6489465/93‑1 (MQ=255)

              aaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCaa                                                                           <  2:2390643/93‑1 (MQ=255)

              aaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCaa                                                                           <  1:3343842/93‑1 (MQ=255)

              aaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCaa                                                                           <  1:190662/93‑1 (MQ=255)

              aaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCaa                                                                           >  1:120146/1‑93 (MQ=255)

              aaCAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCACAGTGCCGACGCCCTACTGGCaa                                                                           >  2:2354972/1‑93 (MQ=255)

              aaCAAGCCATGGCAGCGGCATTTGGTCGGGCATCCGCACCTTATGAGCAACATGCAGATCTACAGCCCCAGAGTGCCGACGCCTTACTGGCaa                                                                           <  2:5236362/93‑1 (MQ=255)

                 aaGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGc                                                                        >  2:4955973/1‑93 (MQ=255)

                  aGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCt                                                                       <  1:655873/93‑1 (MQ=255)

                  aGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCt                                                                       >  2:9494291/1‑93 (MQ=255)

                   gCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCtt                                                                      >  2:5065884/1‑93 (MQ=255)

                   gCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCtt                                                                      >  2:126493/1‑93 (MQ=255)

                     cATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGGGCCGACGCCTTACTGGCAATACTTcc                                                                    >  2:8859595/1‑93 (MQ=255)

                       ttGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCcac                                                                  >  1:9379741/1‑93 (MQ=255)

                       ttGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCcac                                                                  >  1:2426299/1‑93 (MQ=255)

                       ttGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCcac                                                                  >  2:779650/1‑93 (MQ=255)

                         gCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAg                                                                >  1:845035/1‑93 (MQ=255)

                          cAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGc                                                               >  1:8205535/1‑93 (MQ=255)

                           aGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCTAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCg                                                              >  1:9716503/1‑93 (MQ=255)

                           aGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAggg                                                              <  1:7301593/93‑3 (MQ=255)

                           aGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCg                                                              >  2:3814786/1‑93 (MQ=255)

                            gCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGt                                                             >  2:2736663/1‑93 (MQ=255)

                            gCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGt                                                             <  1:8902874/93‑1 (MQ=255)

                              ggCATTTGGTCGGGCAGCTGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTaa                                                           >  1:878235/1‑93 (MQ=255)

                              ggCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTaa                                                           >  1:5589442/1‑93 (MQ=255)

                               tcctttGGCCGGCCACCCGCCCACTATGAGAAACATGCAGATCTCCAGCGCGAGGTTGCCGCCCCTTTAGGGGCAATGCTTCCACAGCGTaaa                                                          <  2:8330058/90‑1 (MQ=255)

                               gCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTaaa                                                          >  1:7354349/1‑93 (MQ=255)

                               gCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTaaa                                                          <  2:5802756/93‑1 (MQ=255)

                                cATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAAt                                                         <  1:7228608/93‑1 (MQ=255)

                                cATTTGGGCGTGCAGCCGCTCACTATGAGCAACGTGCAGATCTACAGCGCCAGAGTGCAGACGCCTTACTGGAAATGCTTCCACAGCGTAAAt                                                         <  2:4388490/93‑1 (MQ=255)

                                  tttGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTAc                                                                                 <  2:3648533/67‑1 (MQ=255)

                                  tttGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTAc                                                                                 >  1:3648533/1‑67 (MQ=255)

                                  tttGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATac                                                       >  1:9795698/1‑93 (MQ=255)

                                  tttGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATac                                                       >  2:2121773/1‑93 (MQ=255)

                                  tttGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATac                                                       <  2:7640877/93‑1 (MQ=255)

                                  tttGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATac                                                       >  1:1054885/1‑93 (MQ=255)

                                  tttGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATac                                                       >  1:1382423/1‑93 (MQ=255)

                                     ggTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACAcc                                                    <  1:9189026/93‑1 (MQ=255)

                                       tCGGGCAGCCGCACACTATGAGCAACAGGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCa                                                  <  1:7254418/93‑1 (MQ=255)

                                         gggCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACg                                                >  2:6550985/1‑93 (MQ=255)

                                               ccGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTgg                                          <  1:9467816/93‑1 (MQ=255)

                                               ccGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTga                                          >  1:7498259/1‑92 (MQ=255)

                                                 gCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGAc                                        >  1:8735667/1‑93 (MQ=255)

                                                   acacTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGAcgc                                      >  1:2500071/1‑93 (MQ=255)

                                                    cacTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGAcgcg                                     >  2:7524608/1‑93 (MQ=255)

                                                     acTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCgg                                    >  1:5922451/1‑93 (MQ=255)

                                                      cTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCggg                                   >  1:2252582/1‑93 (MQ=255)

                                                      cTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCggg                                   >  2:4004618/1‑93 (MQ=255)

                                                       tATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGt                                  >  1:4779983/1‑93 (MQ=255)

                                                       tATGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGt                                  <  1:6932192/93‑1 (MQ=255)

                                                         tGAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTtg                                <  2:1008845/93‑1 (MQ=255)

                                                          gAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTtgt                               >  1:4251014/1‑93 (MQ=255)

                                                          gAGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTtgt                               <  2:6145638/93‑1 (MQ=255)

                                                           aGCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTtgtg                              >  2:8502960/1‑93 (MQ=255)

                                                            gCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTgg                             >  2:9005516/1‑93 (MQ=255)

                                                            gCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTgg                             >  2:9183776/1‑93 (MQ=255)

                                                            gCAACATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTgg                             <  1:4943720/93‑1 (MQ=255)

                                                            gCAACATGCAGATCTACAGCGCCACAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTgg                             >  1:5087086/1‑93 (MQ=255)

                                                              aaCATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGAc                           <  2:9379740/93‑1 (MQ=255)

                                                               aCATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGAcc                          >  1:4089413/1‑93 (MQ=255)

                                                               aCATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGAcc                          <  2:9724842/93‑1 (MQ=255)

                                                                cATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACct                         >  1:4018036/1‑93 (MQ=255)

                                                                cATGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACct                         <  2:5740554/93‑1 (MQ=255)

                                                                  tGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACctgg                       >  2:6666807/1‑93 (MQ=255)

                                                                  tGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACctgg                       <  2:1788811/93‑1 (MQ=255)

                                                                  tGCAGATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACctgg                       >  2:5587388/1‑93 (MQ=255)

                                                                      gATCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACctggctgg                   >  2:2969786/1‑93 (MQ=255)

                                                                       aTCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGa                  >  1:5639880/1‑93 (MQ=255)

                                                                        tCTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGAt                 >  1:7646518/1‑93 (MQ=255)

                                                                         ctgcAGCGCCAGAGTGCCCACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATg                <  2:7989034/90‑1 (MQ=255)

                                                                         cTACAGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATg                >  1:6000152/1‑93 (MQ=255)

                                                                             aGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAgcc            <  1:6210062/93‑1 (MQ=255)

                                                                             aGCGCCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAgcc            >  2:3587086/1‑93 (MQ=255)

                                                                              gcgcCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAgccg           >  2:4628763/1‑93 (MQ=255)

                                                                              gcgcCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAgccg           <  2:3396137/93‑1 (MQ=255)

                                                                               cgcCAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAgccgc          <  1:3536712/93‑1 (MQ=255)

                                                                                 ccAGAGTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCa        <  1:7926081/93‑1 (MQ=255)

                                                                                    gagTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTg     >  2:458224/1‑93 (MQ=255)

                                                                                     agTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTgg    <  2:7862412/93‑1 (MQ=255)

                                                                                     agTGCCGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTgg    >  1:1261786/1‑93 (MQ=255)

                                                                                      gTGCCGACGCCTTACTGGCAATGCTTCCACAGAGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCAGCCCCTgca   >  1:825370/1‑91 (MQ=255)

GGCAACGGTTAATAAACAAGCCATTGCAGCGGCATTTGGTCGGGCAGCCGCACACTATGAGCAACATGCAGATCTACAGCGCCAGAGTGCTGACGCCTTACTGGCAATGCTTCCACAGCGTAAATACACCCACGTACTGGACGCGGGTTGTGGACCTGGCTGGATGAGCCGCCACTGGCG  >  NC_000913/810747‑810926

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