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.5%	626.8 / ‑4.5	199	G98G (GGT→GGC)	yfjX	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (0/0); major base C (98/98); minor base A (2/1); total (100/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 218 total aligned reads displayed.

TACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC  >  NC_000913/2774144‑2774322

tACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGcgcg                                                                                        >  1:6145950/1‑93 (MQ=255)

tACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGcgcg                                                                                        >  2:1714432/1‑93 (MQ=255)

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

   acGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGaa                                                                                     >  2:5686126/1‑93 (MQ=255)

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

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

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

        cAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAg                                                                                <  2:4361190/93‑1 (MQ=255)

         aGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAgg                                                                               >  2:946778/1‑92 (MQ=255)

          gCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGGAATGGCGCGGGAATCAgcg                                                                              >  2:4716110/1‑93 (MQ=255)

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

          gCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAgcg                                                                              <  2:383890/93‑1 (MQ=255)

           cAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAgcgc                                                                             <  2:5329217/93‑1 (MQ=255)

            aaTGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGAAGATAAATGGCATTTGCTCAACGGCATGAATGGGAATGGAGCGGAAATGAGcgcg                                                                            >  2:406635/1‑93 (MQ=255)

             aTGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGAGCgg                                                                           >  1:1406972/1‑93 (MQ=255)

              tGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGa                                                                          <  1:3770995/93‑1 (MQ=255)

               ggCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGaa                                                                         <  2:3559598/93‑1 (MQ=255)

               ggCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGaa                                                                         <  2:2615942/93‑1 (MQ=255)

                gCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGaac                                                                        >  1:4152348/1‑92 (MQ=255)

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

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

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

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

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

                       gTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAgg                                                                 >  1:3042954/1‑93 (MQ=255)

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

                       gTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCAATTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAgg                                                                 >  2:3087422/1‑93 (MQ=255)

                       gTTGATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCCCGGAAGCCGCAgg                                                                 <  2:2823820/93‑1 (MQ=255)

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

                        tttATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGt                                                                <  2:2477955/93‑1 (MQ=255)

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

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

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

                           aTGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATc                                                             <  1:131703/93‑1 (MQ=255)

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

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

                                tgcTGATGCTGACGGTGACGATAAGTGGCTTTTGCGCAACGGCATGAATGGCAGTGGAGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGt                                                        <  2:1188223/91‑1 (MQ=255)

                                  tggATGCTGACGGTGACGATAAATGCCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCt                                                      <  2:2767599/91‑1 (MQ=255)

                                     aTGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGcc                                                   >  2:5115238/1‑93 (MQ=255)

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

                                      tGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCt                                                  <  2:111949/93‑1 (MQ=255)

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

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

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

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

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

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

                                          gACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCGGCCTGAtt                                              <  2:6152145/93‑1 (MQ=255)

                                           aCGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTg                                             >  1:2433669/1‑93 (MQ=255)

                                            cGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGa                                            >  2:6286679/1‑93 (MQ=255)

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

                                              gTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAat                                          <  1:4225644/93‑1 (MQ=255)

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

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

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

                                                 aCGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAAtata                                       <  2:4271200/93‑1 (MQ=255)

                                                 aCGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAAtata                                       <  2:1280565/93‑1 (MQ=255)

                                                   gATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGc                                     <  1:3726565/93‑1 (MQ=255)

                                                     tAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCa                                   <  1:113373/93‑1 (MQ=255)

                                                     tAAATGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGAAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCa                                   <  1:4757905/93‑1 (MQ=255)

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

                                                      aaaTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCAt                                  <  2:786769/93‑1 (MQ=255)

                                                      aaaTGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGCATCGCTGTCTGCCTGATTGAATATAGccct                                  >  2:335621/1‑91 (MQ=255)

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

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

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

                                                         tGGCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCAc                               >  1:5018068/1‑93 (MQ=255)

                                                           gCATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgc                             >  1:4476012/1‑93 (MQ=255)

                                                            cATTTGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcc                            >  1:2004495/1‑93 (MQ=255)

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

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

                                                               ttGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcctgc                         <  2:3777055/93‑1 (MQ=255)

                                                               ttGCTCAACGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACgcctgc                         <  1:5218169/93‑1 (MQ=255)

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

                                                                     aaCGGCATGAATGGCAATGGAGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCAcc                   <  2:187498/93‑1 (MQ=255)

                                                                      aCGGCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCg                  <  1:3838788/93‑1 (MQ=255)

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

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

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

                                                                         gCATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAAt               <  1:9862/93‑1 (MQ=255)

                                                                          cATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATg              >  1:3677500/1‑93 (MQ=255)

                                                                          cATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATg              >  1:1617122/1‑93 (MQ=255)

                                                                          cATGAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATg              <  2:2112018/93‑1 (MQ=255)

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

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

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

                                                                             gAATGGCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGa           >  1:4341776/1‑93 (MQ=255)

                                                                             gAATGGAAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGa           <  1:4941579/93‑1 (MQ=255)

                                                                             gAAGGGCAATGCCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGa           <  1:5751371/93‑1 (MQ=255)

                                                                             gAAGGGAAATGGCGCGGAAATGGGCGCGGAGGCGGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGa           <  1:5369595/93‑1 (MQ=37)

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

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

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

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

                                                                                 ggCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGcc       <  2:2533183/93‑1 (MQ=255)

                                                                                  gCAATGGCGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCa      >  2:5728512/1‑93 (MQ=255)

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

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

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

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

TACACGCTCAGCAATGGCGGAGCGTTTATGGCTCCTGATGCTGACGGTGACGATAAATGGCATTTGCTCAACGGCATGAATGGCAATGGTGCGGAAATGAGCGCGGAAGCCGCAGGTATCGCTGTCTGCCTGATTGAATATAGCCATCACGCCTGCCTCACCGAATGCGATGCCATGAC  >  NC_000913/2774144‑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: `‑`