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	99.6%	752.9 / NA	242	G98G (GGT→GGC)	yfjX	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (116/125); minor base A (1/0); total (117/125)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 270 total aligned reads displayed.

TTACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC  >  NC_000913/2774143‑2774322

ttACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGGGACGATAAATGGCATTTGCTCAACAGCATGAATGGCAATGgcgc                                                                                         >  2:6801280/1‑93 (MQ=255)

 tACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGcgcg                                                                                        <  1:2465223/93‑1 (MQ=255)

     cGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGaaa                                                                                    <  2:4120541/93‑1 (MQ=255)

      gCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAAt                                                                                   >  1:1221810/1‑93 (MQ=255)

        tCAGCAATGGCGGAGCTTTTATGGCTCCTGATGCTGACAGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGa                                                                                 <  2:7068352/93‑1 (MQ=255)

        tCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGa                                                                                 >  1:5406199/1‑93 (MQ=255)

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

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

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

        tCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGa                                                                                 >  1:3192530/1‑93 (MQ=255)

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

         cAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAg                                                                                >  1:5924864/1‑93 (MQ=255)

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

           gCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATTGCAATGGCGCGGAAACGAgcg                                                                              >  1:4157798/1‑93 (MQ=255)

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

                  cGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGc                                                                       <  1:4932756/93‑1 (MQ=255)

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

                   ggAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGcc                                                                      >  2:6580103/1‑93 (MQ=255)

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

                     aGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGc                                                                    <  1:1793460/93‑1 (MQ=255)

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

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

                      gCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCa                                                                   <  1:3547715/93‑1 (MQ=255)

                        gTTTATGGCTCCTGCTGCTGACGGTGACGATAAATGGCATTTGCTCAACGGGATGAACGGAAATGACGCGGTAATGAGGGCAGAAGCCGGAgg                                                                 >  2:3685516/1‑93 (MQ=255)

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

                         tttATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCAGCCGGt                                                                >  1:4428730/1‑93 (MQ=255)

                          ttATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAggga                                                               >  2:526646/1‑91 (MQ=255)

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

                          ttATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTa                                                               >  2:760665/1‑93 (MQ=255)

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

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

                             tGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCg                                                            >  1:2662762/1‑93 (MQ=255)

                              ggCTCCTGATGTGGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGc                                                           <  1:1894191/93‑1 (MQ=255)

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

                              ggCTCCTGATGCTGACGGTGACGATAAATGGCATCTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGc                                                           <  1:4834991/93‑1 (MQ=255)

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

                                    tGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTg                                                     >  1:6361208/1‑93 (MQ=255)

                                      aTGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGcc                                                   <  2:3982552/93‑1 (MQ=255)

                                      aTGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGcc                                                   <  2:2294281/93‑1 (MQ=255)

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

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

                                         cTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGa                                                >  2:4856173/1‑93 (MQ=255)

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

                                          tGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGAt                                               >  1:3235340/1‑93 (MQ=255)

                                          tGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGAt                                               <  2:4706276/93‑1 (MQ=255)

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

                                              ggTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGaa                                           <  2:2108184/93‑1 (MQ=255)

                                                tGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAata                                         <  2:6276704/93‑1 (MQ=255)

                                                tGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAata                                         <  2:781143/93‑1 (MQ=255)

                                                 gACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAatat                                        <  2:6153716/93‑1 (MQ=255)

                                                 gACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAatat                                        >  1:663330/1‑93 (MQ=255)

                                                   cGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAg                                      >  1:6985285/1‑93 (MQ=255)

                                                     aTAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAgac                                    >  2:3445639/1‑91 (MQ=255)

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

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

                                                       aaaTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGGATCGCTGGCTGCCTGATTGAATATAGCCac                                  >  1:2347188/1‑92 (MQ=255)

                                                        aaTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATc                                 <  2:7212279/93‑1 (MQ=255)

                                                        aaTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATc                                 <  1:7782294/93‑1 (MQ=255)

                                                         aTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCa                                <  1:5808454/93‑1 (MQ=255)

                                                         aTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGAAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCa                                <  2:1832120/93‑1 (MQ=255)

                                                          tGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCAc                               <  1:5085472/93‑1 (MQ=255)

                                                          tGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGAGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATccc                               >  2:3355697/1‑91 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                                                                     cAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCAc                    <  1:4225669/93‑1 (MQ=255)

                                                                      aaCGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCAcc                   <  1:7232789/93‑1 (MQ=255)

                                                                       aCGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCg                  >  1:7551261/1‑93 (MQ=255)

                                                                       aCGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCg                  >  2:5566441/1‑93 (MQ=255)

                                                                       aCGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCg                  >  2:7441235/1‑93 (MQ=255)

                                                                         ggCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGaa                <  2:2530220/93‑1 (MQ=255)

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

                                                                         ggCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGaa                <  2:1411246/93‑1 (MQ=255)

                                                                           cATTAAGGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATg              <  1:1851950/93‑1 (MQ=255)

                                                                            aTGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGc             <  2:7369084/93‑1 (MQ=255)

                                                                            aTGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGc             <  2:1188037/93‑1 (MQ=255)

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

                                                                              gAATGGCAAGGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGa           <  1:6837591/93‑1 (MQ=255)

                                                                               aatggcaatggcGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGAt          >  1:2587092/1‑93 (MQ=255)

                                                                               aatggcaatggcGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGAt          <  2:3732106/93‑1 (MQ=255)

                                                                                 tGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGc        >  1:4306244/1‑93 (MQ=255)

                                                                                 tGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGc        <  2:3265770/93‑1 (MQ=255)

                                                                                 tGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGc        >  2:4419504/1‑93 (MQ=255)

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

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

                                                                                  ggCAATGGCGCGGAAAGGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGcc       <  1:1469223/93‑1 (MQ=255)

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

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

                                                                                    cAATGGCGCGTAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCATGATTGAATATAGCCATCCCGCCTGCCTCTCCGAATGCGATGCCAt     <  2:1957385/93‑1 (MQ=255)

                                                                                    cAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCAt     <  1:5485056/93‑1 (MQ=255)

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

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

                                                                                       tGGCGCGGAAATGAGCGTGGAAGTCGCGGGTCTCGCTGTCTGCCTGATTGAATATAGGCCTCAGCCCTGCCTCACCGAATGCGATGCCATgca  >  2:4735906/1‑91 (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` `< 29 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`