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	99.6%	819.4 / ‑4.0	235	coding (302/408 nt)	phnQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (0/1); new base C (126/108); total (126/109)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 254 total aligned reads displayed.

GAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCA‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCG  >  NC_000913/4311403‑4311582

gggTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCa                                                                                          >  1:5725564/3‑93 (MQ=255)

gAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCa                                                                                          <  1:3320176/93‑1 (MQ=255)

   tAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCA‑ccg                                                                                      >  1:5844484/1‑91 (MQ=255)

   tAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCA‑ccg                                                                                      >  1:629552/1‑91 (MQ=255)

    aaTGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCA‑ccgt                                                                                     <  2:6000899/93‑4 (MQ=255)

       gATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTcc                                                                                   >  2:5230834/1‑93 (MQ=255)

       gATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTcc                                                                                   <  1:6272786/93‑1 (MQ=255)

       gATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTcc                                                                                   <  1:1044572/93‑1 (MQ=255)

        aTTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCg                                                                                  <  1:4211145/93‑1 (MQ=255)

        aTTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCg                                                                                  >  2:3536429/1‑93 (MQ=255)

            aTGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGccc                                                                              <  1:1704078/93‑1 (MQ=255)

             tgtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCa                                                                             <  1:3524309/93‑1 (MQ=255)

               tgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAgg                                                                           >  1:5982219/1‑93 (MQ=255)

               tgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAgg                                                                           >  2:787099/1‑93 (MQ=255)

                ggATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGc                                                                          >  2:4416073/1‑93 (MQ=255)

                ggATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCACTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGc                                                                          >  2:4527222/1‑93 (MQ=255)

                  aTAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCg                                                                        <  1:563374/93‑1 (MQ=255)

                   tAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGa                                                                       <  2:1399528/93‑1 (MQ=255)

                   tAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGa                                                                       >  2:2280278/1‑93 (MQ=255)

                   tAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGa                                                                       >  2:5656855/1‑93 (MQ=255)

                   tAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGa                                                                       >  1:6033366/1‑93 (MQ=255)

                    agagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGat                                                                      <  1:3401391/93‑1 (MQ=255)

                     gggTTTCTTTTTGAGGTTGGCTAATGACACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGata                                                                     >  1:6175621/3‑93 (MQ=255)

                     gagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGata                                                                     >  2:852673/1‑93 (MQ=255)

                     gagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGata                                                                     <  1:4397103/93‑1 (MQ=255)

                      agTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGatat                                                                    >  2:5292702/1‑93 (MQ=255)

                      agTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGatat                                                                    >  2:583724/1‑93 (MQ=255)

                       gTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGAtatt                                                                   >  2:204563/1‑93 (MQ=255)

                       gTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGAtatt                                                                   >  2:5136815/1‑93 (MQ=255)

                       gTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGAtatt                                                                   >  2:2129600/1‑93 (MQ=255)

                        tttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGAtatta                                                                  >  2:473459/1‑93 (MQ=255)

                         ttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGAtattat                                                                 >  1:4038014/1‑93 (MQ=255)

                         ttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGAtattat                                                                 >  1:1870834/1‑93 (MQ=255)

                             ttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATcgac                                                             <  2:3148699/93‑1 (MQ=255)

                             ttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATcgac                                                             <  1:2428378/93‑1 (MQ=255)

                             ttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATcgac                                                             <  1:1335947/93‑1 (MQ=255)

                              tttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATcgacg                                                            <  2:6193377/93‑1 (MQ=255)

                                   ggTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACAcc                                                       >  1:469833/1‑93 (MQ=255)

                                      tGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCAtt                                                    <  2:421225/93‑1 (MQ=255)

                                      tGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCAtt                                                    <  2:5007394/93‑1 (MQ=255)

                                      tGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCCATATTATCGACGACACCAtt                                                    <  1:6219486/93‑1 (MQ=255)

                                       ggCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTg                                                   >  1:1842950/1‑93 (MQ=255)

                                       ggCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTg                                                   <  2:1289781/93‑1 (MQ=255)

                                        gCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTgg                                                  <  1:1435293/93‑1 (MQ=255)

                                           aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCaa                                               >  2:281608/1‑93 (MQ=255)

                                           aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCaa                                               >  1:2623644/1‑93 (MQ=255)

                                            aTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAat                                              <  2:3034391/93‑1 (MQ=255)

                                             tGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAata                                             >  2:4295405/1‑93 (MQ=255)

                                               aaaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATAtt                                           >  1:1949044/1‑93 (MQ=255)

                                               aaaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATAtt                                           >  2:5940716/1‑93 (MQ=255)

                                               aaaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATAtt                                           >  1:907894/1‑93 (MQ=255)

                                                aaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTc                                          >  2:6073692/1‑93 (MQ=255)

                                                  cGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcag                                        >  1:3465113/1‑93 (MQ=255)

                                                   gCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagc                                       <  1:2684675/93‑1 (MQ=255)

                                                   gCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagc                                       >  1:2238197/1‑93 (MQ=255)

                                                     tttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagcag                                     >  1:1112549/1‑93 (MQ=255)

                                                     tttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagcag                                     <  1:2837899/93‑1 (MQ=255)

                                                       tCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGc                                   <  1:2062925/93‑1 (MQ=255)

                                                       tCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGc                                   <  1:768117/93‑1 (MQ=255)

                                                          gCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAt                                <  2:2733702/93‑1 (MQ=255)

                                                          gCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAt                                >  1:2000185/1‑93 (MQ=255)

                                                          ccTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAt                                <  1:2602587/92‑1 (MQ=255)

                                                             aTTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTaa                             <  2:2464364/93‑1 (MQ=255)

                                                             aTTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTaa                             >  2:1516434/1‑93 (MQ=255)

                                                             aTTTCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTaa                             >  1:3880410/1‑93 (MQ=255)

                                                              tttCTTTTATTCACCCCGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAAt                            <  2:5735549/93‑1 (MQ=255)

                                                                tCTTTTATTCCCCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGa                          <  2:5564563/93‑1 (MQ=255)

                                                                tCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGa                          <  2:2008904/93‑1 (MQ=255)

                                                                tCTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGa                          <  1:5024281/93‑1 (MQ=255)

                                                                 cTTTTATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGAc                         <  1:3805007/93‑1 (MQ=255)

                                                                   tttATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGc                       <  2:176110/93‑1 (MQ=255)

                                                                   tttATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGc                       <  1:4533902/93‑1 (MQ=255)

                                                                    ttATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTACTGACGCt                      <  1:6359073/93‑1 (MQ=255)

                                                                     tATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCtt                     >  1:1579154/1‑93 (MQ=255)

                                                                     tATTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCtt                     >  1:4453370/1‑93 (MQ=255)

                                                                      aTTCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTta                    >  2:5318276/1‑93 (MQ=255)

                                                                       ttCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTtat                   <  2:668444/93‑1 (MQ=255)

                                                                       ttCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTtat                   >  1:3045558/1‑93 (MQ=255)

                                                                        tCACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTtata                  >  2:4669700/1‑93 (MQ=255)

                                                                         cACCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATaa                 <  1:5592821/93‑1 (MQ=255)

                                                                          aCCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAc                >  2:1865232/1‑93 (MQ=255)

                                                                          aCCCTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAc                >  2:1595422/1‑93 (MQ=255)

                                                                           cccTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAcc               >  2:3924162/1‑93 (MQ=255)

                                                                           cccTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAcc               >  2:3491102/1‑93 (MQ=255)

                                                                            ccTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAccc              >  1:4142034/1‑93 (MQ=255)

                                                                            ccTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAccc              >  2:3713882/1‑93 (MQ=255)

                                                                             cTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCt             >  1:4895843/1‑93 (MQ=255)

                                                                             cTGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCt             <  2:6311276/93‑1 (MQ=255)

                                                                              tGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTg            <  1:691006/93‑1 (MQ=255)

                                                                              tGCTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTg            <  1:3671577/93‑1 (MQ=255)

                                                                                cTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc          <  2:3731042/93‑1 (MQ=255)

                                                                                cTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc          <  1:5519449/93‑1 (MQ=255)

                                                                                cTCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc          <  1:309909/93‑1 (MQ=255)

                                                                                 tCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAcc         <  2:2554007/93‑1 (MQ=255)

                                                                                 tCACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAcc         >  1:9964/1‑93 (MQ=255)

                                                                                 tCACCCTGTCCCC‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAcc         <  2:2572299/93‑1 (MQ=255)

                                                                                  cACGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCg        >  1:5301653/1‑93 (MQ=255)

                                                                                   aCGCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGt       >  2:906596/1‑93 (MQ=255)

                                                                                     gCTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTgg     >  1:3957897/1‑93 (MQ=255)

                                                                                      cTGT‑CCACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTggg    >  1:5176854/1‑93 (MQ=255)

GAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGT‑CCA‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCG  >  NC_000913/4311403‑4311582

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