BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,786,097	A→G	M124V (ATG→GTG)	b3614 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,786,097	0	A	G	99.8%	1522.1 / ‑5.9	472	M124V (ATG→GTG)	b3614	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (229/242); total (230/242)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 487 total aligned reads displayed.

CTACGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTA  >  NC_000913/3786007‑3786186

ctacgcCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTg                                                                                         <  2:7577691/93‑1 (MQ=255)

 tacgcCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTgg                                                                                        >  2:489154/1‑93 (MQ=255)

     cccGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGgcag                                                                                    >  1:6625012/1‑93 (MQ=255)

       cGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGgcagca                                                                                  >  1:1147455/1‑93 (MQ=255)

       cGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGgcagca                                                                                  >  1:3297571/1‑93 (MQ=255)

        ggAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCaa                                                                                 >  2:6934114/1‑93 (MQ=255)

         gagaTGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAg                                                                                >  2:8847531/1‑93 (MQ=255)

          agaTGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGa                                                                               <  2:6500364/93‑1 (MQ=255)

          agaTGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGa                                                                               <  1:6877695/93‑1 (MQ=255)

            aTGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATg                                                                             >  1:4590813/1‑93 (MQ=255)

            aTGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATg                                                                             >  1:197117/1‑93 (MQ=255)

            aTGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATg                                                                             >  1:7532363/1‑93 (MQ=255)

             tGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGa                                                                            >  1:4203511/1‑93 (MQ=255)

             tGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGa                                                                            <  1:5058332/93‑1 (MQ=255)

              gAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGAc                                                                           <  2:1792985/93‑1 (MQ=255)

              gAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGAc                                                                           >  2:6265087/1‑93 (MQ=255)

               agcagcGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGAcc                                                                          >  1:1519457/1‑93 (MQ=255)

                gcagcGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACct                                                                         >  2:9730998/1‑93 (MQ=255)

                gcagcGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACct                                                                         <  1:9756072/93‑1 (MQ=255)

                 cagcGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACctc                                                                        <  1:1989863/93‑1 (MQ=255)

                  agcGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACctct                                                                       >  1:9363170/1‑93 (MQ=255)

                    cgacgaAATTGAGCGCATTATTCGTAGTGCGGGCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTaa                                                                     >  1:8268577/1‑93 (MQ=255)

                      acgaAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAAcc                                                                   <  2:2341257/93‑1 (MQ=255)

                       cgcaaTAGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGAAAGATGACCTCTAAcca                                                                  <  2:8983164/90‑2 (MQ=255)

                       cgaAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTCACCt                                                                  >  2:5178774/1‑93 (MQ=255)

                        gaAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTg                                                                 >  2:670591/1‑93 (MQ=255)

                         aaaTTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGt                                                                >  2:1482772/1‑93 (MQ=255)

                          aaTTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGtt                                                               >  2:679595/1‑93 (MQ=255)

                          aaTTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGtt                                                               >  1:8636114/1‑93 (MQ=255)

                          aaTTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGtt                                                               >  1:1838299/1‑93 (MQ=255)

                          aaTTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGtt                                                               <  1:4824099/93‑1 (MQ=255)

                           aTTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGttt                                                              <  1:6666484/93‑1 (MQ=255)

                           aTTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGttt                                                              <  1:2967890/93‑1 (MQ=255)

                           aTTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGttt                                                              >  2:3620603/1‑93 (MQ=255)

                             tGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTgg                                                            <  1:3409027/93‑1 (MQ=255)

                              gAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGt                                                           <  2:3387530/93‑1 (MQ=255)

                               aGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTa                                                          >  1:3220340/1‑93 (MQ=255)

                               aGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTa                                                          <  1:3999228/93‑1 (MQ=255)

                               aGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTa                                                          <  1:1668318/93‑1 (MQ=255)

                                  gcATTATTCGTAGTGCGGTCAATAACGTGCCCTGTGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGc                                                       <  2:3547860/93‑1 (MQ=255)

                                  gcATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGc                                                       >  1:2767311/1‑93 (MQ=255)

                                   cATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCa                                                      <  1:7693635/93‑1 (MQ=255)

                                    attattCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAg                                                     >  2:5102566/1‑93 (MQ=255)

                                    attattCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAg                                                     >  1:3605518/1‑93 (MQ=255)

                                      tattCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGaa                                                   <  1:8434269/93‑1 (MQ=255)

                                      tattCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGaa                                                   >  2:5759217/1‑93 (MQ=255)

                                      tattCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGaa                                                   >  1:3532135/1‑93 (MQ=255)

                                        ttCGTAGGGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAg                                                 >  2:5827825/1‑93 (MQ=255)

                                         tCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGt                                                >  1:7473317/1‑93 (MQ=255)

                                           gTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGa                                              <  2:2106701/93‑1 (MQ=255)

                                            tAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGAt                                             >  2:8249661/1‑93 (MQ=255)

                                             aGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATg                                            <  1:7036154/93‑1 (MQ=255)

                                             aGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATg                                            >  2:184039/1‑93 (MQ=255)

                                               tGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCa                                          <  1:9552497/93‑1 (MQ=255)

                                                gCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAg                                         <  1:9013070/93‑1 (MQ=255)

                                                 cGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAgg                                        <  2:989015/93‑1 (MQ=255)

                                                    tCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGgcgc                                     >  1:638016/1‑93 (MQ=255)

                                                     cAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCt                                    >  2:1281330/1‑93 (MQ=255)

                                                     cAATAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCt                                    >  1:1215478/1‑93 (MQ=255)

                                                      aaTAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTg                                   <  1:6097702/93‑1 (MQ=255)

                                                      aaTAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTg                                   >  2:3043714/1‑93 (MQ=255)

                                                      aaTAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTg                                   >  2:7004380/1‑93 (MQ=255)

                                                       aTAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTgg                                  <  1:7346540/93‑1 (MQ=255)

                                                       aTAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTgg                                  >  1:3761264/1‑93 (MQ=255)

                                                        tAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGa                                 >  2:7388996/1‑93 (MQ=255)

                                                        tAACGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGa                                 >  2:4704123/1‑93 (MQ=255)

                                                        aaccGTGCCCTATGCCGTGGGCATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGa                                 <  2:156875/90‑1 (MQ=255)

                                                          aCGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGc                               <  1:5727782/93‑1 (MQ=255)

                                                          aCGTGCCCTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGc                               >  2:5316452/1‑93 (MQ=255)

                                                               cccTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGGGCCGGAGCGttca                          >  2:8672326/1‑91 (MQ=255)

                                                               cccTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTAc                          <  2:5397034/93‑1 (MQ=255)

                                                                ccTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACa                         >  1:3901549/1‑93 (MQ=255)

                                                                ccTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACa                         <  1:4374183/93‑1 (MQ=255)

                                                                ccTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACa                         >  1:8162382/1‑93 (MQ=255)

                                                                ccTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATCCAGGCGCTGGAGCGTTACa                         >  1:2245314/1‑93 (MQ=255)

                                                                ccTATGCCGGGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACa                         <  1:6358602/93‑1 (MQ=255)

                                                                 cTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACaa                        <  2:3171009/93‑1 (MQ=255)

                                                                 cTATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACaa                        <  2:6365065/93‑1 (MQ=255)

                                                                  tATGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAAt                       <  2:7272406/93‑1 (MQ=255)

                                                                   aTGCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATc                      <  1:2549981/93‑1 (MQ=255)

                                                                     gCCGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCtt                    >  2:399474/1‑93 (MQ=255)

                                                                     gCCGTGGGGATAAACAACCACGCGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACACTCtt                    >  1:7908472/1‑93 (MQ=255)

                                                                      ccGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCttt                   >  2:436877/1‑93 (MQ=255)

                                                                       cGTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTa                  <  2:8018898/93‑1 (MQ=255)

                                                                        gTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTAc                 <  2:4810221/93‑1 (MQ=255)

                                                                        gTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTAc                 >  2:9436596/1‑93 (MQ=255)

                                                                        gTGGGGATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTAc                 <  2:9089676/93‑1 (MQ=255)

                                                                          ggggATCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACtt               <  1:8900758/93‑1 (MQ=255)

                                                                              aTCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTc           >  2:4204499/1‑93 (MQ=255)

                                                                              aTCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTc           <  1:4478468/93‑1 (MQ=255)

                                                                              aTCAACAACCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTc           >  2:6068186/1‑93 (MQ=255)

                                                                                 aacaacCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGAc        >  2:7781308/1‑93 (MQ=255)

                                                                                 aacaacCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGAc        <  2:4635021/93‑1 (MQ=255)

                                                                                 aacaacCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGAc        <  1:1380911/93‑1 (MQ=255)

                                                                                 aacaacCACGGGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGAc        <  2:2739339/93‑1 (MQ=255)

                                                                                  acaacCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACa       >  1:5179386/1‑93 (MQ=255)

                                                                                  acaacCACGTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACa       >  2:9283874/1‑93 (MQ=255)

                                                                                   caacCACGTGGTCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAg      <  2:10033905/93‑1 (MQ=255)

                                                                                      ccACGTGGGTAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGt   >  1:2399399/1‑93 (MQ=255)

                                                                                       cacgTGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTa  >  1:1029441/1‑93 (MQ=255)

CTACGCCCGGAGATGAGCAGCGACGAAATTGAGCGCATTATTCGTAGTGCGGTCAATAACGTGCCCTATGCCGTGGGGATCAACAACCACATGGGCAGCAAGATGACCTCTAACCTGTTTGGTATGCAGAAAGTGATGCAGGCGCTGGAGCGTTACAATCTTTACTTCCTCGACAGCGTA  >  NC_000913/3786007‑3786186

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