BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,774,233	T→C	G98G (GGT→GGC)	yfjX →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,774,233	0	T	C	98.6%	430.2 / ‑2.5	140	G98G (GGT→GGC)	yfjX	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/1); new base C (78/60); total (79/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

Only 100 of 151 total aligned reads displayed.

TTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC  >  NC_000913/2774143‑2774322

ttACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGgcgc                                                                                         <  1:4393679/93‑1 (MQ=255)

ttACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGgcgc                                                                                         <  1:2606462/93‑1 (MQ=255)

  acacGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCgg                                                                                       >  2:1258251/1‑93 (MQ=255)

   cacGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGa                                                                                      >  1:3738595/1‑93 (MQ=255)

    acGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGaa                                                                                     <  1:331835/93‑1 (MQ=255)

     cGCTCAGCTATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGaaa                                                                                    <  2:4413516/93‑1 (MQ=255)

      gCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAAt                                                                                   <  2:861318/93‑1 (MQ=255)

      gCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAAt                                                                                   <  2:1839141/93‑1 (MQ=255)

      gCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAAt                                                                                   <  1:1384726/93‑1 (MQ=255)

       cTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATg                                                                                  <  2:1889791/93‑1 (MQ=255)

       cTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATg                                                                                  <  2:463565/93‑1 (MQ=255)

        tCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGa                                                                                 <  2:2019783/93‑1 (MQ=255)

          aGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAgc                                                                               <  1:4491519/93‑1 (MQ=255)

           gCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAgcg                                                                              >  2:3282991/1‑93 (MQ=255)

           gCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAgcg                                                                              >  1:3410625/1‑93 (MQ=255)

            cAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAgcgc                                                                             >  2:3431717/1‑93 (MQ=255)

            cAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAgcgc                                                                             >  2:3204713/1‑93 (MQ=255)

            cAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAgcgc                                                                             >  2:493508/1‑93 (MQ=255)

             aaTGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGcgcg                                                                            <  2:3300896/93‑1 (MQ=255)

             aaTGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGcgcg                                                                            >  1:1076548/1‑93 (MQ=255)

             aaTGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGcgcg                                                                            <  2:4372310/93‑1 (MQ=255)

                  cGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGc                                                                       >  1:2092301/1‑93 (MQ=255)

                   ggAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGcc                                                                      <  1:339991/93‑1 (MQ=255)

                    gAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCg                                                                     >  2:3349606/1‑93 (MQ=255)

                    gAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCg                                                                     >  1:3190105/1‑93 (MQ=255)

                    gAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATAGCGCGGAAATGAGCGCGGAAGCCg                                                                     >  2:3524300/1‑93 (MQ=255)

                     aGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGc                                                                    >  1:3169176/1‑93 (MQ=255)

                      gCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCa                                                                   >  1:4414345/1‑93 (MQ=255)

                      gCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCa                                                                   >  1:1999288/1‑93 (MQ=255)

                        gTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAgg                                                                 <  1:3689204/93‑1 (MQ=255)

                         tttATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGt                                                                >  2:3105262/1‑93 (MQ=255)

                         tttATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAACCCGCAGGt                                                                >  1:1264762/1‑93 (MQ=255)

                          ttgtGGCTCCTGATGCTGACGGTGATGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTa                                                               <  1:946061/90‑1 (MQ=255)

                          ttATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTa                                                               <  1:106130/93‑1 (MQ=255)

                           tATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTAt                                                              <  2:1253416/93‑1 (MQ=255)

                           tATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTAt                                                              >  2:1332345/1‑93 (MQ=255)

                            aTGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATc                                                             >  1:1442767/1‑93 (MQ=255)

                              ggCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGc                                                           >  1:116949/1‑93 (MQ=255)

                              ggCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGc                                                           <  2:429807/93‑1 (MQ=255)

                               gCTCCCGATGCTGACGGTGACTATAAATGGCATGTGCTCAACGGCATGAATGGCAATGTCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCt                                                          <  2:10472/93‑1 (MQ=255)

                                cTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTg                                                         >  2:4131745/1‑93 (MQ=255)

                                cTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTg                                                         <  1:2953963/93‑1 (MQ=255)

                                cTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTg                                                         <  2:313969/93‑1 (MQ=255)

                                 tCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGt                                                        >  2:1043841/1‑93 (MQ=255)

                                  ccTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATTAGCGCGGAAGCCGCAGGTATCGCTGTc                                                       >  2:1644440/1‑93 (MQ=255)

                                  ccTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATTAGCGCGGAAGCCGCAGGTATCGCTGTc                                                       >  2:1644441/1‑93 (MQ=255)

                                  ccTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTc                                                       <  2:3176265/93‑1 (MQ=255)

                                     gATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGc                                                    <  1:3320764/93‑1 (MQ=255)

                                       tGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCt                                                  >  2:300281/1‑93 (MQ=255)

                                       tGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCt                                                  >  1:544989/1‑93 (MQ=255)

                                       tGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCt                                                  >  1:1793189/1‑93 (MQ=255)

                                        gCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTg                                                 >  1:1158790/1‑93 (MQ=255)

                                         cTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGa                                                <  1:1682584/93‑1 (MQ=255)

                                           gACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGAtt                                              <  1:841166/93‑1 (MQ=255)

                                           gACGGTCACGATAAATGCCATTTCCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGAtt                                              <  2:926881/93‑1 (MQ=38)

                                           gACGGGGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGAtt                                              <  1:841165/93‑1 (MQ=255)

                                               gTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAat                                          >  2:1050038/1‑93 (MQ=255)

                                               gTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAat                                          >  2:1063122/1‑93 (MQ=255)

                                               gTGACGATAAAGGGCGTTTGCTCAACGGAATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAat                                          <  1:2635782/93‑1 (MQ=38)

                                                tGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAata                                         >  2:2039610/1‑93 (MQ=255)

                                                  aCGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAAtata                                       >  2:1519800/1‑93 (MQ=255)

                                                    gATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGc                                     >  1:4046135/1‑93 (MQ=255)

                                                    gATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGc                                     >  2:693842/1‑93 (MQ=255)

                                                     aTAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGcc                                    <  1:2966739/93‑1 (MQ=255)

                                                     aTAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGcc                                    >  2:4148878/1‑93 (MQ=255)

                                                     aTAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGcc                                    <  1:3938628/93‑1 (MQ=255)

                                                      tAAATGGCATTTGTTCAACGCAAAAAATGGAAATGGCGCGGAAATGAGCGTGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCa                                   <  2:634606/93‑1 (MQ=25)

                                                      tAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCa                                   >  2:3065098/1‑93 (MQ=255)

                                                      tAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGCAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCa                                   >  2:4065656/1‑93 (MQ=255)

                                                       aaaTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCAt                                  >  1:704641/1‑93 (MQ=255)

                                                       aaaTGGCATTAGCTCAACGGCATGAATGGCAATGGCTCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCAt                                  <  2:2775477/93‑1 (MQ=255)

                                                        aaTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATc                                 >  1:2188713/1‑93 (MQ=255)

                                                           ggCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCAcc                              >  1:65177/1‑92 (MQ=255)

                                                           ggCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACg                              >  2:371044/1‑93 (MQ=255)

                                                             cATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcc                            <  1:1248852/93‑1 (MQ=255)

                                                              aTTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcct                           <  2:468880/93‑1 (MQ=255)

                                                              aTTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcct                           >  1:2422136/1‑93 (MQ=255)

                                                               tttGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcctg                          >  1:3783906/1‑93 (MQ=255)

                                                               tttGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcctg                          >  1:421061/1‑93 (MQ=255)

                                                                 tGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcctgcc                        <  1:1536353/93‑1 (MQ=255)

                                                                 tGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcctgcc                        >  1:3970205/1‑93 (MQ=255)

                                                                    tCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCa                     >  2:375573/1‑93 (MQ=255)

                                                                     cAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCAc                    >  2:293075/1‑93 (MQ=255)

                                                                         ggCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGaa                >  1:2791687/1‑93 (MQ=255)

                                                                           cATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGAATGAATATAGCCATCACGCCTGCCTCACCGAATg              >  1:4572859/1‑93 (MQ=255)

                                                                            aTGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGc             >  2:3188988/1‑93 (MQ=255)

                                                                            aTGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGc             >  1:101670/1‑93 (MQ=255)

                                                                            aTGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGc             >  1:538595/1‑93 (MQ=255)

                                                                             tGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCg            >  1:1078793/1‑93 (MQ=255)

                                                                             tGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCg            >  2:3542049/1‑93 (MQ=255)

                                                                              gAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGa           <  2:4264382/93‑1 (MQ=255)

                                                                                aTGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATg         >  1:425810/1‑93 (MQ=255)

                                                                                  ggCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGcc       >  2:4060532/1‑93 (MQ=255)

                                                                                   gCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCa      <  1:1974048/93‑1 (MQ=255)

                                                                                    cAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCAt     >  1:2397197/1‑93 (MQ=255)

                                                                                     aaTGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATg    >  2:3221107/1‑93 (MQ=255)

                                                                                     aaTGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATg    >  2:4687268/1‑93 (MQ=255)

                                                                                     aaTGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATg    >  1:2528095/1‑93 (MQ=255)

                                                                                      aTGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGa   <  2:4476872/93‑1 (MQ=255)

                                                                                       tGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAc  <  1:1618398/93‑1 (MQ=255)

TTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC  >  NC_000913/2774143‑2774322

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