BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,311,496	+C	coding (301/408 nt)	phnQ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,311,495	1	.	C	98.2%	544.8 / 1.6	170	coding (302/408 nt)	phnQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/2); new base C (76/91); total (77/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 193 total aligned reads displayed.

AGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCG  >  NC_000913/4311404‑4311582

aGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑c                                                                                        <  1:402603/93‑1 (MQ=255)

aGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑c                                                                                        >  2:3140414/1‑93 (MQ=255)

 gTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑cc                                                                                       <  2:621475/93‑2 (MQ=255)

 gTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑cc                                                                                       >  2:3831820/1‑92 (MQ=255)

  tAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑ccg                                                                                      >  1:3407104/1‑91 (MQ=255)

   aaTGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑ccgt                                                                                     >  1:553796/1‑90 (MQ=255)

   aaTGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑ccgt                                                                                     >  1:1488048/1‑90 (MQ=255)

    aTGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑ccgtt                                                                                    >  1:277867/1‑89 (MQ=255)

      gATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTcc                                                                                   <  1:844881/93‑1 (MQ=255)

       aTTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCg                                                                                  >  2:3500519/1‑93 (MQ=255)

       aTTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCg                                                                                  <  1:2214778/93‑1 (MQ=255)

         tAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGc                                                                                <  1:2866775/93‑1 (MQ=255)

          aaTGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGcc                                                                               <  1:1726579/93‑1 (MQ=255)

           aTGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGccc                                                                              <  1:1558411/93‑1 (MQ=255)

           aTGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGccc                                                                              >  2:418836/1‑93 (MQ=255)

           aTGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGccc                                                                              <  1:2494195/93‑1 (MQ=255)

           aTGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGccc                                                                              <  1:4449511/93‑1 (MQ=255)

            tgtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCa                                                                             <  2:1119824/93‑1 (MQ=255)

            tgtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCa                                                                             >  1:1250878/1‑93 (MQ=255)

            tgtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCa                                                                             >  2:4340931/1‑93 (MQ=255)

            tgtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCa                                                                             <  1:4349696/93‑1 (MQ=255)

             gtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAg                                                                            <  1:2276072/93‑1 (MQ=255)

               ggATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCGCCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGc                                                                          <  2:2175632/93‑1 (MQ=255)

               ggATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGc                                                                          >  2:512730/1‑93 (MQ=255)

               ggATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGc                                                                          >  2:2731386/1‑93 (MQ=255)

                gATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGcc                                                                         >  2:2074409/1‑93 (MQ=255)

                  tAGAGTTTCTTTTTGAGGTTGGCTAATGATACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGa                                                                       >  2:4608525/1‑93 (MQ=255)

                   agagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGat                                                                      >  2:335806/1‑93 (MQ=255)

                    gagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGata                                                                     >  1:1826464/1‑93 (MQ=255)

                    gagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGata                                                                     <  1:911706/93‑1 (MQ=255)

                     agTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTAACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGatat                                                                    >  1:3191471/1‑93 (MQ=255)

                        ttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCCCGCTGTCC‑ACCGTTCCGGCCCAGGCCGAtattat                                                                 <  2:1786231/93‑1 (MQ=255)

                        ttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGAtattat                                                                 >  1:2750056/1‑93 (MQ=255)

                        ttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGAtattat                                                                 >  1:2548294/1‑93 (MQ=255)

                         tCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATc                                                                <  1:3581269/93‑1 (MQ=255)

                         tCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATc                                                                <  2:4170868/93‑1 (MQ=255)

                         tCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATc                                                                <  2:2060091/93‑1 (MQ=255)

                         tCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATc                                                                <  1:4737391/93‑1 (MQ=255)

                           tttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATcga                                                              >  1:3294151/1‑93 (MQ=255)

                            ttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATcgac                                                             <  2:774384/93‑1 (MQ=255)

                              ttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATcgacga                                                           <  1:3279280/93‑1 (MQ=255)

                                 aGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGacac                                                        <  1:1022273/93‑1 (MQ=255)

                                 aGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGacac                                                        >  1:1419361/1‑93 (MQ=255)

                                 aGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGACGATATTATCGACGacac                                                        >  2:265673/1‑93 (MQ=255)

                                    ttGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCAt                                                     <  2:1931038/93‑1 (MQ=255)

                                    ttGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCAt                                                     >  1:4720110/1‑93 (MQ=255)

                                     tGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑C‑CGTTTCGGCCCAGGGCGATATTTATGACGACACCATTg                                                   >  2:217464/1‑93 (MQ=255)

                                      ggCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTg                                                   <  2:2605126/93‑1 (MQ=255)

                                      ggCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTg                                                   <  1:4540856/93‑1 (MQ=255)

                                      ggCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTg                                                   <  2:4122479/93‑1 (MQ=255)

                                       gCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTgg                                                  <  2:3349959/93‑1 (MQ=255)

                                       gCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTgg                                                  <  1:2021053/93‑1 (MQ=255)

                                       gCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTgg                                                  <  2:128285/93‑1 (MQ=255)

                                        cTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGc                                                 <  2:3701360/93‑1 (MQ=255)

                                          aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCaa                                               >  2:2358571/1‑93 (MQ=255)

                                          aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCaa                                               >  2:3743312/1‑93 (MQ=255)

                                           aTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAat                                              >  1:3718038/1‑93 (MQ=255)

                                           aTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAat                                              >  1:1506077/1‑93 (MQ=255)

                                             gAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAatat                                            <  1:2171093/93‑1 (MQ=255)

                                             gAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAatat                                            >  1:1906800/1‑93 (MQ=255)

                                              aaaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATAtt                                           <  2:4261508/93‑1 (MQ=255)

                                               aaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTc                                          <  2:3126454/93‑1 (MQ=255)

                                               aaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTc                                          >  2:1398525/1‑93 (MQ=255)

                                               aaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTc                                          >  2:4038801/1‑93 (MQ=255)

                                                aCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTca                                         <  2:1730859/93‑1 (MQ=255)

                                                   cTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagca                                      <  1:67164/93‑1 (MQ=255)

                                                    tttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagcag                                     >  2:932383/1‑93 (MQ=255)

                                                            aTTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTaa                             >  1:3284896/1‑93 (MQ=255)

                                                             tttCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAAt                            >  1:1071509/1‑93 (MQ=255)

                                                             tttCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAAt                            <  1:3219339/93‑1 (MQ=255)

                                                              ttCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATg                           <  2:1820513/93‑1 (MQ=255)

                                                              ttCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATg                           <  1:4092154/93‑1 (MQ=255)

                                                               tCTTTTATTCGCCCTGCTCACGCTGTCCTC‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGa                          <  1:1146286/93‑1 (MQ=255)

                                                               tCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGa                          <  2:2953553/93‑1 (MQ=255)

                                                               tCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGa                          >  1:821756/1‑93 (MQ=255)

                                                                tttttATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGAc                         <  1:214669/92‑1 (MQ=255)

                                                                  tttATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGc                       <  2:1555707/93‑1 (MQ=255)

                                                                   ttATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCt                      <  2:2876005/93‑1 (MQ=255)

                                                                    tATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCtt                     >  1:298801/1‑93 (MQ=255)

                                                                     aTTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTta                    <  2:166982/93‑1 (MQ=255)

                                                                     aTTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGCCGATTAATGACGCTta                    <  2:4508102/93‑1 (MQ=255)

                                                                      ttCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTtat                   <  2:2133128/93‑1 (MQ=255)

                                                                      ttCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTtat                   >  1:3843900/1‑93 (MQ=255)

                                                                        cACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATaa                 <  2:3064627/93‑1 (MQ=255)

                                                                         ccccTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAc                <  2:1238489/92‑1 (MQ=255)

                                                                         aCCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAc                >  2:4050024/1‑93 (MQ=255)

                                                                           ccTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAccc              >  1:4569941/1‑93 (MQ=255)

                                                                              gCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGa           <  2:59846/93‑1 (MQ=255)

                                                                               cTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc          >  1:1361259/1‑93 (MQ=255)

                                                                                tCACGCTGTCC‑ACCGTTCCGGCCCAGGGCGATATTATCGACGACACCATTGGCAACATTCAGCAGGCGATTAATGACGCTTATAACACTGAcc         >  2:1395625/1‑93 (MQ=255)

                                                                                tCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAcc         <  2:4448616/93‑1 (MQ=255)

                                                                                tCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCAGAcc         >  2:576143/1‑93 (MQ=255)

                                                                                   cGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTg      <  1:1101254/93‑1 (MQ=255)

                                                                                   cGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTg      >  1:2846567/1‑93 (MQ=255)

                                                                                   cGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTg      <  2:1302770/93‑1 (MQ=255)

                                                                                    gCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTgg     >  2:1142604/1‑93 (MQ=255)

                                                                                     cTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTggg    <  1:235079/93‑1 (MQ=255)

                                                                                      tGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGc   >  2:3112389/1‑93 (MQ=255)

                                                                                      tGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGc   >  2:85243/1‑93 (MQ=255)

                                                                                       gTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCg  >  2:2826097/1‑93 (MQ=255)

AGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCG  >  NC_000913/4311404‑4311582

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