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.2%	735.5 / ‑1.3	235	coding (302/408 nt)	phnQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (1/1); new base C (132/101); total (133/102)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 257 total aligned reads displayed.

GAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCG  >  NC_000913/4311403‑4311582

gAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑a                                                                                          <  1:3996972/93‑1 (MQ=255)

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

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

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

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

     aTGATTAATGTGGATAGAGTTTCTTTTTGGGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑ccgtt                                                                                    <  1:951649/93‑5 (MQ=255)

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

       gATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTcc                                                                                   >  2:2225891/1‑93 (MQ=255)

          tAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGc                                                                                >  2:2911789/1‑93 (MQ=255)

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

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

             tgtgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCTGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCa                                                                             >  1:2384524/1‑93 (MQ=255)

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

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

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

               tgGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAgg                                                                           >  2:2408302/1‑93 (MQ=255)

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

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

                 gATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGcc                                                                         <  1:4093449/93‑1 (MQ=255)

                  aTAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCg                                                                        >  1:1850322/1‑93 (MQ=255)

                   tAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGa                                                                       >  1:333489/1‑93 (MQ=255)

                    tgagTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCCCCCTGCTCTCGCTGTCC‑ACCGTTCCGGCCCAGGCCGat                                                                      <  2:2619374/92‑1 (MQ=255)

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

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

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

                       gTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGAtatt                                                                   <  1:1047276/93‑1 (MQ=255)

                       gTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGAtatt                                                                   <  1:975463/93‑1 (MQ=255)

                       gTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGAtatt                                                                   >  2:4593036/1‑93 (MQ=255)

                       gTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGAtatt                                                                   >  2:2762445/1‑93 (MQ=255)

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

                         ttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGAtattat                                                                 <  1:3698955/93‑1 (MQ=255)

                          tCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATc                                                                >  1:169253/1‑93 (MQ=255)

                          tCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATc                                                                >  2:1112909/1‑93 (MQ=255)

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

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

                            tttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATcga                                                              <  1:1645132/93‑1 (MQ=255)

                                 gAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGaca                                                         <  1:4429820/93‑1 (MQ=255)

                                      tGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCGC‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCAtt                                                    <  1:3183802/93‑1 (MQ=255)

                                      tGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCAtt                                                    <  1:2259296/93‑1 (MQ=255)

                                          tAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCa                                                >  2:3237223/1‑93 (MQ=255)

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

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

                                            aTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTCGCAat                                              >  1:2474090/1‑93 (MQ=255)

                                             tGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAata                                             >  2:40971/1‑93 (MQ=255)

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

                                                aaCGCTTTCCTCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTc                                          <  1:3876528/93‑1 (MQ=255)

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

                                                  cGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcag                                        <  2:1165254/93‑1 (MQ=255)

                                                  cGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcag                                        <  1:3008110/93‑1 (MQ=255)

                                                   gCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagc                                       <  1:683541/93‑1 (MQ=255)

                                                   gCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagc                                       >  2:2286353/1‑93 (MQ=255)

                                                   gCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagc                                       <  1:835974/93‑1 (MQ=255)

                                                   gCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagc                                       <  1:2842582/93‑1 (MQ=255)

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

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

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

                                                     tttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagcag                                     <  2:3544207/93‑1 (MQ=255)

                                                      ttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCCGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAgg                                    >  1:2953933/1‑93 (MQ=255)

                                                       tCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGc                                   <  2:923203/93‑1 (MQ=255)

                                                        ccGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCg                                  >  1:1473778/1‑93 (MQ=255)

                                                         cGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGa                                 >  1:462744/1‑93 (MQ=255)

                                                          gCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAt                                <  1:2171455/93‑1 (MQ=255)

                                                          gCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGACATTATCGACGACACCATTGGCAATATTCAGCAGGCGAt                                >  2:3392351/1‑93 (MQ=255)

                                                           cTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCTGGCGAtt                               >  1:3522561/1‑93 (MQ=255)

                                                            tATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTa                              >  2:3195224/1‑93 (MQ=255)

                                                             aTTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTaa                             <  2:2152553/93‑1 (MQ=255)

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

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

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

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

                                                               ttCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCAATTAATg                           >  2:2353925/1‑93 (MQ=255)

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

                                                                  ttttATTCCCCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGCCAATATTCAGCAGGCGATTAATGACg                        <  1:2610791/93‑1 (MQ=255)

                                                                  ttttATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACg                        >  2:1727345/1‑93 (MQ=255)

                                                                  ttttATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACg                        >  2:4147398/1‑93 (MQ=255)

                                                                   tttATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGGCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGc                       >  1:1328199/1‑93 (MQ=255)

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

                                                                    ttATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCt                      >  1:4026827/1‑93 (MQ=255)

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

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

                                                                        tCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTtata                  <  1:596755/93‑1 (MQ=255)

                                                                        tCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTtata                  >  2:3040136/1‑93 (MQ=255)

                                                                         cccccTGCTCACGCTGTCC‑ACCGTGCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATaa                 <  2:4407086/91‑1 (MQ=255)

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

                                                                          aCCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAc                <  1:370277/93‑1 (MQ=255)

                                                                           cccTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAcc               <  2:2053948/93‑1 (MQ=255)

                                                                             cTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCt             >  1:2885529/1‑93 (MQ=255)

                                                                                cTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc          <  2:2567069/93‑1 (MQ=255)

                                                                                cTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc          <  1:3509624/93‑1 (MQ=255)

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

                                                                                cTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAc          <  2:2231984/93‑1 (MQ=255)

                                                                                   aCGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGt       >  2:3748741/1‑93 (MQ=255)

                                                                                   aCGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGt       >  1:2079761/1‑93 (MQ=255)

                                                                                    cGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCTATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTg      <  1:3296593/93‑1 (MQ=255)

                                                                                     gCTGTCC‑ACCGTTCCGGCCCAGGCCGATATCATCGACGACACGATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTgg     <  2:4537539/93‑1 (MQ=255)

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

                                                                                       tGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGc   <  1:638040/93‑1 (MQ=255)

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

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

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

GAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCG  >  NC_000913/4311403‑4311582

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