BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,577,071	A→G	*114Q (TAA→CAA)	yjiW ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,577,071	0	A	G	99.5%	737.6 / ‑5.2	209	*114Q (TAA→CAA)	yjiW	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (1/0); new base G (94/114); total (95/114)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 225 total aligned reads displayed.

AAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCG  >  NC_000913/4576981‑4577160

aaaaTCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGcg                                                                                         <  2:1670625/93‑1 (MQ=255)

aaaaTCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAAACGGGAGTAAACGCTGTTGcg                                                                                         >  1:1087940/1‑93 (MQ=255)

 aaaTCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGGAAACGCTGTTgcgc                                                                                        >  1:5393023/1‑93 (MQ=255)

   aTCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGa                                                                                      <  1:3421068/93‑1 (MQ=255)

    acAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGAc                                                                                     <  2:4866083/92‑1 (MQ=255)

      aGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGGGCGACtt                                                                                   >  1:4675376/1‑93 (MQ=255)

      aGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACtt                                                                                   <  1:5140967/93‑1 (MQ=255)

        aCTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTTTTGCGCGACTTTc                                                                                 <  1:4843478/93‑1 (MQ=255)

        aCTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACtttt                                                                                 >  2:5374535/1‑92 (MQ=255)

        aCTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTc                                                                                 <  1:3084810/93‑1 (MQ=255)

        aCTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGCGGTAAAAGCTGTTGCGCGACATTc                                                                                 >  2:1815115/1‑93 (MQ=255)

         cTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCt                                                                                >  1:1956165/1‑93 (MQ=255)

         cTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCt                                                                                >  2:4156267/1‑93 (MQ=255)

          tGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTg                                                                               <  1:209597/93‑1 (MQ=255)

           ggTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGt                                                                              >  2:5238530/1‑93 (MQ=255)

           ggTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGGGGTTAACCGGGAGGAAAAGCTGTTGCGCGACTTTCTGt                                                                              >  2:323894/1‑93 (MQ=255)

            gTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGAGAGTAAACGCTGTTGCGCGACTTTCTGtt                                                                             <  2:4273442/93‑1 (MQ=255)

               aaTATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGGTCCCGGTGGTTAACCGGGAGGAAACGCCGTTGCGCGACTTTCTGTTTAc                                                                          >  2:3521822/1‑93 (MQ=255)

                atatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTAcc                                                                         >  1:4058235/1‑93 (MQ=255)

                atatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTAcc                                                                         <  1:5268423/93‑1 (MQ=255)

                atatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTAcc                                                                         >  2:5711354/1‑93 (MQ=255)

                 tatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTGACCGGGAGTAAAAGCCGTTTCGCGACTTTTTGTTTAACg                                                                        >  1:3514737/1‑93 (MQ=255)

                 tatTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCg                                                                        >  2:5849284/1‑93 (MQ=255)

                  atTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTAACgg                                                                       >  2:3871074/1‑93 (MQ=255)

                  atTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGGGGGTAACCGGGAGGAAACGCTGGTGCGCGACGTTCTGTTTACCgg                                                                       >  2:237075/1‑93 (MQ=255)

                   tttACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAAAGCTGTTGCGCGACTTTCTGTTTACCGGc                                                                      >  2:3018250/1‑93 (MQ=255)

                     tACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGGGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTTTGTTTACCGGCaa                                                                    >  1:5002514/1‑93 (MQ=255)

                       cTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATc                                                                  <  2:885080/93‑1 (MQ=255)

                       cTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATc                                                                  <  1:6138310/93‑1 (MQ=255)

                        tGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAAttc                                                                 >  1:3376789/1‑91 (MQ=255)

                        tGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCa                                                                 <  2:2914341/93‑1 (MQ=255)

                        tGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCa                                                                 <  1:5727249/93‑1 (MQ=255)

                            ggggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTcc                                                             <  1:2939248/93‑1 (MQ=255)

                             gggAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGCCTTTCTGTTTACCGGCAATCACTCCa                                                            >  2:902807/1‑93 (MQ=255)

                               gAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAAt                                                          <  1:2371269/93‑1 (MQ=255)

                                aaaGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATa                                                         <  1:3599628/93‑1 (MQ=255)

                                 aaGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGCCTTTCTGTTTACCGGCAATCACTCCAATaa                                                        >  1:298483/1‑93 (MQ=255)

                                 aaGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTTTGTTTACCGGCAATCAATCCAATaa                                                        >  1:1468030/1‑93 (MQ=255)

                                 aaGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATaa                                                        >  2:2465032/1‑93 (MQ=255)

                                 aaGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATaa                                                        <  1:4901117/93‑1 (MQ=255)

                                  aGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATaaa                                                       <  2:2258821/93‑1 (MQ=255)

                                  aGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATaaa                                                       <  1:2564316/93‑1 (MQ=255)

                                     ttCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGc                                                    <  1:1242288/93‑1 (MQ=255)

                                      tCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCCATAAACGcc                                                   >  1:5840047/1‑93 (MQ=255)

                                      tCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGcc                                                   >  2:3586747/1‑93 (MQ=255)

                                       cTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCt                                                  >  1:2176767/1‑93 (MQ=255)

                                        tATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTg                                                 >  1:1881239/1‑93 (MQ=255)

                                         aTTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGc                                                <  2:4363637/93‑1 (MQ=255)

                                         aTTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGc                                                <  2:2597133/93‑1 (MQ=255)

                                         aTTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGc                                                >  1:4243976/1‑93 (MQ=255)

                                         aTTGAGTCAGTGGAAGGCTCCCGGGGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGc                                                >  1:1516558/1‑93 (MQ=255)

                                          ttGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCa                                               >  2:402240/1‑93 (MQ=255)

                                          ttGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCa                                               >  1:942427/1‑93 (MQ=255)

                                           tGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCAc                                              <  1:4951362/93‑1 (MQ=255)

                                           tGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCAc                                              <  2:5077812/93‑1 (MQ=255)

                                            gAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCAcc                                             <  1:1430193/93‑1 (MQ=255)

                                                 aGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAACCGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCtt                                        <  2:269865/93‑1 (MQ=255)

                                                 aGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTGCCGGCAATCACTCCAATAAACGCCTGCACCTGCtt                                        <  1:1615996/93‑1 (MQ=255)

                                                  gTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCCCCTGCttt                                       >  2:5363572/1‑93 (MQ=255)

                                                   tGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCtttt                                      >  1:3870787/1‑93 (MQ=255)

                                                      aaGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGtt                                   <  1:4371810/93‑1 (MQ=255)

                                                        ggCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTa                                 <  2:4988357/93‑1 (MQ=255)

                                                        ggCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTa                                 >  1:3619010/1‑93 (MQ=255)

                                                        ggCTCCCGGTGGTTAACCGGGAGTAAAAGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTa                                 >  2:6042114/1‑93 (MQ=255)

                                                          cTCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcg                               <  2:4687551/93‑1 (MQ=255)

                                                           tCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcgc                              <  1:5316389/93‑1 (MQ=255)

                                                           tCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcgc                              >  2:2945488/1‑93 (MQ=255)

                                                           tCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcgc                              >  1:1746164/1‑93 (MQ=255)

                                                           tCCCGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcgc                              <  2:771252/93‑1 (MQ=255)

                                                            cccGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcgcg                             >  2:3817112/1‑93 (MQ=255)

                                                            cccGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTAcgcg                             >  2:4507367/1‑93 (MQ=255)

                                                             ccGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACgcgc                            >  2:5596274/1‑93 (MQ=255)

                                                              cGGTGGTTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGcc                           <  2:2794024/93‑1 (MQ=255)

                                                               ggtggtTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCg                          >  1:4030039/1‑93 (MQ=255)

                                                                gtggtTAAGCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGa                         <  1:2791732/93‑1 (MQ=255)

                                                                 tggtTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGAc                        <  2:3154001/93‑1 (MQ=255)

                                                                  ggtTAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACa                       >  1:1834482/1‑93 (MQ=255)

                                                                    ttAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGt                     >  1:5546365/1‑93 (MQ=255)

                                                                    ttAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGt                     <  1:4505726/93‑1 (MQ=255)

                                                                     tAACCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGtt                    <  1:75708/93‑1 (MQ=255)

                                                                      aaCCGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGttt                   <  2:4403941/93‑1 (MQ=255)

                                                                         cGGGAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGca                <  2:193093/93‑1 (MQ=255)

                                                                          gggAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGcac               >  2:1895305/1‑93 (MQ=255)

                                                                           ggAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGcaca              <  1:479262/93‑1 (MQ=255)

                                                                           ggAGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGcaca              <  2:1471362/93‑1 (MQ=255)

                                                                             aGTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACAcc            <  1:4861099/93‑1 (MQ=255)

                                                                             aGTAAACCCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAATCGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCATAcc            <  2:1333969/93‑1 (MQ=255)

                                                                              gTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCt           <  2:5262248/93‑1 (MQ=255)

                                                                              gTAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCt           <  1:4919243/93‑1 (MQ=255)

                                                                               tAAACGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTg          <  1:5959537/93‑1 (MQ=255)

                                                                                  aCGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCg       >  1:3835744/1‑93 (MQ=255)

                                                                                  aCGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCg       >  2:1752071/1‑93 (MQ=255)

                                                                                   cGCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGt      >  2:272219/1‑93 (MQ=255)

                                                                                    gCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTa     <  1:2156440/93‑1 (MQ=255)

                                                                                    gCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTa     <  2:3662910/93‑1 (MQ=255)

                                                                                    gCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTa     <  1:187432/93‑1 (MQ=255)

                                                                                    gCTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTa     >  2:118375/1‑93 (MQ=255)

                                                                                     cTGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAg    >  1:1858099/1‑93 (MQ=255)

                                                                                      tGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAgg   >  2:3076167/1‑92 (MQ=255)

                                                                                      tGTTGCGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGc   >  2:5806786/1‑93 (MQ=255)

AAAATCAGACTGGTAAATATTTACTGAGGGGGAAAGTTTCTATTGAGTCAGTGGAAGGCTCCCGGTGGTTAACCGGGAGTAAACGCTGTTACGCGACTTTCTGTTTACCGGCAATCACTCCAATAAACGCCTGCACCTGCTTTTGTTTACGCGCCGACAGTTTGCACACCTGGCGTAGCG  >  NC_000913/4576981‑4577160

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