BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	628,816	A→G	V98V (GTA→GTG)	ybdB →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	628,816	0	A	G	99.6%	814.5 / ‑5.7	261	V98V (GTA→GTG)	ybdB	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (125/134); total (126/135)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 271 total aligned reads displayed.

GGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCAC  >  NC_000913/628726‑628905

ggCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGcg                                                                                         <  1:6535560/93‑1 (MQ=255)

ggCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGcg                                                                                         <  2:1924391/93‑1 (MQ=255)

  ccGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGcgcg                                                                                       >  1:3370072/1‑93 (MQ=255)

  ccGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGcgcg                                                                                       >  2:2395245/1‑93 (MQ=255)

    ggATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggc                                                                                     <  1:3845761/93‑1 (MQ=255)

     gATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCgcggcg                                                                                    <  1:2337883/93‑1 (MQ=255)

       tttATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTc                                                                                  <  1:2402543/93‑1 (MQ=255)

       tttATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGCCGTc                                                                                  >  2:2665621/1‑93 (MQ=255)

       tgtATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTc                                                                                  <  1:607849/91‑1 (MQ=255)

         tATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTg                                                                                >  1:2022586/1‑93 (MQ=255)

          atgatgACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGc                                                                               >  2:605949/1‑93 (MQ=255)

           tgatgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGcc                                                                              >  1:911187/1‑93 (MQ=255)

           tgatgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGcc                                                                              <  1:6371378/93‑1 (MQ=255)

           tgatgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGcc                                                                              <  2:6361265/93‑1 (MQ=255)

            gatgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCa                                                                             <  1:528876/93‑1 (MQ=255)

              tgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGc                                                                           >  1:2543880/1‑93 (MQ=255)

              tgaCCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGc                                                                           <  1:3475547/93‑1 (MQ=255)

                 cccGCGGCGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCGCGGTGCCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGc                                                                        <  1:1622119/93‑1 (MQ=255)

                 cccGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGc                                                                        <  1:1123335/93‑1 (MQ=255)

                 cccGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGc                                                                        <  1:428782/93‑1 (MQ=255)

                   cgcgACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTg                                                                      >  1:5641253/1‑93 (MQ=255)

                     cgACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCa                                                                    >  1:6228998/1‑93 (MQ=255)

                     cgACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCa                                                                    >  2:5814506/1‑93 (MQ=255)

                      gACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCAt                                                                   >  1:6779821/1‑93 (MQ=255)

                       aCGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACCCCATCCCCCGGTGTCTGAGGGAAAGGTTCGCGGCGTCTGCCAGCCGCTGCATc                                                                  <  2:1504586/93‑1 (MQ=255)

                       aCGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATc                                                                  <  2:1141253/93‑1 (MQ=255)

                       aCGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATc                                                                  >  2:5436089/1‑93 (MQ=255)

                        cGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCt                                                                 <  1:6809022/93‑1 (MQ=255)

                         ggACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCtt                                                                <  2:3125720/93‑1 (MQ=255)

                            cAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGt                                                             >  1:4697153/1‑93 (MQ=255)

                              gtgtgtGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCg                                                           >  2:5539868/1‑93 (MQ=255)

                              gtgtgtGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCg                                                           >  1:2629294/1‑93 (MQ=255)

                              gtgtgtGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCg                                                           <  1:2812405/93‑1 (MQ=255)

                              gtgtgtGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCg                                                           >  2:6300606/1‑93 (MQ=255)

                                 tgtgGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCa                                                        >  2:4965215/1‑93 (MQ=255)

                                    ggTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCaaaa                                                     <  1:1681420/93‑1 (MQ=255)

                                     gTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAAt                                                    >  2:2807993/1‑93 (MQ=255)

                                     gTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAAt                                                    >  2:3657885/1‑93 (MQ=255)

                                        agCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCag                                                 <  2:1693391/92‑1 (MQ=255)

                                         gCACAGATCTTAATGCAACACCCGATAGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCaga                                                <  2:5618678/93‑1 (MQ=255)

                                         gCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCaga                                                <  1:4574809/93‑1 (MQ=255)

                                          cacaGAACTTAATGCAACACCCCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCagag                                               <  1:1305155/93‑1 (MQ=255)

                                          cacaGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCagag                                               <  1:7852182/93‑1 (MQ=255)

                                            caGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCt                                             >  2:5911451/1‑93 (MQ=255)

                                            caGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCt                                             >  1:1835608/1‑93 (MQ=255)

                                             aGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTg                                            <  2:2146747/93‑1 (MQ=255)

                                              gAACTTAATGCACCACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTgg                                           <  1:501296/93‑1 (MQ=255)

                                              gAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTgg                                           <  2:585848/93‑1 (MQ=255)

                                              gAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTgg                                           <  2:7092354/93‑1 (MQ=255)

                                               aaCTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTggg                                          >  2:5021058/1‑93 (MQ=255)

                                                 cTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaa                                        <  2:7208078/93‑1 (MQ=255)

                                                  ttAATGCAACTCACCATCGCCCTGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaaa                                       <  2:6796853/93‑1 (MQ=255)

                                                  ttAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaaa                                       <  1:4949269/93‑1 (MQ=255)

                                                  ttAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGaaa                                       >  1:696164/1‑93 (MQ=255)

                                                   tAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAt                                      >  1:1389596/1‑93 (MQ=255)

                                                    aaTGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtc                                     >  1:3534598/1‑93 (MQ=255)

                                                     aTGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtcg                                    <  2:3077276/93‑1 (MQ=255)

                                                     aTGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAAtcg                                    <  2:3331142/93‑1 (MQ=255)

                                                         aaCACACCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATcgtcgt                                >  1:7731499/1‑93 (MQ=255)

                                                          acacacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGtt                               >  2:1795962/1‑93 (MQ=255)

                                                          acacacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGtt                               <  1:2156862/93‑1 (MQ=255)

                                                           cacacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGttt                              >  2:1384538/1‑93 (MQ=255)

                                                           cacacCATCGCCCCGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGttt                              >  2:7769507/1‑93 (MQ=255)

                                                             cacCATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGttttg                            >  2:6235121/1‑92 (MQ=255)

                                                               ccATCGCCGGGTGTCGGAGGGAAAGGGGCGGGGCGTCTGCCAGCCGCGGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGa                          <  2:1641030/93‑1 (MQ=255)

                                                               ccATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCgt                          >  1:4499045/1‑92 (MQ=255)

                                                                cATCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGAt                         <  1:2987720/93‑1 (MQ=255)

                                                                cATCGCCAGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGAt                         >  2:4594001/1‑93 (MQ=255)

                                                                 aTCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATg                        >  2:6454221/1‑93 (MQ=255)

                                                                  tCGCCCGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGa                       >  1:3128542/1‑93 (MQ=255)

                                                                     cccGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACa                    <  2:3761387/93‑1 (MQ=255)

                                                                      ccGGTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAg                   <  2:5319299/93‑1 (MQ=255)

                                                                       cGGTGTCTGAGGGAAAGATGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgg                  <  1:7025135/93‑1 (MQ=255)

                                                                        ggTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAggg                 >  1:7230659/1‑93 (MQ=255)

                                                                        ggTGTCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATCAACAggg                 <  2:7520258/93‑1 (MQ=255)

                                                                         gtgtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgggg                >  2:3532004/1‑93 (MQ=255)

                                                                         gtgtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgggg                <  1:3747497/93‑1 (MQ=255)

                                                                         gtgtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgggg                <  1:440089/93‑1 (MQ=255)

                                                                         gtgtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCGTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAgggg                <  1:7738838/93‑1 (MQ=255)

                                                                          tgtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggc               >  2:7885981/1‑93 (MQ=255)

                                                                           gtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcg              >  1:5471309/1‑93 (MQ=255)

                                                                           gtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcg              >  1:2312479/1‑93 (MQ=255)

                                                                           gtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcg              >  1:7647930/1‑93 (MQ=255)

                                                                           gtCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcg              >  1:6640610/1‑93 (MQ=255)

                                                                            tCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcgg             <  1:1018712/93‑1 (MQ=255)

                                                                            tCTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcgg             >  2:1308329/1‑93 (MQ=255)

                                                                            gcTGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcgg             <  2:7401234/92‑1 (MQ=255)

                                                                             cggAGGGAAAGGGGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGggcggc            <  2:5666148/91‑1 (MQ=255)

                                                                              tGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGgcggcg           >  2:2006744/1‑93 (MQ=255)

                                                                              tGAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGgcggcg           >  1:6552742/1‑93 (MQ=255)

                                                                               gAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGt          <  1:3245114/93‑1 (MQ=255)

                                                                               gAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGt          >  1:4153180/1‑93 (MQ=255)

                                                                               gAGGGAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGt          <  2:2603443/93‑1 (MQ=255)

                                                                                 gggAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtg        >  1:7625379/1‑93 (MQ=255)

                                                                                 gggAAAGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtg        >  1:5899799/1‑93 (MQ=255)

                                                                                     aaGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTtgctgc    >  2:4816834/1‑93 (MQ=255)

                                                                                      aGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGCGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCa   >  2:3422312/1‑93 (MQ=255)

                                                                                      aGGTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCa   <  1:1395893/93‑1 (MQ=255)

                                                                                       ggTGCGCGTCGTGTGCCCGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCAc  <  1:5557054/93‑1 (MQ=255)

                                                                                       ggTGCGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCAc  >  2:5189874/1‑93 (MQ=255)

GGCCGGATTTATGATGACCCGCGACGGACAGTGTGTGGTAGGCACAGAACTTAATGCAACACACCATCGCCCGGTGTCTGAGGGAAAGGTACGCGGCGTCTGCCAGCCGCTGCATCTTGGTCGGCAAAATCAGAGCTGGGAAATCGTCGTTTTCGATGAACAGGGGCGGCGTTGCTGCAC  >  NC_000913/628726‑628905

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