BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,609,699	A→G	V68V (GTA→GTG)	b3474 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,609,699	0	A	G	99.4%	966.2 / ‑3.9	309	V68V (GTA→GTG)	b3474	hypothetical protein
Reads supporting (aligned to +/- strand): ref base A (2/0); new base G (164/143); total (166/143)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

Only 100 of 355 total aligned reads displayed.

ACTGGCGCTGTTTTTTGCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTAGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTATGCTGCCGAAGTTTAA  >  NC_000913/3609609‑3609788

tcTGGCGCTGTTTTTTGCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGgt                                                                                         <  2:1630335/92‑1 (MQ=255)

 cTGGCGCTTTTTTTTGCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTgtggtg                                                                                        <  1:2613054/93‑1 (MQ=255)

 cTGGCGCTGTTTTTTGCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTgtggtg                                                                                        <  1:2302543/93‑1 (MQ=255)

 cTGGCGCTGTTTTTTGCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGACCCGTTGCCATTACGATTgtggtg                                                                                        >  1:7078469/1‑93 (MQ=255)

   ggCGCTGTTTTTTGCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGtggtggg                                                                                      >  1:7761055/1‑92 (MQ=255)

   ggCGCTGTTTTTTGCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTggtggt                                                                                      <  2:1076420/93‑1 (MQ=255)

     cgcTGTTTTTTGCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGtggtggtgg                                                                                    <  2:4869355/93‑1 (MQ=255)

       cTGTTTTTTGCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTgtggtggtggtg                                                                                  >  1:7387337/1‑93 (MQ=255)

        tGTTTTTTGCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGa                                                                                 >  1:4741112/1‑93 (MQ=255)

        tGTTTTTTGCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGa                                                                                 <  1:1943722/93‑1 (MQ=255)

        tGTTTTTTGCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGa                                                                                 >  1:141755/1‑93 (MQ=255)

         gTTTTTTGCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGAt                                                                                >  1:7227450/1‑93 (MQ=255)

         gTTTTTTGCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGGGGTGCTGAt                                                                                >  1:2060185/1‑93 (MQ=255)

                gCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATc                                                                         <  1:672501/93‑1 (MQ=255)

                   aTCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTTCGATTGTGGTGGTGGTGATGCTGATCGcc                                                                      >  2:6638635/1‑92 (MQ=255)

                   aTCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCa                                                                      >  2:3911897/1‑93 (MQ=255)

                   aTCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCa                                                                      >  1:1223119/1‑93 (MQ=255)

                    tCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGGCATTACGATTGTGGTGGTGGTGATGCTGATCGCAc                                                                     <  2:4534793/93‑1 (MQ=255)

                     cGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACt                                                                    <  1:4922524/93‑1 (MQ=255)

                       tCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGTGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTaa                                                                  >  2:7689369/1‑93 (MQ=255)

                       tCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGGACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGTTGATGGCACTac                                                                  >  1:5545409/1‑92 (MQ=255)

                        ccTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAAc                                                                 <  1:1136724/93‑1 (MQ=255)

                        ccTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAAc                                                                 <  2:5732332/93‑1 (MQ=255)

                         cTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTTCGATTGTGGTGGTGGTGATGCTGATCGCACTAAcc                                                                <  2:5071346/93‑1 (MQ=255)

                         cTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAAcc                                                                >  2:1560998/1‑93 (MQ=255)

                          tCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTTGCGATGCTGATCGCCCTAAccc                                                               >  2:4236317/1‑92 (MQ=255)

                          tCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACcg                                                               >  2:7399207/1‑93 (MQ=255)

                           cAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACcgc                                                              >  1:2104575/1‑93 (MQ=255)

                             aCCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCgcgc                                                            <  2:3329937/93‑1 (MQ=255)

                             aCCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCgcgc                                                            >  2:4966997/1‑93 (MQ=255)

                             aCCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCgcgc                                                            <  1:6259000/93‑1 (MQ=255)

                              cccGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCt                                                           >  1:2671803/1‑93 (MQ=255)

                              cccGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGGGATGCTGATCGCACTAACCGCGCt                                                           >  2:4667251/1‑93 (MQ=255)

                               ccGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTg                                                          <  2:4342311/93‑1 (MQ=255)

                                cGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTgg                                                         <  2:1714381/93‑1 (MQ=255)

                                cGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTgg                                                         >  2:590879/1‑93 (MQ=255)

                                 gCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGt                                                        <  2:5771492/93‑1 (MQ=255)

                                  cTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTc                                                       <  1:443107/93‑1 (MQ=255)

                                  cTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTc                                                       <  1:6135293/93‑1 (MQ=255)

                                    cGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCgg                                                     <  2:5223045/93‑1 (MQ=255)

                                    cGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCgg                                                     >  1:2031392/1‑93 (MQ=255)

                                     gTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGc                                                    <  1:4839703/93‑1 (MQ=255)

                                     gTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGc                                                    <  1:7650733/93‑1 (MQ=255)

                                     gTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGAGATGCTGATCGCACTAACCGCGCTGGTCGGc                                                    >  1:6406832/1‑93 (MQ=255)

                                      tCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCg                                                   >  1:4156614/1‑93 (MQ=255)

                                        cccTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTa                                                 <  1:7125760/92‑1 (MQ=255)

                                        aCCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTa                                                 >  2:3640177/1‑93 (MQ=255)

                                        aCCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTa                                                 >  1:3480900/1‑93 (MQ=255)

                                         ccTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTAc                                                <  1:1661148/93‑1 (MQ=255)

                                          cTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACt                                               >  1:7899719/1‑93 (MQ=255)

                                           tGGTTTAGTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTg                                              <  1:5149796/93‑1 (MQ=255)

                                            ggTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTgg                                             <  2:2928583/93‑1 (MQ=255)

                                            ggTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTgg                                             >  2:1132662/1‑93 (MQ=255)

                                              tttATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTggcg                                           <  2:6790524/93‑1 (MQ=255)

                                              tttATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTggcg                                           >  1:4470366/1‑93 (MQ=255)

                                                  ttCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCAt                                       <  2:5125000/93‑1 (MQ=255)

                                                  ttCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCAt                                       <  2:3436040/93‑1 (MQ=255)

                                                  ttCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCAt                                       <  1:3832964/93‑1 (MQ=255)

                                                    cgtcgGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCg                                     <  2:2574207/93‑1 (MQ=255)

                                                     gtcgGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGt                                    <  1:4522397/93‑1 (MQ=255)

                                                     gtcgGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGt                                    >  1:1380109/1‑93 (MQ=255)

                                                      tcgGGGAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGtt                                   <  1:5714199/93‑1 (MQ=255)

                                                        ggggAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTa                                 >  1:519687/1‑93 (MQ=255)

                                                        ggggAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTa                                 <  2:2718451/93‑1 (MQ=255)

                                                        ggggAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTa                                 >  2:3181007/1‑93 (MQ=255)

                                                          ggtgTACAACGCCCCGGTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGGACTGGCGGCATCGTTTAAc                               >  2:7025540/4‑93 (MQ=255)

                                                          ggAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAAc                               >  1:5088986/1‑93 (MQ=255)

                                                           gAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACg                              >  2:6424935/1‑93 (MQ=255)

                                                           gAGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACg                              <  2:1889959/93‑1 (MQ=255)

                                                            aGTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGa                             >  2:7336894/1‑93 (MQ=255)

                                                             gTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGaa                            >  2:2232339/1‑93 (MQ=255)

                                                             gTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGaa                            >  1:3579343/1‑93 (MQ=255)

                                                             gTACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGaa                            >  2:3950887/1‑93 (MQ=255)

                                                              tACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAAt                           <  2:6114921/93‑1 (MQ=255)

                                                              tACAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAAt                           <  1:2968819/93‑1 (MQ=255)

                                                               aCAACGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAAtt                          <  1:2468660/93‑1 (MQ=255)

                                                                 aaCGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTa                        >  2:7315967/1‑93 (MQ=255)

                                                                 aaCACCCCGTTGCCATTACTATGGTTGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTa                        <  2:6322204/93‑1 (MQ=255)

                                                                   cGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTAtc                      <  1:2358545/93‑1 (MQ=255)

                                                                   cGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTAtc                      <  1:2832147/93‑1 (MQ=255)

                                                                   cGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTAtc                      <  1:1617395/93‑1 (MQ=255)

                                                                   cGCCCCGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTAtc                      <  2:5508309/93‑1 (MQ=255)

                                                                     ccccGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTAtctc                    >  1:6039810/1‑93 (MQ=255)

                                                                      cccGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATctct                   <  1:4246177/93‑1 (MQ=255)

                                                                       ccGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTa                  <  1:5816218/93‑1 (MQ=255)

                                                                       ccGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTa                  >  1:4386388/1‑93 (MQ=255)

                                                                       ccGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTa                  <  2:2170233/93‑1 (MQ=255)

                                                                       ccGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTa                  <  1:1809177/93‑1 (MQ=255)

                                                                        cGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTAt                 >  2:1775328/1‑93 (MQ=255)

                                                                        cGTTGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTAt                 <  2:7467449/93‑1 (MQ=255)

                                                                          ttGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTAtgc               >  1:4642565/1‑93 (MQ=255)

                                                                           tGCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTAtgct              >  1:6862904/1‑93 (MQ=255)

                                                                            gCCATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTAtgctg             >  1:3194727/1‑93 (MQ=255)

                                                                              cATTACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTATGCTGcc           <  1:4560339/93‑1 (MQ=255)

                                                                                ttACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTATGCTGCCGa         >  1:4009865/1‑93 (MQ=255)

                                                                                ttACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTATGCTGCCGa         >  2:2415145/1‑93 (MQ=255)

                                                                                ttACGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGCTTAACGCATTTATCTCTATGCTGCCGa         >  1:5902051/1‑93 (MQ=255)

                                                                                   cGATTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTATGCTGCCGAAGt      >  2:5414512/1‑93 (MQ=255)

                                                                                     aTTGTGGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTATGCTGCCGAAGttt    <  2:3558237/93‑1 (MQ=255)

                                                                                       tgtgGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTATGCTGCCGAAGTTtaa  <  1:2245247/93‑1 (MQ=255)

ACTGGCGCTGTTTTTTGCCATCGTCCTCAACCCGCTCGTCACCTGGTTTATTCGTCGGGGAGTACAACGCCCCGTTGCCATTACGATTGTAGTGGTGGTGATGCTGATCGCACTAACCGCGCTGGTCGGCGTACTGGCGGCATCGTTTAACGAATTTATCTCTATGCTGCCGAAGTTTAA  >  NC_000913/3609609‑3609788

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