BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	527,625	T→C	G151G (GGT→GGC)	b0499 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	527,625	0	T	C	88.0%	692.4 / 76.5	276	G151G (GGT→GGC)	b0499	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (15/18); new base C (118/124); total (134/142)
Fisher's exact test for biased strand distribution p-value = 8.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

Only 100 of 279 total aligned reads displayed.

TGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATGTACGGCAGAAAAGTAAAG  >  NC_000913/527535‑527715

tgaagatatgaagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGgcgt                                                                                          <  2:5496582/93‑1 (MQ=255)

 gaagatatgaagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGtt                                                                                         >  1:808852/1‑93 (MQ=255)

  aagatatgaagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTa                                                                                        <  2:3003255/93‑1 (MQ=255)

   agatatgaagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAg                                                                                       <  1:7407149/93‑1 (MQ=255)

      tatgaagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCg                                                                                    >  2:5054545/1‑93 (MQ=255)

        tgaagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGAt                                                                                  >  1:4604089/1‑93 (MQ=255)

         gaagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATc                                                                                 <  2:3377731/93‑1 (MQ=255)

         gaagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATc                                                                                 <  1:7140583/93‑1 (MQ=255)

         gaagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATc                                                                                 >  1:6879291/1‑93 (MQ=255)

         gaagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATc                                                                                 <  1:3455841/93‑1 (MQ=255)

         gaagaGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATc                                                                                 >  2:7345687/1‑93 (MQ=255)

           agagATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCgg                                                                               <  2:2360414/93‑1 (MQ=255)

           agagATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCgg                                                                               >  1:3874385/1‑93 (MQ=255)

            gagaTTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGc                                                                              <  2:3886584/93‑1 (MQ=255)

            gagaTTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGc                                                                              <  1:2272987/93‑1 (MQ=255)

              gaTTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGa                                                                            <  1:1155993/93‑1 (MQ=255)

                ttAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGaaa                                                                          >  1:282453/1‑93 (MQ=255)

                ttAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGaaa                                                                          <  1:3223841/93‑1 (MQ=255)

                ttAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGaaa                                                                          <  1:2091026/93‑1 (MQ=255)

                ttAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGaaa                                                                          <  1:4643530/93‑1 (MQ=255)

                ttAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATCGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGaaa                                                                          <  1:5817467/93‑1 (MQ=255)

                 tAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAg                                                                         <  1:4095348/93‑1 (MQ=255)

                 tAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAg                                                                         <  1:1964561/93‑1 (MQ=255)

                 tAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAg                                                                         <  2:6259231/93‑1 (MQ=255)

                   aaTTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAggg                                                                       >  2:3636754/1‑93 (MQ=255)

                   aaTTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAggg                                                                       >  1:65466/1‑93 (MQ=255)

                    aTTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGc                                                                      <  1:2708093/93‑1 (MQ=255)

                    aTTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGc                                                                      <  2:6795456/93‑1 (MQ=255)

                    aTTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGc                                                                      <  1:5565593/93‑1 (MQ=255)

                    aTTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGc                                                                      <  2:2651540/93‑1 (MQ=255)

                       tAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTgg                                                                   <  1:6533254/93‑1 (MQ=255)

                        aGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTAGCGTTAGCCGATCGGCGAAAGGGCTggg                                                                  <  1:1531711/93‑1 (MQ=255)

                         ggTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGt                                                                 <  2:5640342/93‑1 (MQ=255)

                         ggTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGt                                                                 <  2:4455853/93‑1 (MQ=255)

                          gTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGtt                                                                <  1:4007941/93‑1 (MQ=255)

                          gTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGtt                                                                >  2:1818213/1‑93 (MQ=255)

                            gggACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTAt                                                              <  1:2411023/93‑1 (MQ=255)

                              gACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGa                                                            >  1:1878/1‑93 (MQ=255)

                               aCTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGaa                                                           >  2:663237/1‑93 (MQ=255)

                                cTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGaaa                                                          >  1:3725690/1‑93 (MQ=255)

                                 tGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGaaaa                                                         >  2:6126061/1‑93 (MQ=255)

                                 tGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGaaaa                                                         >  2:2516772/1‑93 (MQ=255)

                                 tGATCAATTTTTCCATTGTATGGCATTTTGTCGAGCGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGaaaa                                                         >  2:484783/1‑93 (MQ=255)

                                  gATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGaaaaa                                                        >  2:5256823/1‑93 (MQ=255)

                                  gATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGaaaaa                                                        >  2:3557874/1‑93 (MQ=255)

                                   aTCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAag                                                       >  2:2202529/1‑93 (MQ=255)

                                    tCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAaga                                                      <  1:3806060/93‑1 (MQ=255)

                                    tCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAaga                                                      >  2:654506/1‑93 (MQ=255)

                                      aaTTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagaga                                                    >  1:1010285/1‑93 (MQ=255)

                                         ttttCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagagatta                                                 <  1:300139/93‑1 (MQ=255)

                                         ttttCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagagatta                                                 <  1:1664055/93‑1 (MQ=255)

                                         ttttCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagagatta                                                 <  2:4280935/93‑1 (MQ=255)

                                          tttCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagagattag                                                <  2:5679892/93‑1 (MQ=255)

                                          tttCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagagattag                                                >  1:4750890/1‑93 (MQ=255)

                                           ttCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagagattaga                                               <  1:5092268/93‑1 (MQ=255)

                                             ccATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagagattagaga                                             <  1:1729231/93‑1 (MQ=255)

                                             ccATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAagagactagaga                                             <  2:6987810/93‑1 (MQ=255)

                                                ttGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAgagattagagatta                                          <  1:2240540/93‑1 (MQ=255)

                                                ttGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAgagattagagatta                                          >  2:3285732/1‑93 (MQ=255)

                                                ttGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAgagattagagatta                                          >  2:1440866/1‑93 (MQ=255)

                                                 tGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTAc                                         <  2:846504/93‑1 (MQ=255)

                                                 tGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTAc                                         >  1:60341/1‑93 (MQ=255)

                                                  gTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACg                                        <  1:1467266/93‑1 (MQ=255)

                                                  gTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACg                                        <  1:539547/93‑1 (MQ=255)

                                                  gTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACg                                        <  1:5749119/93‑1 (MQ=255)

                                                  gTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACg                                        >  2:4591031/1‑93 (MQ=255)

                                                  gTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACg                                        <  2:3732205/93‑1 (MQ=255)

                                                   tATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACgg                                       <  2:1647379/93‑1 (MQ=255)

                                                   tATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACgg                                       <  2:4714457/93‑1 (MQ=255)

                                                    aTGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACggg                                      >  2:5891358/1‑93 (MQ=255)

                                                     tGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGt                                     >  2:636926/1‑93 (MQ=255)

                                                       gCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTa                                   >  1:4252739/1‑93 (MQ=255)

                                                        cATTTTGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTaa                                  >  2:2479632/1‑93 (MQ=255)

                                                          ttttGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAAt                                <  2:3197114/93‑1 (MQ=255)

                                                          ttttGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAAt                                >  1:4368288/1‑93 (MQ=255)

                                                            ttGTCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCt                              >  2:15572/1‑93 (MQ=255)

                                                               tCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGtt                           <  1:4901320/93‑1 (MQ=255)

                                                               tCGAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGtt                           >  1:4148322/1‑93 (MQ=255)

                                                                ggggtgtCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTc                          <  1:6099047/90‑1 (MQ=255)

                                                                cGAGTGTCTAAATTAAAGGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTc                          <  1:6250270/93‑1 (MQ=255)

                                                                 gAGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCg                         >  1:3733738/1‑93 (MQ=255)

                                                                  aGTGTCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCgg                        >  1:7531048/1‑93 (MQ=255)

                                                                    tgtCTAAATTAAATGACGCTGGCGTTGGCCGATCGGCGAAAGGCCTGGGTTATGAAAAAGAGATTAGCGATTACGGGTTAAATCTGTTCGGTa                      >  1:871928/1‑93 (MQ=255)

                                                                    tgtATAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTa                      <  2:3931031/93‑1 (MQ=255)

                                                                     gtCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTAt                     >  2:2918116/1‑93 (MQ=255)

                                                                      tCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATg                    <  2:1332536/93‑1 (MQ=255)

                                                                      tCTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATg                    >  2:1119541/1‑93 (MQ=255)

                                                                       cTAAATTAAATGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATGt                   >  2:6203340/1‑93 (MQ=255)

                                                                        taaattaaatGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATGTa                  >  1:5974004/1‑93 (MQ=255)

                                                                        taaattaaatGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATGTa                  <  1:7237321/93‑1 (MQ=255)

                                                                          aattaaatGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATGTACg                <  1:5489830/93‑1 (MQ=255)

                                                                           attaaatGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATGTACgg               >  2:1289611/1‑93 (MQ=255)

                                                                            ttaaatGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATGTACgac              >  2:1765754/1‑91 (MQ=255)

                                                                              aaatGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATGTACGGCAg            >  2:7521436/1‑93 (MQ=255)

                                                                                atGACGCTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATGTACGGCAGaa          >  2:7022173/1‑93 (MQ=255)

                                                                                      cTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGATATGTACGGCAGAAAAGTaa    <  2:5390324/93‑1 (MQ=255)

                                                                                      cTGGGGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATGTACGGCAGAAAAGTaa    >  1:3035786/1‑93 (MQ=255)

                                                                                      cTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATGTACGGCAGAAAAGTaa    >  1:5648628/1‑93 (MQ=255)

                                                                                      cTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATGTACGGCAGAAAAGTaa    <  2:5118402/93‑1 (MQ=255)

                                                                                      cTGGCGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATGTACGGCAGAAAAGTaa    <  2:4950833/93‑1 (MQ=255)

TGAAGATATGAAGAGATTAAATTTAGGTGGGACTGATCAATTTTTCCATTGTATGGCATTTTGTCGAGTGTCTAAATTAAATGACGCTGGTGTTAGCCGATCGGCGAAAGGGCTGGGTTATGAAAAAGAGATTAGAGATTACGGGTTAAATCTGTTCGGTATGTACGGCAGAAAAGTAAAG  >  NC_000913/527535‑527715

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