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.3%	1325.1 / 0.5	420	coding (302/408 nt)	phnQ	hypothetical protein
Reads supporting (aligned to +/- strand): ref base . (2/1); new base C (238/179); total (240/180)
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.30e-01

Only 100 of 450 total aligned reads displayed.

GAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑A‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCG  >  NC_000913/4311403‑4311582

gAGTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑a                                                                                          >  2:7597763/1‑93 (MQ=255)

  gTAATGATTAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGACACGCTTTCCGCTATTTATTTTATTCACCCTGCTCACGCTGTCC‑A‑cc                                                                                       >  1:7599334/1‑92 (MQ=255)

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

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

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

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

         ttAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCgg                                                                                 >  2:790578/1‑93 (MQ=255)

         ttAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCgg                                                                                 >  1:1807917/1‑93 (MQ=255)

         ttAATGTGGATAGAGTTTCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCgg                                                                                 >  2:3192450/1‑93 (MQ=255)

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

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

                        tttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGAtatta                                                                  >  1:3646965/1‑93 (MQ=255)

                        tttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGAtatta                                                                  >  2:1394667/1‑93 (MQ=255)

                        tttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGAtatta                                                                  >  2:2118758/1‑93 (MQ=255)

                        tttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGAtatta                                                                  >  2:5435657/1‑93 (MQ=255)

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

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

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

                         ttCTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCCATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGAtattat                                                                 >  2:6278686/1‑93 (MQ=255)

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

                           cTTTTTGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATcg                                                               >  2:3325039/1‑93 (MQ=255)

                             ttttGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATcgac                                                             >  1:4356189/1‑93 (MQ=255)

                                tGAGGTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCgacgac                                                          <  1:5821707/93‑1 (MQ=255)

                                   ggTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACAcc                                                       <  2:7000934/93‑1 (MQ=255)

                                    gTTGGCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCa                                                      >  1:2890233/1‑93 (MQ=255)

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

                                       ggCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTg                                                   >  2:7743737/1‑93 (MQ=255)

                                       gacTAATGAAGCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTg                                                   <  2:150860/91‑1 (MQ=255)

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

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

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

                                        gCTAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTgg                                                  >  1:6288553/1‑93 (MQ=255)

                                          tAATGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCa                                                <  1:4545635/93‑1 (MQ=255)

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

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

                                           aaTGAAACGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGACGATATTATCGACGACACCATTGGCaa                                               >  1:5328722/1‑93 (MQ=255)

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

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

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

                                               aaaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATAtt                                           >  1:6518414/1‑93 (MQ=255)

                                               aaaCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATAtt                                           >  1:7414166/1‑93 (MQ=255)

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

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

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

                                                 aCGCTTTCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCCGC‑CGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTca                                         >  1:1465860/1‑93 (MQ=255)

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

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

                                                    cTTTCCGCTATTTCTTTTATTCACCCTTCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTcagca                                      <  2:604974/93‑1 (MQ=255)

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

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

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

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

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

                                                      ttCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAgg                                    <  2:2901536/93‑1 (MQ=255)

                                                       tCCGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGc                                   >  2:149510/1‑93 (MQ=255)

                                                         cGCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCgt                                 >  2:2260693/1‑92 (MQ=255)

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

                                                          gCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAt                                >  1:2668682/1‑93 (MQ=255)

                                                          gCTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAt                                >  1:548175/1‑93 (MQ=255)

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

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

                                                           cTATTTCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGAtt                               >  1:7398201/1‑93 (MQ=255)

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

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

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

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

                                                               ttCTTTTATTCTCCCTGCTCACGCTGTCC‑ACCGTACCGGCCCAGGCCGATATTATAGACGACACCATTGGCAATATTCAGCAGGCGATTAATg                           <  1:5941956/93‑1 (MQ=255)

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

                                                               ttCTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATg                           >  2:204763/1‑93 (MQ=255)

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

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

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

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

                                                                 cTTTTATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGAc                         <  1:244565/93‑1 (MQ=255)

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

                                                                   tttATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGCCGACACCATTGGCAATATTCAGCAGGCGATTAATGACGc                       >  1:1991231/1‑93 (MQ=255)

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

                                                                     tATTCACCTTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCtt                     >  1:6222913/1‑93 (MQ=255)

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

                                                                     tATTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCtt                     <  2:3952437/93‑1 (MQ=255)

                                                                      aTTCACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTta                    >  1:5102259/1‑93 (MQ=255)

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

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

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

                                                                         cACCCTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATaa                 >  2:559692/1‑93 (MQ=255)

                                                                           cccTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAAcc               >  1:6995250/1‑93 (MQ=255)

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

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

                                                                             cTGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCt             <  1:3011929/93‑1 (MQ=255)

                                                                              tGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTg            >  2:472270/1‑93 (MQ=255)

                                                                              tGCTCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGAGATTAATGACGCTTATAACCCTg            >  2:4283954/1‑93 (MQ=255)

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

                                                                                 tCACGCTGTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGAcc         >  1:4782695/1‑93 (MQ=255)

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

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

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

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

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

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

                                                                                        gTCC‑ACCGTTCCGGCCCAGGCCGATATTATCGACGACACCATTGGCAATATTCAGCAGGCGATTAATGACGCTTATAACCCTGACCGTGGGCg  >  2:2618683/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` `< 24 ≤` `ATCG/ATCG` `< 30 ≤` `ATCG/ATCG` `< 35 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`